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NUMERICAL PROCEDURES -EOR TEE CALCULATION OP THE 


^ r . .. ..STRESSES ..IN. MONOCO^UES - ;• . - . 

Ill - CALCULATION OP THE .BENDING MOMENTS IN PU SB LAGS 'PRitME S 
By N. J. Hoff, Paul A. Libby, and Bertram Klein 


■ . .SUMMARY 


This repo-r.t deals with the -carculat i on- "o f the bending 
moments in and the distortions of fuselage rings upon which 
known concentrated, and distributed 1-oadB are acting. In the 
procedure suggested, the ring i.s divided into -a. number of 
beams each, having- a constant -radius of curvature",- -The forces 
and moments caused in the end sections of the beams by indi- 
vidual unit displacements of. t-h.e. ,en.d- -sect ions. ar-er r li'st ed in 
a table designated as. the oper at i t abl e in-;. conformity ■ with 
Southwell 1 s nomenclature-. • . - ; . ...* '■ *. - -- i 

• • . ■_ '■ . t “* 

The operations table and; the.- e:Rtte.fn-al.: loads' are “equiva- 
lent to a set of linear equations, For. their, solution'- the' 
following three procedures are presented: 


1. Southwell's method o f sy g-tbm&t I c relaxations. This 
is a step-by-step approximation procedure guided by the phye^ 
ical interpretation of the changes in the values of the uni—'"' 
known s- - . .. r . ; -> > - ».~r — - ,-r 


• ■ ; ; j >- ' I *. ■ - 

2. The ; .g^ pwing ;up-it procedure in which the individual 
beams ar.e --c onMnad gu c. he ^ s^ve ly ~,i.nt o beams of- increasing 
length , 'unti 1 ' finally .the . fj,n.hi.r e. ring -becomes .a single be 

T n 1 on rtVi "U « a _ — _ j.i _ -C* i. 'IZ. T 


am , 
ihr e e 


m . . > — : - T yt**. »• WWW w*i*v s . o. D x 

In ’each ,.s.t ep ^of- ..the -pr oc,edure. a? ...s.et' -.o f not, more t han 
simultaneous linear .aq.ua,ti,Qns is solved, 

3. .Solution of t-he e.ntir.ei , set of simultaneous equations' 
by the methods of ,the matrix, calculus, ■ ■ : 
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In order to demonstrate the manner in which the calcula- 
tions may be carried out, the following numerical examples 
are worked out: ' V" .. ; T I 


1. Curved beam with booth its end sections rigidly fixed. 

The load i s'-a concentrated . force . 

2. Egg-shape ring with" symmetric concentrated loads 

3. Circular ' rin-g -with antisymmetric c’once nitrated loads 

and shear flow (torsion of the fuselage) 

4. Same with V-braces incorporated in the ring 

5. Egg-shape ring with antisymmetric concentrated loads 

and shear flow (torsion of the fuselage) 

6. Same with - ; V-braces incorporated in the ring ■: • 

' ‘ . * . ' ' ' , > 

The results of these calculation's are checked, whenever 
possible, by calculations .carried out according to known 
methods of analysis... .The' agreement is found to be good. 

The amount o f ■ wor k- n e ce s sar y « for the solution of 'ring 
problems by the me th od g de scr i bed • i n> the-presqnt report ib' ‘ 
practically independent of the degree of -redundancy of the ” " 
structure. For this reason the methods are recommended for 
use particularly in "py oblems • of rings having one ' or more in- 
ternal brae i ng element s . i • • 1 * ‘ * -■* T 


: ", T NTH 0 DU C 1 1 ON • : : 

' • ■ <- * 

The methods and the formulas used in the analysis of 
monocoquo aircraft structures have been developed almost in- 
variably -for cylinders.- of circulate.', • ->or' possibly "Sl'llpt i c , 
cr 08s section and - - of uniform me'cha'n-i'CaT pr’op'ert itfs "Yet , in 
actual aircraft' such structural: e.lem'en-ts 'are - . Seldom/ if- ever, 
found. -1 nf ortunately ,.- the .dir ect ■ m.elhocbs - "Cf analysis are 
little suited to cope with. ‘ptr.ob laelm® involving- ’C otaiplerx croBS- 
sectional shapes, irregular distribution of reinforcing ele- 
ments., concentrate^- load'e.-, and 1 .cut-outs, > It Us believe-d 
that the indirect methods racfentliy.' advanced' by Hardy Gross, 
and particularly by B. 7. Southwell, (references 1 and 2) 
promise a solution of such problems. 
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The ■ fundamental ■ jus t if .lcat ion ? f pr.- this ■ indirect approach 
lies in the comparat iv.e-; stease ■ w$.-th ■ vh-i ch stresses in a com- 
plex structure -can he- 2 .a f lc 4 ila.'t..ed 1 .if •; they arp cau-s'pd-’ by - 5 ope ■ 
stipulated simple displacement , ■ while the determination of ■ . 
the stresses and' displacements' caused by. known, ext ernal loads 
is often very diff ieult or-' even impossible;.: The ■ calculat ion." 

of the stresses in a small 'unit of the structure corresponding 
to given displacements is known' as the "unit pr oblem": in« " 
Southwell's relaxation method. • .In this method the individual' 
displacements involved - in - the'unlt ' problems - ' are -bet>inbined“ by a 
trial-and-err or procedure until a displacement pattern of the 
entire structure • is "found' that cbfr espohds to'^an equilibrium 
of the internal stresses and the'- ejfct ernal load's over the en- 
tire structure!- ■' ■' ' • ’ *• 1 " r - 

•’ * * * ‘ ' " 

In parts I’And II (references 3 and 4) of "the present in- 
vestigation. a convenient unit problem was devigfed and solved 1 
for the pufpdse -of 'Calculating the stresses in sheet-and- 
stringer combinat i oiii * furthermore, 'a systematic procedure- of 
combining the individual' Operation^ was developed which re- 1 ' 
suits in a rapid approach to the actual state of distortions 
caused by the known- loads acting upon the structure. ThS:‘ 
stresses calculated by this procedure were found' to >he In 
reasonable agreement -.with thos e .meabur ed'. in experiments/ , ■"* 

In the present report the indirect method is applied to 
the calculation of the' b adding ' moments in rings (frames) upon 
which known external loads ar 6 Tact iag>. It 1 b planned.-to -dis- 
cuss in the futureithe s tr es s ’-.problem * of monoc.bque fuselages, 
which are combinations of -r ihgs and ' s tiff ened - curved panels, - 

' J * * v / „ : I .. .• f t ~ ‘ *. _ J i 

In the course of the pr es ent . invest igat ions it- was /founds 
that the calculations necessary for the solution of the unit 
problems are rather -[labor ipus , ~ For this- reason-, the solution” 
was worked out- numer ically .-and-.-the reaklts are present ed in * 
graphs and tables ip-part -I^ trefeSencS ! 5~); -Moreover,' it ‘ 
turned out that in many ring problems it is difficult to pre- 
dict -the -mpst'likely displacement patterns , and 'the order of. r -' 
magnitude .-of-' the -displacements , This' fact was "not aiitiei.-^'. ' - " 
pated- at ,t-he-- beginningt-cf the -pres ent invest Igationsf since ifl* 
the earlier work it was found that the displacements of rein- 
forced panels could be predicted with a reasonable degree of 
accuracy. It appears that with rings' the difficulty lies in 
visualizing the amounts of rotation of the end sections of 
the individual arcs, and to some extent the tangential dis- 
placement-^, ‘while -it: -is' -pos s ible t o' ant icipate 'comparat ively' 7 ' 
well the radial dis-p’taee-meht- pattern , 1 7 0 n the' other hand, the 
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interaction between - r otat i one, -tangential 'displacements, .and 
radial displacements is very : strong,-' "and 'the "final bending 
moments, shear ' force s , and- normal forces in the sections are 
often the small differ ence s ; of large .value s . caused by the in- • 
dividual distortions. Unfortunately, 'the Southwell procedure ■■ 
becomes very -slowly c onver gent ;• '-that is; it yields the correct 
final ro suit's ' oflly aft er ' a -very great number of- individual- 
operat i ons i f the ’steps in the relaxation. - or . t.he . individual 
di splacei'meht s - are undertaken -’at 'random, -without -a pr. eqon-r . 
ceived p-’i ctur^’ Of ■ the final pat.tern of deformations. 

In- ofddr ■ to • overcome this difficulty, a pr §$edure has 
been wo'f ked ' out whi ch 'mi ght be'termed the "proc§f.ure of the 
growing unit." In it the structure is broken up intp-units, 
the unit problems are solved, and the operations table is set 
up in the same. manner as was suggested by Southwell, On the 
other hand , ’ the • solut i on of the ’stress problem, by the trial- 
and-error -procedure of the method of ■ systematic relaxations 
is r epla 5-ed • by'a procedure involving . the. combinati on • of the 
individual iiaits into utiits of increasing g.ijze, -The external 
loads are balanced’ theft by solving sets of. two or three simul-- 
taneous li near * eijuat i oh s ih^olving the r influence coefficients ■ 
of the lar ge' tini t s . The -new procedure ib . rea sonably simple 
and fast, a "6* ihay ■ be Bees' from- tlie -numerical, examples 'contained 
in the body of the paper, 

i . • • • 

Finally* 1 the examples' present set;. ih thi s report '.were'.aibo, 
calculated’' 'by considering the so-called operations table/®JT 
the Sout h Veil"' me tlib'd' as’ a set of simultaneous linear equations, 
and solving ' i t ; ' by “■t'he use b‘f the’ matrix, calcults. In the<ap-. . 
pendix a simple explanation is given of the matrix calcula- 
tions necessary -for- the solution ’ . ■<' .. - - 


For the' under st andi’ng -of* the’; pre¥§'h : t report familiarity 
with the Southwell and the H'A'r'dy Cho-sfe' method's', or with parts - 
I and II ; b-f thi s' i hv e s t i gafi'd'h"' f-S ’• ft’bt reqUired'i- ■ •* 

r ' • i '■ l-‘ * » ,•*; * :-4X ;. - f .: ‘ .«? 

Thl’s 'inv'e'st i gat i bn , * 'c'bndUc t ed at-’the Polytechnic" Inst itute 
of Brookiynv + % f hS' ! sponsor-ed by' and conducted witty- the financial 
as si s-VancV- of t-hb' Nal-i-bha’l ’ Advi s’bry - 1 Commi’t t e e for- Aeronautics, ' 
‘•r ! ■> ? .? k.: rt'-, ■ > -s , fc . q .■ ’< .. > • 4 ... . - * ; 


*• i 


i ■ 

i ' 4 * 
\ ± . 


SYMBOL'S t( + i- 

’ 1 *. ■* .f 


cto s'S-fe'CC-t-.i onal area '.o.f : . b a r. ; o.i'.-'la point* o.jx? a., ring;, or 
. th’ar included area .of ft. r ing { ; or.- fit: matrix; - 
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A* 

effective shear ar efit (-based, on, tension) 

_ 


a 

lever 

arm of shear flow 



a, b, 

c, d, 

e, f, g matrix coefficients. 



B. C, 

D , E 

points on a ring ’ 

j — ' 


E 

Young* 

s modulus of elasticity 



& 

shear 

modulus •• 




I moment of inertia of cross -section; or" i’den't ity matrix 

j matrix coefficient * ** *: •' - ; -- 

L developed length. of ring segment 

M bending moment; or a matrix 

, ‘ i 

N end moment reacting pn bar or on constraint 

q •./ shear flow acting along bp.r 

R ,, ead. r adi al reaction acting on bar or oh constraint 
r , radius of curvature of a ring segment 

T end' tangential reaction a.cting on bar or on constraint; 
dr applied torque * : * 

U strain energy 

u " displacement of a Jpo^int. . in' tangential direction, 

v ■ ‘ displacement’ 'of a", point in radial direction 

V ■ shear force acting on a cross section 

w rotation of a section of a bar; or an unknown quantity 

x,y rectangular coordinate; or an unknown quantity 

... - • - ~ 1 * *■ •% ,i. * * ■ 

z an unknown quantity’ 

£ angle subtended hy ring segment j i o -ma r ■--- 

~ _ r ? - > i .. 
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y s ect ion— length parameter' (AL s /l) ‘ 3 

A determinant 

k = 1 - 1 - ( 8 8 /v ) [( I/O , + £ 'ii " r ' • '■ ’ ; , ; ' h r> ' h ' " 

\ = l - (8 2 /y) [(i/O i- i] ■' ' ' • • 1 

u = (G/S) ‘ ‘ . ; ' :i - ’ ' ' 

_ ... .. ty 

l ( A !|t / a ) ; or angular coordinate 

■ » E £ • * 1 * !?*•• * f ' «* * » 

cp ' *- angular ' cootdinat e 

. . , - ,■ f -M ^ 

(u rotational coordinate 

The symbols used to denote influence c'flBfficietts 'are de- 
fined in the following manner:. _ ... 

The term ('iTb ) stands for the force or moment a caused 
by a unit movement' in the direction' of b (whidh direction 
is^that- of the force H or T, or of the^moraent ^_N). Thus 
( nn ) is the moment due to a unit rotation! while' ( tr ) is the 
tangential force arising from a, unit radial^displacement . 

Fur thfer { o :i dis t ingu’ish the reactions at the fixed end from 
those at the mova.bl^ en.i,. the, subscr ipts . 3? and M are om- 

ployed. Oons equent ly '( nt )p is the moment arising at the fixed 

en<jL .. of „the purved 4 bargas., a. result ofa^unit tangential dis— ■’ 
placement of the movable end, while ( 1 1 J m - stands for the tan- 
gential force at the movable end due to a unit tangential dis- 
placement of that end. V 

The moment T _> ra&ia} ; for ce , and—tangent ial . f or ce caujaed by. 
a constant shear flow are denoted by the symbols nq , rq, and 
tq, respectively^ -The' values are ;valid for -the end of the 
bar toward which the shear flows. At the opposite end the 
reactions considered her e , are, act ing/ from- the support , -upon the 
curved bar. 

V , 1 ;• . - i ; .. • !<■ i * /« ? t "t »/ 

DEVELOPMENT.;} 0? .THE RE LAX AT ION If $ T^I,0|> W| TH V 5§E A jJDv _QF ^ AN '.'EXAMPLE 

The Structur ^ ^ a ASJ ^ C yoc;<ai‘ an e 

_ i t; e ,r y j > h i z v d ho A r - ■ ‘ d us 5 £■; n & ^ 

The structure used as an example in the development of 

the procedure is shown in fi-gure 1, It is a curved beam of 
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square cross section which may he thought of as an assembly 
of two cantilever quarter circular beams AB and AC joined 
rigidly at point A'.' • f As ea'ch of- :the, two cury.e.d, .b. earns .could 
support by its elf the 100 pound.' load applied , ; a, t- ’A’,, ' the dis—] ;' 
tribution of the load to the two beams cannot be determined' 
by the laws of statics ; alone *• ._.The problem is^ three times re- 
dundant, since the unknown stresses in "the cross section at 
A add up, in general, to ajfi unknown normal f or c e , an unknown 
shear force, and an unknown moment ih the 'plane of the circles. 
There cannot be any force and momedt 'resultants perpendicular 
to the plane, since the external load acts in the plane of the 

curved beam. ' " 

/ r • 

The problem of the load distribution can be solved only 
by taking into account the deformations of the structure. 

This can 'be done conveniently if- the structure is considered 
cut at A apd the deformations of each quarter— circle 'are 
calculated i-rEdepe.nde-ut.ly . . When- thes e .calqulat\ions are com- 
pleted, the continuity of the actual Btrb.c’tb.re can be re— es— " 
tablished without difficulty through a proper choice of the 
unknown moment- and forces : id Be'ction --.A-,' as .will be shown 
lat er . • i : 


i Tin i t P>r..o’b lem i> 

The s trui'turV is,- broken up into- t.wo .-.ufn i.t S’ ; ( namely, the. 
quarter circles AB and AC, . The unit pr.o-blem .consists, in 
finding the tangential (normal) force T, the radial (shear) 
force R , aup. the moment 'N ; (see fig. 2 ) und-er. the action 
of which the f r 3e : ''ehd 'p'o int A (of the unit undergoes. any pre- 
scribed displacement u in the x direction, v in the y 
direction, and any prescribed rotation w, which latter is 
considered positive if it is'- Counter-clockwise like ,u>. • The 
unit pr oblem„ pan he solved best by., calculat ing first u, v, 
arid w caused 'by- unknown forces T : and'-'R, : iandl ah \ unknown- 

moment N , and determining the unknowns afterwards so as to 
obtain -'the prescribed values.rofsithe displacements and^the ro- 
tation.' In conformity with, s.tandar d pract ice-, e.xjtfus .ional and 
shearing deformations will be; neglect ed , s ins e ■ ithpy ar piimu'ch . 
smaller than the bending deformations. 

' -s 

The bending moment M caused by T, R, and N in a sec- 
tion of the euirv e d •* b e am defined by the- angle ■ 9 .. (pee., f $.g ; , 2 ) 
is given by . ” 
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The i'elkfive' rot.a^ ion:. , dw of two normal s ect ions, o.f tbp.. . 
curved, beam an.^pf-Ju-lt es imal distance, d'g '/apart . • ; .. - 

■ ■■ •' aw .j= Mds/( El)' = Sfrd'cp / (SI)' ' ’ (Z) v:_„ 

-r.ii *} ■ . . ' ' ' ‘ ~ ' •• ••• . v* i I 

■ • i . . 

Since "point "B- 3 i’s* -'f-'ioced , r t.be -rotation *w^_ of point A can’ 
be calculated frdnr'the integral ' * , ' /" "' 

’ *' { . . ; „v ^ i ■ . 

;Tf/ 


. / _ . p n / ■* • ■ ' 

w_^ = / C Mr /(E 1)3 dxp>.- . 


Subst itut iob- n-f.-f •• M from equation (l) and integration yields 
W A ~ (r/3Sl)/C(»r/a). -- 1] Tr +- Hr + ( u/2)n|- ( 3 ) 

Wire n- -the infinitesimal element of the beam at cp under- 
goes a. ‘•r.ot-at ion dw , point A is displaced on infinitesimal 
distance du in the x' direction, where 


du 


r(l— coscp)dw 


(4) 7 


With B fixed, the toatl displacement u A of point A in 
the x direction becomes. 


n r 


J 

l p Tr/s . . p*?f 3 

iA = / du . = ( rfftl) 7 ' M( 1 — coscp) d cp 


Substitution, and integration gives 
u< 


= (r®/®-0. ^[(3;W/^) ^ 2] Tr + (1/2) Er + [(rr/2) -l] Kj" •( 6 ) 

•■i ■ « . 

When -the - inf init es,i.m,al element .of the beam at cp under- 
goes: an Infinitesimal rotation dw, point A is displaced an 
infinitesimal distanqe ’ ,,’dv in the y direction: 


'£Lv = r s in cp dw 

Cons ider g.t i ons s_im ijar t o those stated before give 
. rr / a p rtf 2 


( 6 ) 


-L 


dv = ( 


7®D 


M s in cp d cp 
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* 




f 


and 

Jt • 

v A = ( : r 2 /Sl) ^('l/2)''Sr + (tt/ 4) :Rr V* ' (?) 

f «... * 

* t ■ , . , f , r 

: ; . The task 'df determining t.he, di splacementrs of ‘the : end 
section A caused by the f.orc’e.6 ■ T. . and R and "the' moment 
H has thus teen completed. ' The problem must now be inverted, 
and the values -of T., • R ,. and IT. -- must be calculated that 
cause.' pr e scribed' distortions u A , v A , and w A . This calcula- 
tion can be carried out, .by-. solving the t hr eh: 'simultaneous 
equations (3)"-, (5), and' (7) for the unknowns T , R, and IT 

when the values of u A , v A ;,..and. w A . are given. In computa- 
tions that follow later in this report it will be found con- 
venient to have the solutions of three unit problems charac- 
terized by the following groups of prescribed distortions: 


(a)_ u A =..l 
( b ) u A = 0 


v A •» ! O: ; ■” r -'V A ”» 0 

.• • ' 1 '' 

v A = 1 w A -= 0 


( c) u A = 0 y A 0 . 

. . ■; > I, i 1 ) J- A '’ ’ 


v A = 1- 


The numerals 1 ■ ; depot e ; - unit 's inany convenient and consistent 
system of distances and angles. , ^ 

• - __ _ _ ^ - 

The solution of equations (3), .(5)., and (?) corresponding 
to values, (a) of t-h.a distortions is 

T = 43.87 (El/r 3 ) R = -39.37 (E_l/r 3 ) . H = 9.484 (El/r a ) 

When values (b) apply f the solution is 


T = -39.37 (El/r 3 ) R = 42.87 (El/r 3 ) N^=._13.03 (El/ r a ) 

■■ -r-n : •' - ' ' •' 

When values (c) apply, the solution is 

T = 9.484 (El/r 3 ) R = -13.03 (El/r s ) N = 5 ..44R (Ellrl ?■" 

In c 6'mput'i hg‘ the above value s a calculating mach,i.n% hqd; rt ?;0 . b &. :■ 
used,^ sgnce. ^malh _di f f.qr ancsis of large number s 'had to be de- 
termined.' The values can be verified by subs'ti-fjutijon in .equai-:: 
tiops (.3)., ,(_5;)..,.. and “respectively. 

... - -i * C ! >- C-5 ' ■_ " 

In subsequent calculations * ; use will . be,, mag. e ; of ; quant i-' ■ 
ties denoted.a-s "influence coefficient s." Influence so.effi^ ' 
cient tl; is defined-, as the' tangential force! $. caused, by k 

p 1 . ... v w j. ’I _> .j - “ - " 
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unit displacement in the tangential (t) direction; influence 
coefficient t? is the tangential force -caused' by a unit 
displacement ln ; ’the radial (r) direction; rn denotes the 
radial force B caused by a unit rotation (in the n direc- 
tion); and so forth. With this not at i on'- the solution pf the 
unit .. pro bl em s” may be stated concisely aa 

tt ; ,= 4S.87- (BI/r 3 ) tf = -39.37 (Sl/r 3 ) tn * 9.484 (ET/r 3 ) 

£r -tit ,4-3. 87 (Sl/r 5 ) fa = -13.03 (El/r 3 ) nn => 6.449" (El/r) (8) 

Influence, coefficients rt , nt , nr. are; not listed. The rea- 
son,- ;f. or the omission lieB in the ’bctfialities 

— ^ t tr ; v-‘ ; nl = €tL ‘ ' ' ’ nr = r-n, ... (9} . 

The validity of equations (9) follows directly from tko expres- 
sions given earlier as the solutions of equations (3), ( 5 ), 
and (7). 


The numerical values of. tjie , ; i nf luen ce coefficients will 
now be c omput ed ' f or the example ‘"‘at hand. The moment of iner- 
tia of the cross section of the curved beam is . 

,, .• • * ■ ? .0 r r • • .. . 

" I = (1/ 2 ) 4 / 1 2, = "<J .-00 62 1 no : h 4 

m * . • ^ , r $• t • f Vj * 1 . •« A 

The modulus o-f- 24 §t*T aluminum alloy’ ; is ; giveri =. JL0*.&-X 10 
£*1V '&4hce ' » • ' tv '• * ‘ 


i ■ .... .. El « -5^ i 800 ■-^euhd-i nch es 3 ! ' 

With beam AB • ’ 1 !f } 


(10) 


EI/r .546.0 -inch-pounds . E.I/r 3 
Consequently, 


= 546 poinds E“I /r = 54.6 pounds 
. { - r" . . j a f,.- pe.r .innh 


tt = 3340 pounds per inch 


~ : ' -fV I •: ^ 

tr = -3146 pounds per inch 


t;n = 5170 pounds . 


rs .f. --o ■ 
•r, v-s ~ - 


& 

3 S 'l 


: v. 


rr J =*'■’ 3.340 , p oupd-g' per = in<?h>(ll) 

ft -I . ' 't fi '! ' “ • / 


rn •*••-7110 pQu.^i^a: .;g. v.a i.XM-’' = ~§9 , 700 ‘inch-pound s 

!"•' " ' : •_ v ’ ( o o , .r r > o •: , ■ - , ■ - j 

With beam AC it is convenient to assume Ijhe positive 
senses of the .tangent i^l ,-f-q.r.ce ,, ^rad'ia'l force the moment 
as -eh’owh ^in figur.e., .3, i Tlie'.-p,'o's'lti'y'4' s'enses of"!, the .displacements 
and' the rotation corjre'oprond-' to th'bse : 'hf the" 'forces and the 
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i 


moment ; Since" the radius of this. beam is 20 inches, the in- 
fluence! coefficient e.. have ,the„ yq4jies^ . ^ ‘ ~ . , '* ' " 


tt = 292.5 pounds per inch 

^ O. <• sh.s r 

tn = 1292 pounds per radian 


tr = ' -268 .'5 r pound pen' inch" 
r"r = 292'. 5 p'd'u'n'd'fe 'pet:' 1 '! noh ; |Kl2) 


rn .=—1-777 ' p ounds per radian nn 14', 850 inch— pounds per 

• • ' ’ radian • > 


, • .( • : The Operations Table . _ .. , ^ 

* * * *■ •«.. U . ‘ . 7- 

' It, i b- possible. - now? t.o consider the effect of displacing 
point A upon the. pnm.pl at e .structure . S.ince in the complete 
structure there is no cut through the curved' “beam at point A, 
end point A of beam AO is always displaced in the same 
direction and through the pame die t-anc.e as i.s. end .p« 5 _int A of 
beam AB . A simple "unit operation" is defined as a' unit dis- 
placement in the direction -of one, of . „t-he.. ihre.e displacement 
coordinates (tangential, .radial,, r ot.at" i oiial ). while the other 
two displacements are held unchanged. The forces arid the mo- 
ment caused at point A during any unit operation can be 
easi-ly calculated- with the. aid _■ .of- /the. influence coefficients 
given in the preceding article. 


In order to displace point A 1 inch in the tangential 
direction and 'in .the positive- _sen.ee .indicated in ..figure 2, 
the following forces and moment have to be appl’led’to point 
A of -bar AB ; .... >. _. .... „ , . . ; ... .. e .. , 

’T = 2340 pounds' : B =. -2146 pounds ,N = 5170 inch-pounds^ 

In a similar manner the forces and moments can be ealeu- ' 
lated that -are -ne ce s Bary .t q displace end point A, ..of beam AC 
1 inch in the tangential direction. It is advantageous to 
adopt the . same , sign convent! 9, n ..for the forces and moments at .„ 
A independently of whether they are derived from beam AB or 
AC. If the sign convention of figure 2 is used, the pCeitive 
unit tangential -.displacement, of point A of beam AB corre- 
sponds to a negative unit 'tangential' displacement of point A 
of beam AC in the system of coordinates shown in figure 3. 
Mor e over 


® AC 



’Eac = Has ’ ' ’ ' ^ac = -Nab 


. ’! 1 - V ’ - " ■ ; •- • r v * * j f • * - • * . • l * . . * * * 

wher e , the ; . sub script ; AC,., refers to beas. -A^ and. the coordi- . 
nate system of figure-.?^ and subscript .■ ,...AB to b$g.m _AB and 

. . *> -* * f* i * f— • . ^ * f 


- i . ’* r 
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the coordinate system of figure 2. Hence th6 following forces 
and moments are ineceBsafy in order to cause poi-nt-’ A of beam 
AC to displace 1 inch t-o the rights 

T := 292.5 pounds R' = 268.5 pounds N = 1292 inch-pounds 

Altogether, the foilowing forces need he applied to point 
A to displace it (and consequently the ends of both bars AB 
and AC) 1 inch to the right: 

■T =: 2632.5 pounds R = -1877.5 pounds N = 6462 inch-pounds 


The effect of a unit displacement in th;e positive radial 
direction can be obtained in a similar manner. Beam AB re- 
quires 

T = -2146 pounds .R = .2340 pounds. H = -7110 inch-poundB 

With the System of coordinates' of figure 2 the forceB 'and the 
moments necessary to move point 'A .'of beam AC 1 inch upward 
are 

, '■•I .-•*••• ' ' " ' 

T = 268.' 6 pounds R = B92 :-5 ' pott'nAa.., N = 1777 inch-pounds 

Altogether, 


T = -1877,5 pounds Rj. = SGS&l 5^ ^oU'ndB.. if* .•*•6333 inch-pounds 

. . 3 fl e > ' * ' “ ’ j, r, « r "■ ' ■' J \ n ' 

Finally, in orde.r t q-.rutalt end A of bar. AB - in the 
positive N direction through an angle equal to a radian, the 
following forces and moment must be applied : to > point" A of 
bar AB: .. r: ’’ 

' ■ '' 4 .. ... . 

T = 5170 pounds R = -7110 pounds If, = 29^700'' i"nch T pqunds 

With the 'sign convention of figure 2 bar AC ; ; r eqttires fq.r the 
same rotation ' ' ' ! ' 


t 


. T =.'1292 pound.,8' <5. = 1777 pound s , N = ■■ 14 , 8.50' i n ch-p ound.p 

-• • . ... . "• 1 t. ‘ 

Altogether n • ‘ . ' c ..... 

^ — • * . . , ' > ‘ * 

. I; -■ ! 6'462 pounds R = -5333 pounds H = 44- r 550 inch-pounds 


It might be mentioned that' all' the values just calculate^ 
are f ict_.it i'ou s in the sense that they would., be. tke-Actual 'val- 
ues of the forces a$d moment s ■ only ' i f- tJ^e limit ?f proportioh- 
ality „o.f .jt he -mat er ia'T wer.e -higher than 'the 'stresses caused by 
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une i u ore e ana moment 8 ■, * ana ■ 1 1 TjUB-ittfgB u. j. s u ujt u j. uu b v. 

affsefcst'hfc , -'However; "both the forces 

di splafie&eit S' ahe^dlvi&ed 'by, ■ say ,' 100, no' 'ob jWtet i on s 

■ ”■ v . ..rf S h 4 : 4 ■ ( ' •* V ' . P P. 1 I ^ 


thelf opcSs" and * moment s £ ;: and "■ if • t he - ; lar ge'distO^tions did not 

and "the 
can he 

raised against ' the ‘tfdl idit'y 6f ‘‘the dat&~ ;}uet‘ dfilciilat ed . 

■ ■■ ) ? ;;fb i st ...;,. .? . .r i I ,- 

-The ^result’s ■ of "t'he'‘dalcniations afe BummafizS'd in tlie *;* 
tattler '.-t&et- follower '••'*' - ; v -* " " t ’ • * V ,I ’ ' Ifc - ” ‘ r •' w 

1 • •- ;• „ 1 ?-v* t . no r i J - * . 

it * :’..f s > „ . ‘ ■ 

-file- ; i : 


jj *! o o r* 1 


^ - 1 


* jt ;* 


: r j L .,-v. 

... '..r arable 

“ 1 . *' (Jper^Ct 1 ohs 

•Table z ‘ : y : 

■*. * »• ■ «• *> 

' .. ^ : ■ i( 

ic'-di: r 0 ■; 0 1 6' 

- - ■x 

\i V 1 * i£* i 

■v »*, »c *' r. ^ 

, ■ ' - • ’ b u. lqmsl: 

' R ' ■ 

nos..- . 

■ Operation i : ‘- 


- ’(’lb) v 

_ . r r •. 

; . t . ... .. 

» . * ,, • A *V 

I.* i, # . ^ V* ^ 

. 1-.. 

c *i 1 ■ 

u = -1-: inch '• • 

*hu2g3-g ;s"»‘ ” 

* +18 77 ;^ ' 

[2]' - 

v w - inch ■ • .- 

• y +1877.-6 

t* 2’6&2 

[3] ■ ■ 

w =■ rl' radian’- 

-'646 2" ‘ " ' 

+‘536‘3* 


- ■ ■ m. . . ■ . 

s * i - .J. * •'* 


*5 ’r 9j •* 

f i- ” , ? < - a 

i- i ! •* _ : . 

.* 1 r f ? 


tii 


•^(iv.yih) 


S-gWs'2. 


*- ,s * : + ; 5^33 f 

• :* . 3 . ^ r? * »r .,vf; » 

'.f'4>66tk 
, • - .•! r o(". . 


tiV ri 


The table is known as the "oper 
forces and moments previously calculated appear multiplied by 
(-1) ..Ifl^pt.-h^jp . words , • tha‘- forces and- mo*hnt"s lisbd&^af-e not 
those! jj o,-ith#r -applied to the" “Structure ‘ at' : p-&InV A ■“ in r ^'6r cl&r 't : o'' 
c au s e • . ,t h e - pip ,e;g e-r i b« d - & i s t o r t i 0 n s ; bu t £ & e i f rea'e'tloS'b/-' 

r- — r • er:-* . • r/rirru- • .„■ 


r> ~ n i 


t. e 


■j*s \ 


.•? i"T» 


i* on. 


. Eelazati on-'Tahlh ■ ■ J - 

• < • ; * l:.-. - ? ft' 


y If H h £ 3 I r O I 3 O I: if 

V ; / 2 o /j y r :■* 

c> J ^ > 5 r [ c - JL\ -■*!.. * V ; * I'C.-r. “ ? *^1 vt ^V u: t ^ 

?h§ operations.'-, table. . will-b'e ^tfSed'nstf-' l'n : ofdef to 'estab- 
lish. ftlje -di-St art ed. Ahape :th'e ' dt nn-c t'dr A -‘df ^ f igdfe' 1 i' a as sume s 
under ■‘the -Lactfclon i< 

It is 'coiiv anient: v 
hind -the v-sb.w&ct.tfr.e 

and that-vbijr a •:c-.lSffip-i , hg "d e vi'ce’ th'S hed > m-"^s 'ri'gf&'fy att'aclaad' f’o 

- - ■ 1 ** ' *• Sdo aSu 


the wall aft, iprO'inKt 



Moreover-', - t'-he'- fz e xt'e^naT ‘i<5ad ^majP^ft 


thought of as being suspended from the clamping device. 


The structure in its original, nondistorted form and the 
external load are thus in equilibrium, as long as the clamp 
and the wall are there. It is the purpose of the ."relaxation 
procedure 11 to transfer the load from the "constraints" - that 
is, from the wall and the clamping device — to the curved 
beam in a number of successive operations during the course 
of which the beam gradually assumes its final distorted shape.” 
A record of the individual operations is kept in the so-called 
"relaxation table" (table S) . 
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. , ^ . ? - *. i *'- ■ i 

In the .first row - of! tils r elaxat i on ' t able the external 
loads are listed. In the present example the only external 
load is the 100-pound "lc>a:& ivhi ch acts in ,th.e negative? rad'fal 
direction if, as before, the sign convention of figure 2 is 
adopted. Inspection of the operations table reveals that op- 
eration 2 is best suited for balancing the vertical load. 
Consequently, the clamping device is loosened in a manner to 
permit a jvertical translation of section A without allowing 
it to translate horizontally or to rotate. Since, according 
to the operations .table, 1-incb displacement upward would 
cause a vertical downward force of 2632.5 p.ound.s to act from, 
the structure upon .the clamping device, a downward ■■( negat iy.ej 
displacement y = , lip, Q/ (-2632.5 ) » -0,038 inch- .must 'be under- 
taken. The secbn'd -,rpw of the relaxation table is now filled 
in-with the values' of the second row of the operations table 
multiplied by (-0.038). It may be seen that f operation Just 
undertaken; namely., ,t.he vertical downward displacement;' of ',’. 
Q...-038 1 Inch caused t.he curved' ‘beam’ t o .exert-. upoa the : b'on- 
stjradnt' f a' Vert ical. upward force' "erf 100 p'ounds;., : -This foircV 

balances ' the 1 applled c load v.- ‘ 1 

. •# •* *. *.! * , . \ * #.* ” ■ 

- “ , . » «. .i '» *\ ) .' « • * * 

Unf ortunat ely , however, the displacement also caused 
the curved beam to exert upon the constraint a horizontal 
force of 71 pounds to t ha tleft , • c aVd' a clockwise moment of 
202 inch-pounds. Thus, the vert ical. force s are' -bafLaifce'd/' b^.tj ( 
a new unbalanced hop isp-nt-al fordoe- and a new unba-lahrcfe’d foment , 
are intr oduc'edV The-, sy-jetem-,' i ef' 1 ebu'r se , ' i s ■ in equilibrium,^ V' * '•, 
since . the clamping device and thb' ;; gi d.^n^l,.- ~ that.-, i e V' the _ \ , 
constraint- take care of the : hdri ^orjltaj ^.^or.ce and of 'the mo- 
ment . The purpose .of the .relaxaH f on ;p[f qcedu-r e , however ,^‘i s 
to a.pproach gradually a -state of- 'di at dr’t'i ons in which? ''artifi- 
cial constraints are not-, needed .for' 'equilibrium. It fa im- 
perative, therefore, to .balance out the new unbalanced for.ee 
and moment. However, before any further steps are undertaken, 
the clamping device must again be tightened completely^ ■ *"■ r: * ** j‘. 
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. • • . 

•- -Sable-. Z-i. r Relaxation 

Sable • 




‘ ; _ 

, ;■ 

*' • ** ' ■ * N 



K*V . \ 

t it 5 - r 

.J ( in .-lb ) 


E-xt ern^l • force 

• ■ iff.- ,....0 : : • :■ 7 

-ioo. : ■ 

- '--.-O 


-0.038 x [2] 

" .‘,71 

; v •* ■;* 

100 

_’"l;_-gPS 


-71 . • •«. 

l iQ‘ i. 

e.i;. *202 

, .. 

-0 . 027 x fl] 
< ■ 

- ‘ 1 *>*? * 

* ■ . 71 > ; • 

I 51 

‘ r, -.; r . ;:175 . 

- 



— 

: * . ^ tj * O ’. 1 '' 

1-51 

= :u:t. T! -- g? 

.. =- j. . t i - • - — ff ■ 


L-V" 0 ' 1 * 

-0.0194 *X .tsT . 

* i ‘‘ «• .•» 

■ ' ; ■ -36 

.-36 • • . ‘ 

_51 ' _ ' 
0" 

• ' \ * ,; Vl03-' T: 


-Q:0^37 x [1] > 

' l ’ ’ V V‘_36._ - 

-26 

*' 88 

* / 

0 

-26 

-43 


-0.0099 rx- ['2j 

*■ ; . i •* ,* • • • 

<’ * •* . -19 ’} ■ - 

_26 _ 

-53 



.*• * t ■ ^lS’ r* 

o • 

-95 


-•0.00722 -X [1] " 

4 ' 19 • ■ 

'-14 

^47 


^ • 

. - • 0 

-14 

48 


-0.00532 X ' [2 ] 

*• • -io 

.14 _ 

- - - =sa_ 



-10 

0 

20 


-0.0038 X [1] 

10 

7 

25 


0 

-7 

• • . 45 


-0.00266 X [2 3 

- --5 - - - - i 

-5 ■ 

2 _ 
6 

-14 

: 3i; 

r _ 

0.000696 X [3 ] 

t-4 

’ll'- 

' -31' ‘ 



■ -9 

' 4_ . • 

o' 


-0.00342 X [1] 

9 

_-6 

22 


- • _ i * * * 

0. 

- f ‘ ■: . - : • - 1 - 

-2 

... , , . -; 22 


, » *■ l 

r . , .. - ^ . 


■ ■■ J .» , 


O.Q00493 X ;.t33.. 

: JL .■* ■■■—! ,T— 0— — > J— . “— 

_ 3 

- • . - -22 

- t _ 

, . * . /-f 

1 -- J, - • -r 3 < ; v C 

iB-i’ : 

1 ■■■■•■■ *'.fr 
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I lie unbalances are the algebraic sums of corresponding 
values in rows 1 and 2. They are listed in row 3, under the 
dashed line. The question arises now whether T or N . „ • v 

should'be balanced first. Because it is expected' that:..th& 
final distorted shape of the structure will involve consider^ 
able horizontal and vertical translations of sectioh A,' r but 
only a slight rotation, it is advantageous to concentrate 
first on 'the cancellation of. the unbalanced forces rather 
than of the moment. For this reason, operation 1 is now un- 
dertaken. The clamping device is again loosened, but now in 
such a manner as to permit a horizontal displacement, but at 
the same time to disallow any vertical translation and rota- 
tion. The ensuing horizontal motion will come to a stop when 
u m 7l/(-2632.5) e -0.027 inch. Multiplication of the valufe'B ' 
listed in the first row of the operations table by -0.027 and^_ 
subtraction yields the remaining unbalanced, or n re sldual , 11 
quantit L e s "recorded in the fifth row of the relaxation table. 

According to the fifth row there 'are no horizontal un- 
balanced or ce s acting now up'on the constraint, but an unbal- 
anced vertical force of -51 pound appears again, which 1 b ac- 
companied by an unbalanced moment of -27 inch-pound. Conse- 
quently^’ the .tw.o displacements undertaken succeeded in reduc- 
ing the • unbalanced vertical force to 51 percent of its initial 
value, but only at the expense of introducing an unbalanced 
moment at the same time. The relaxation must be continued, f 

therefore, until all the residual quantities become small 
enough to be neglected in engineering calculations. 

. The! clamping device , is tightened and loosened again suc- 
cessively, but always in a manner to permit only one single 
type of motion at a time. The effect of, these motions upon 
the forces and moment is calculated in the relaxation table.. - - 
After nine operations the residual forces are T » -5 pound, 

E = 0, N a 31 inch-pounds. It is thought that now it is time '* 
to permit section A to rotate. Thus, operation 3 is under- . . 

taken and recorded. ThiB is followed, by one more operation [1], 
and again operation [3], after which the residuals are small . •' - 

enough to., be neglected., • - • ' " , 

Because of the great number" of* arithmetic operations in- 
volved in the preceding calculat i on s ... it ' i s- well - t o check the 
results. A'repetitioh of : all the computations would be very 
time conjsuming. Fortunately, this i s-unnece ssary 'since a 
much simpler check' is available. The sum of all the vertical 
displacements can be obtained by simple algebraic additions 


■r a - 1 


t Ot 


* -0.07528 inch 
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1? 


* 


* 


Similarly, 





[ d ! - . rt i 


u t0 ^. * -0.0551,4. inch. 


■ r* - ‘ O'. 001X89" raHT an 
tot 




, t . 


The forces and the^moment caused ..by—ea-e-h'-one—of these total 
distortions can he obtained by multiplying the corresponding 
row in the operations table. “The ,vdlues computed are entered 
in the “check table. “ The algebraic .of .each .column of 
the check table represents the final residual quantities. 


For the present example these' calculations are carried 
out in table 3. In the absence of errors the residual quanti- 
ties of the last row of the check table should be identical 
with the- residuals listed in the last row of the relaxation 
table". This is. obviously not the case with the present exam- 
ple . . : 

j ~ ' ‘ ‘ 

0 Table..3..- Cheek- Table 


■4 •. ' ** 

i ‘i r r <, f 2 , is 

.external force _ . 

u =- -0 . 055'1 : 4- "inch' ' .. 
v = -0.075)88/ inch . ; ' . 
v «= 0.001l8’9 radian 
Residuals 


* ->■ r .. ; -. • 

(lb) ' • ..(lb),;.. ; . (i-.n .-lb) 

O ' .,- 10.0 = . ‘0 

145. . .--104'' ' 356 

. '-14 It' '19 8 -401 

' ' " -8 ’ 6 • ' ~ “ -53 

-4 0 -98 


Fortunately^ /there is no nee<i for hunting for' errors. 

Instead of doing so, a new set of relaxations can'be carried 
out starting out from the residuals of the check table as 
the given esternal. Ipfdingj. i.Iable 4, the second relaxation 
table, presents these calculations. It may be seen that after 
only five operations the residuals are reduced to 1 percent of 
the initial external lo.ed . 
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Table 4.- $;ec.e*nd .Ralaxat i on Table 



T 

u, Mb) 

7 ;• *f * 

H 

‘ (ib) 

N 

(in .-lb) 

External force 

-4 

o' ' 

-98 

-0.0022 x [3] 

V . Liv _ 

-*r. s .10 . 

r -12 
-12 ** 

98 

„ . ’ ' ; b • 

1 . • _ . “ -i-* : . • a 

-0.00456 X [2] ' '' ,,n 

1 r • 

•- 

1 

- - ■ - . ■ jT . - - 

... . •' -j 1 ' ]' d 

_ .1 2 J J’J 
0 

tsc;? o.i ' .» . . 

-■ L ; x ,t 

1 — 1 

to 
1 1 

X 

iO 

0 

0 
0 

0 

1 

^ ~ ■ 3 
4 ’ 

-3' 

V, r ' \5 V ~ 

• ■”> ' 'J t Qt * 

+ 0.00168 X [1] 

- ' _“4 
0 

3 ' 
6 

. "" ■' e#.' 

- 5i#JLjO '4 ^ J- 

~ '"-l.O-f.~n/ 

-.0.000225 X [3] 

1 

I 

-1 . 
-1 

10 

0 


. , f.'-'- 1 .. 

Those results are -now cheeked 'in the second. ( £beck table 
(table _!>•). The residuals of this ar.e sub stan.td^ + l^'* the same' 
as those of the second relaxation table. The ■ .s.ljlgT^b.'.d'evia-'' 
tions can be explained' by the inaccuracy o f. ,s jLtid^ rhlje. palcb-v 
lations, and by the omission of fractions of pounds' a'hd' i-nch->’ 
pounds in all the tables. " - •‘■'■'n • & . 

i , 

The fact that error's need' hot be traced but ,can Jpe „,elim- 
inated in a few additional Qperati one.-i s co^sidere^ on# "of 
the major advantages of the relaxation -prosedure ; r 4 bo * , 3 a j 


» . , ' . lf / i0 e .ft 

Table 5.- Sbcond? Check Tab! e ; ‘ 




* r' ’ . . 

i’X -■ j' :< c 


T 

E 

. .-/■■■ < ■* 

N ' 


(lb) 

(lb) 

(in. -lb) 

External force 

0 

-100 

0 

u a -0.05362 inch 

141 

-101 

346 

v s: -0,07984 inch 

-150 

210 

-426 

v = -0.001776 radian 

11 

-9 

79 

He sidual s 

2 

0 

-1 
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Direct ' Solution- of • the' Operations Table 

The operati ohs ■ table -represent s ' k sAt of- linear equations. 
It is easy to verify- that ■ table- 1 cah be * r'etor it t en in the form' 

. •: • . m . •: c- . ■■ i a:-. ' • • 


.1 . . ; . . • .•■ k : v ' ’ H.-i? 

<- 2632 ; 5 u- + 1677. 5-' v -'-64-62 w + T A"*) 

1877.5 u - 2632 . 5 v + 5333 w + R = 0 
-6462 u + 5333 v - 44,550 v : + N. =>' : 4) •• 


( 13 ) 




£ ' \ 1 » ■ V % * 'j ^ . * ■» * 

where u, v,'w are the unknown displacements add rotation, 

are the given, ex- 
*** e " ‘ the" vklues Of “ 
(13) be- 


respectively, of section A; and T, R, N 
ternal forces and moment, respectively, 
the external loads of the present example , equat i ons 
come : . . 






-2632.5 : ;tL'V 1377 .5 v - 6462 w'f. * 


A 


1877.5 u - 2632.5 v + 5333 w = l66' ; ), 


-6^62.-u .+ ,S333 ^t- - 44,650 w = .0. 


j 


0 v 


(14) 


These equations can be solved direct ly .’ The solution is 

u ='-0.0513 Inch ’ v = -^0.0785 inch Vr '= -0.00193 radian 

V‘‘ . ‘ ' ' * . ’ . ' (15) . 

Substitution in equations (14)' shows- that the ac&tifacjf is 
better than 1 percent. The -agreement ■ of t'hi s rigof ons -.■solu- 
tion with the results obtained in the relaxation procedure is 
al bo good. - 


It might be argued that there is no need for the relaxa- 
tion procedure if the simultaneous equations represented in 
the operations table can be solved directly without difficulty. 
This argument is justified as far as the present example is 
concerned. In many other, more complex problems, however, 
when the number of simultaneous equat i ons - ihv olved i s very 
large, the relaxation procedure may represent the quicker and 
easier solution. ■ :.0 T - 0 . \ ... . 5 • . - ■ 
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CalcuLat i'o'n ': 0 .f the, Banding 'Moment s. . ; 


The -he'n'dihg moment in secti on ;A .can .be easily 
lated VitK ’the- stfd' "of the .influence • c oef.fi ci ents of 
given in equati ons ( 11 ) , or those of tar AC given 
tions (12). If the distortions listed in equations 
used, the following bending moment. . is obtained 

tion A of beam AB ; 


qa 1 ou- 
bar - AB -■ t 
in e qua-- — 
(15) are 
for sec- 


N« * u fit + v nr + w nn 


= rs0,0$1.3 x. 5170 + 0.0785 X 71.10 - 0.00193 X 29,700 
, ' * * 1 ( 

’ ' v - * •» * . ... * - • '* : 

» 235 .4 lnch.-po.unds ’ . 

f ‘ ' J • •* , ti i‘ ' 1 

If the di solacement u and the rotation w are attributed” 
positive signs |n _op.d.er ; t.o comply . with the sign convention of 
figure 3 and equations (l2), the following expression is ob- 
tained for eeotion A of beam A3:. 


N a » 0.0513 x 1392 +''o . 0785 ' X 1 777 +’ , 0t'00193 • x 14,850 
= 234 ,.4. i.nch T pound s „ , .. . .... 


The .rtFo: values differ ^yjjless than„l percent. .The dif- 
ference is considered permissible since the calculations were 
carried out with a slide rule. The normal force T A and the 

shear -fcrrae ‘can be-.. c'al<?ialat ed in a similar • manner , . When 

the 'calculati ons &fe "based • Oflf. bar AB, . • •/.•> ; ■: .. j 

e i. ■: ■■ . . V »■****• r*~ - : -tv 

T . = u 1 1 + v t r + w tn . f ^ 


7 • ~ -i >'■ » "* ' S' 1 -’ .> •.{ ", 0 <s *rf 

-0. 0513-, X.- 2340--+. On p.795. X 2146 - 0.0Q193-.X .5170. : f 


38 , 4 pounds. ■, 


► „■.« nf-i.fi • >. * * ,*,• 

- ■ b-; x'i 1. ; csrr. ? ’ • 

, •; ... ••• O v. 7 ft 1,’w t v: . '• 

®A » v u«-rt.: 4- V\ rf v <*■. $r rn; 7 e r. -.ij 3 ' Tt -. • . - v 

' sr - 1 f.ii T • »•» Me:* - £ • 1 

= +0.0613 X 2146 - 0.0785 X 2340 + 0.00193 x-7110 ‘ 




-60.9 pounds 
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When the calculations are carried out for bar AC, 

T a = 0.0513 X 292.5 + 0.0785 X 268.5 + 0.00193 X 1292 = 38.6 
E a « -0.0513' a 268. 6.^ 0:0t85 X 2^,5;V 0.00193' X 1777 » -40.1 

< ", .r 

The deviation of the 4 tensile ■ force obtained for AB from that, 
obtained for AC is much- less -than .1 pe-rcent of tha load. 

The two shear forces add up "to ’+101 -.'pounds* instead of>the cor- 
rect' -100 pounds. ■■ ' ; * iT . 

••Tlie fnc^beam ' AB ■ can be caiculatefd” fronP~ 

equation (l). In order- -to ■ ge't c oh si’ st ent -result s , ‘the 'follow- 
ing values are used in the equations: 

. . . -• r - i->_ . 

N ss 234 .9 ' Inch-poun'ds T •» 38.’5 pounds ®-aB * -60.4 pounds 

* * - r 

. i ", R ao b -.39:6 podiids 

• . m t 

Bauat I on (l) ' becomes •: • 

f 

JT * 38.5 x 10(1 - cos cp) - 60.4 X 10 sin cp + 234.9 

-:v. j ' • ; < 

- 619.9 - 385 cos cp-- 604* fein. <p ‘ 

Similarly, the bending moment in r tjeam"" ■'AC ' is: 

M w 38.5 30(l - cos :<p^ j- '59 . 6 .3? ’20 sin cp +.234.9, 

» 1004.9 - '770 3^s;cp-’- ; ;’793; sin^ .. ' ‘ 

The bendipg ; mement ..diagram i s .- shown in .figure - 4 . -s s 


si." * 


f . - 


A. — ~ ' . * - 


CALCULATION .05 THE,!. B.END-I NO MOMENTS' IN RiNG-t '6l AR^ITKAR.? - S SAPS 

-■ : i. 1, f-- - , -> ii.: t -i- J f -- * : *. c . “ -I 

... .. ; ,Eevie.w ,:of Methods, of ,,Calculatf on 


j : .' 


When ..the . shape ... of the median line ,,q.f . v i;he .,f raiie. ._devfat e. s ^ . 
considerably ,, from a i-circle, or- ellipse,,- ,or. t> whe.n there ..are, e,ddi-^ 
tional , internal bracing .members incorporate.^ .'in the ’ fra.m.a f . ’ tK© 
classical analytic methods of caiculati on ; .can f .seldom .be used. • 
The numerical or graphic procedures ; sugge. feted. -by Lu,ndqui.st and 
Burke (reference 6) and Hoff (reference 7), and the column 




. * * • 


.t . . * . 


3 ■ 
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analogy developed by Hardy Cross and .recently described *f'or 
the benefit of aeronautical engineers by du Plantier (refer- 
ence 8) are .useful in the former case - that is, when -the 
moments have to be determined in a closed ring without inter- 
nal bracings. In the publication by Lundquist and Burke, 
equations are- also presented for the case "of a rihghaving 
one internal bracing element. It is believed that the pro- 
cedures described, in the present report will be f-ound reason- 
ably simple even. wb,en several internal bracing members qre in- 
corporated in a frame of arbitrary shape. „ . 

It will be shown in the next article how the relaxation 
method .can be employed to calculate the bending .m-cunant s in an 
"egg shape" frame loaded symmetrically with two equa'1; and, -op ; - 
posite forces. The egg shape in question is the combin.a.td, o.n 
of two arcs of circles and two straight lines. It is believed 
that the median line of every fuselage frame can be. ,jrjep.re-. ... 
6ented with sufficient accuracy by arcs of circles and straight 
lines. 

As the determination of the influence coefficients of 
arcs of circles is a cumbersome task because of the small ,dd £rr In- 
ferences of large numbers Involved, tables and graphs of in- 
fluence c oef.fi ci eptp. .haye been prepared and presented for-p.se 
in part .TV of* the' present investigation (reference 5). This 
part IV, entitled Influence Coefficient s. cxf .Curbed- Bar s -/or 
Distortions in Their Own Plane, is intended for use in future 
frame calculations in the same, manner as tabjLes of trlgopo- . 
metric or hyperbolic’ fun'otions' are us^d'. ; ’ ' 1 ' ~ ' 

It was found'tha-t a' - r e suit" of suf*H.c i ent ’ accuracy can be 
obtained by the relaxation method in a reasonably small number 
of steps only if the c omput er ~ ekin' make k” good gues's regarding 
the distorted shape of the structure. Unfortunately, it is 
hard to vi suali ze the atfrount s of ..rotation*- -fthe 6nds" Of’ the* in- 
dividual elements of the frame undergo, Por this reason, a 
second type of numerical approach is presented. This uses 
the epOr a t?.i oh s table of the. relaxation' method as - it s smarting 
point and reaches the answer to the problem through $ number 
of solutions of two or thrdd 'simuli'&nedue -linear' equations. 

The procedure is denoted as the "growing unit" procedure. It 
is applied to the present problem and gives substantially the 
same result s as.^ttfte ; re'lakattdh method 

•- ” .• . n , ; -re h ' 

Por a-- -! check of the results thd^rdblem was 'also solved 
by the grs.phi c 'p'r-ocedure df reference ?r' In addition, the 
system of -si mudt arte ous equations rejrre's-'ehted by the operat i ons 
table wa.s.'S.qlv.e.d ( wi fch ,tihe -aid -of the 'fjfiat-rix cal cuius .' ' "Thi s 
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latter approach Is discussed in the appendix, hut the results 
of the calculations are presented in this chapter for the pur- 
pose of a comparison with the results reached by the various’’ 
other methods. The agreement is found to he good. 

finally, the results of one more comparison are given. 

It may. he seen from reference 5 that the influence coeffi- 
cients phi opiated ' on the' oasis of inextensi onal deformations 
differ considerably from those obtained when- the extensions 
and- the shearing deformations of the curved bars are taken l;i 
into' account . Ih’the last solution o.f' the pro sent problem 
the. operations table was recalculated using the "inextensional" 
influence coefficients. It- turned out that the bending moment 
distribution obtained? -through a matrix solution of this modi-, 
fied operations table was practically the same. &b that calcu- 
lated earlier on the basis of the more accurate operations 
table. 

In all the problems discussed in this report the cro.se- 
sectional (dimensions of the ring are- ■ consider ed small as com- 
pared to the radius of curvature -of the ring... Consequently ,. 
the distortions of the ring elements can be calculated from 
formulas based on the linearity of the stress distribution 
rather than the hyperbolic law valid for curved beams. Simi- 
larly, in problems involving shear flow the lever arm of the 
Bhear flow - that is, the distance of the sheet from the neu- 
tral axis of the ring section - is neglected. 


Basic Data of the Egg-Shape Ring 

•The dimensions of t he . ring.- are t given in figure. 6. The 
6hape of the median line is taken from- reference, 7 and can be 
obtained by drawing the two circles of 20-inch and 26 -inch 
radius, respectively, with their centers 25 Inches apart, and 
the common tangents to the two circles.. It- follows from the 
geometry 1 of the -figure that the angle-, subtended by,, arc' AB is 
78.46°, and that subtended by arc CD is 101. 54*. 

It Is assumed that" the frame i s . .manufactured' by bending 
an. aluminum alloy I-section to the required shape. The area 
of the'’ 1 I- sect ion is 0 ; . 610-. .aquar e -Inc.h, its maximum moment of 
inertia 0.952 inch 4 . The deformations of the r.-ing pnder the 
loads depend upon, the bending, extensional, and shearing rigid- 
ities of its section.-' I -It. was -.shewn in reference 5 that the 
shearing rigidity has a considerable effect upon the magnitude 
of the ’ inf Ihence- "ooeff •ilaie.iwb-a-. . An "-.effective shear area" A* 
was, therefore, .defi-ned. :-i£n --.such a; manner that .the' 'strain 
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energy in shear stored in an element; ^fpth.e-'hea'® of ‘ah infin- 
ite s imal • lengt h* ■ ' &1 1 "could "be . caic]ji<|tgd.-> fv ok . the' simple f or- 
mula - -Y" ~ : ' ' ri ; “ ~ ^ Y. ‘ : 


■■dir 


shear 


(V 2 /2i5A»)dI 


where V is the shear arc - e in the s'ecti/on and B p i ; s Young’s 
modulus. The. .rat.i o.-> f A^/A > where A is the actual .cross- 
sect i on'al’ area’ of the 'beam', was denoted by | and may be 
termed the "shear ri gidity fact of . " It dependj upon.- the shape 
of the cross sectionand fche-value of d-/35 , Vhere 8 is the 

shear mo <jluiu s . . jn.ref ersnee. S', f ormula s , were developed for 
the calculation . of >the shear -ri gidity ' fac’t or for some repre- 
sentative cr o s 8-saat i ©nal..LShape s . ' - 

In the case of t.he .I-sd'ct i on shown in figure 5, the value 
of | can. be calculated from equation (.10) of reference 6. 

The computation gives A* = |A = 0.081 square inch. It may be 
noted that with the -I-section in question the effective shear 
area is almost exactly the area of the web multiplied by the 
rati o G/E . ' 1 ; ... 

In the calculations that follow, E is assumed to be w 

10.5 x 10 6 psi . With this rvalue EX becomes 10 7 pound -inches 3 . 

The ratio G/B was taken -as 0.385. . . - . - 

i 


Calculation of the Influence Coefficients 


Since all computational work was carried out on eight 
bank calculating machines, not more than eight figures were 
kept in the calculati ohs . - 

Arc A3 .- The angle subtended by the arc of circle was 
denoted 8 in reference 5. The parameter £ was discussed 
in the preceding arti<?lej.: In the case of arc AB- .the values, 
are ; • . 

* ... - • . 

P * 1.3694384 radians t * A* / A =* 0.133,1058 

The sect i on-l.eng.th parameter Y was defined in equation ( 1 S ) 
of reference 5. . la, the case of arc AB . - : 

j--.- ■--- . Y . s* AL 1 « 48 0.8 619 , ■ 

where L is the length of the arc. The parameters « and' \ 
were defined in equation (19) of reference 5. Substitution 
yields 
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1C = 1 +• (8 S /Y) + 0 2 /V I), a- 1.0332 
= 1 + (g S /Y) - (.p*/Y|) = 0,9746 


The influence coefficients are obtained as ratios of the nu- 
merator and denominator determinants. First, the denominator 
is computed from equation (21a) of reference 5;. .j, * 

. A - 0.00521084 r* 

•'I ’ ■ I = xu: ■ ' 

Equations (2lb) to (21g) of reference 5 yield the -^numerat or s ; 


vi-imi-r-'i * X « ,r 


A_„ *' 0,0198 6728 ;8tr 
nn 

• A ri i * -0.04079544 ~ - 

‘ =’°V 6273076 2 IS'"-- • 

- : Ay r =■'§. £3958066 (Ef/r) 

_ A^-j. ? 0. 143^024? (81 /r) 

A tt = 'C;198042S4MSI;/£? 


•? » ' d « 'if 


. ~ {•'• li C- ? 

. •} i. ■' • 


"‘Tliei influence coefficients are notf ' dbt ained by division; 


, ■> ». ,“l IJ 

■* *. ' ^ ryn c. 

Tn 

T-f • r* 

- * - 

S 

4»»/4 

• • . - t*f 

: , f K <t : . 

rn 

s 

Apn/ A 


tn 

- 

Atn/^ 


rr 

= 

A rr /A 


tr 

s 

^tr/^ 


tt 

s 

Att/A 

A unit 

displacement is 

rotation as 

10 -4 radian, in 

in the computations. 


Hence 




/A = -7.82896 (Bl/r®) 


2 . 

* 

3 - 


(16) 




.-4 


p 1 i edL by ■ 1 O' 




With 
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5*0’ x ' l-O 4 i'nch-pquhd-s '■• 

- •• qB-J / r 8 ( *' 'S'! 5 1 X - 1 0 4 • pound e -• ■' 

ai/r 3 = 0.125 x 10 4 pounds per inch 


■ .1 i 


•; « * • . c :,*ji * * :■ .■ * '3 

« * 

the final values of the ^influence cOSf f idi ent s 'become: 


^ 1 ‘ - 4 

nn = 190.634 inch-pounds per 10 radian 

rn = -19 . 5724 pounds per 10 radian 

tn = 13.1Q1$ pounds: pep. 10 4 .radian 

r r = 3 . 3 4 8 3 3 j.p ound-s p e r , 1 0 4 i n c h 

tr a -3 . 4352,0 .pounds per 10~ 4 inch 

tt a 4.75,07 pounds per 10“ 4 inch 


(17) 


Straight bar' BO .- For' the straight bar simple formulas 
are given in equations (25), reference 5. Symbol £ was as- 
sumed to be 0.13305785. As long as th/fe'' same value of £ is 
taken consistently in the calculation of the influence coeffi 
cients of any one portion of the frapp., it ,i s , of .11 1 tl.?- impor 
tance whether" the same value is uied'for all the por- 

tions. A change in £ from one portion to another is equiv- 
alent to a si ight . va^i_^t i on in the crosB-s^tional shape. 

The difference inPhe^two values of ■ £ given corresponds to 
a changet in the grp se-r sect i onal shape thpt^fs well within man 
ufacturing t o.lef a3ce s .' The following : dat a of the straight 
bar are needed.: , . . _ 1"., 

L = ?4,.4949 l.nches. V = 385^14588 


Equations (25) ©preference ,5 give: 


\4 


nn » 137.232 lnch-pound3 per 10“ 4 radian 

rn = -8 .10-2.€7 S .pounds per ; liO“- 4 ' : raidd.>ah' -■ r /■ i i 

rriefdo ot -1 i-' ? s.'. ?'<.:*• ■■■'£ ■■ ■' 

tn ®-B as inti :i BX'I uv t ?:\ ; , *«.. ! * <.3.' 

■ f d 1 % . 

rr = 0.66158 pound per 10“ 4 inch 
t r = 0 

tt = 26.2058 pounds per 10“ 4 inch 


( 18 ) 
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an n 151.135 inch-pounds per 10 radian 
rn = -11 . p^gi89 ,pounds per 10 -4 radian 
tn = 5.24667 pounds .per.lCT" 4 radian 
rr = 1.24047 pounds .per . 10 -4 inch 
tr «s -0.88929 pounds per 10” 4 Inch 
tt = 0.8774 pound per 10 -4 inch 


> ( 20 ) 


;; f <? v 

In the example of 
to he established only- 
fixed en.de B and 0 
moments were available 


»t r ' ^ r », •.% r j 

* llbr ium" had 


■ -The - ' I nf luence ^.C.oeff i ci ent s at the. Fiy^d. Epd^ 

the example of the preceding chapter , ‘equil! 

at point A of figure 1, Bince at the 
of the arcs any reaction forces and 
. For this reasoii, only point A waB 
balanced in the relaxation p.ropadure, and the effect of dis- 
placements of point A upon the fixed ends of the area was 
not investigated. The: situation is different in the case of 
the- frame shown in figure 5. There three points A, B, and 0 
are displaced in the r.'el axat.t ■op -procedure . Consequently, the 
forces and moments caused, for example, at A and C by dis- 
placements of point -3 : must be- de-t.ermlned . As in the pre- 
vious example, only a single motion is undertaken at a time 
in the relaxation pr oqedui'e For instance, point B is dis- 

placed radially and at the same time it is prevented from 
moving tangentially and from rotating. Simultaneously, points 
A and 0 are considered rigidly fixed in their positions. 


The influence coefficients calculated up to this point 
represent the forces and. moments necessary at the point that 
moves, while the other end of the curved or straight beam is 
considered rigidly fixed. Simple equilibrium considerations 
suffice for the calculation of the ’ react i on s at the fixed end, 
when the forces and the mdment. at the movable end are known. 
These reactions also can be reduced to the influence coeffi- 
cient form. In order to distinguish between the two groups 
of influence coefficients, the. former are denoted as the in- 
fluence coefficients for the movable end, the latter as'those 
for the fixed end. In the formulas the subscript M refers 
to the former, subscript F to the latter. If there 1 b no 
pdssibility for misunderstanding, as in the preceding example, 
the subscript may be omitted. 


In determining the sign of the forces and moments the 
beam convention shown in figure 6a is used. 



NACA TF No.' 998 


29 


The influence coefficients at the fixed end were- calcu- 
lated with the aid of’ equations (27), ( : 28), and (29) of', refer- 
ence 5. for arc 'AB their values are: ■- • • ■ - - - - ; 


any = 0.33433 (Sl/r) 
rnj = -3 . 56 8 8 7 (El/r 2 ) 
tn F = 8 .71889 (EI/r 2 )"- : - 
rrj = 21 (56 91 (El/r 3 )' 


. 


tr F = -31.7417 (Bl/r ■) 
tt F = 34.5276 (E i/r 3 ) 


Substitution gives 


( 21 ) 


nn F * 16.7165 inch-.pouhds p$r 10 — 4 r&dian 
rnj - -8.92216 ppUndsper 10‘ r ' 4;; radian ” 
tn F * 21.797 pounds per lO" 4 radian 
rrj = 2.69614 pounds per 10~ 4 inch 
. t r F \ '= ' ^3 ( 9 6771 ' p ourids per 10 —4 In6h " r ’ 

' tt F ’;' = 4 ;31595 pounds per 10 —4 'iiich ' ' •*. ''■••• 

Bor the ' $*tr§.igh.t''pbrtibh ' there is obtained! ■■ 

■ t ... 1 1 * -?. , r - • ... ~ i ■ _ ■*. ■ — . . , ... 

, ’ ; V, ' J : * r ' 3 * '* - ^ '-r* -i - S(;.^ 

‘n'ni’s -£I r ; 24-2 •' inch-pounds per 10. ~ 4, radian? ; 

' J* 1 ifC f ~ V .. ,* -- -- -f * - . : . - -a . . I - l . 


> (22) 


r - - •”> r. 


. a- 


, rhg.^'p 0 . 10-2 6 7" pound's ' per ’ 10 ~ 4 ’ radian 

T ■? : J is. r r . ‘ • ■ -■ * J '• ' r ■ 

i- r S = ~ ~V “'*■ ’ — I 

'■ ** -»' ■ •• • -• 

rr- = -0.66158 poundsper 10 ^ 4 inc;h' ' 




>• . i (23) 


L E 

• *F K .r V ; :.; 


s- r dr 


t"t._p’ =~ £j$.,.205,8' pounds per : 1 O'” 4 : ihch ;: •• J- ■ 
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For arc 0D-., the influence. - cbfeff'Ifeifent's-' for tfee - fixed 
end will not fee -needed i nT£iie~'r-elAS:atl on procedure; However, 
it is convenient to have ■ them ; f«r checking phrpo'ses: 


— 4 

nny s 37.649' inch pounds per 10 radian 
rnj =» -7 .'3'524 founds per .-IQ t .radian 
tnj, = 9 . 78 61- \p ound s per 10 • radian 

rr-g, = 1.1T9'416 pounds per :10 4 inch 
try * -1 . 03755 pounds per 1.CT*. inch 
t"ty * 0 . 69 583 9. jpound p - er 10, 4 inch 


\ ( 24) 


Operations Table 


In the,, .operat'd on s tabLe.- -the forces' .and momenta caused 
by the individual ooorations are listed. An individual oper- 
ation consists of a ai splacemen.t,-, or rotation,- of '.one point 
of the frame while all the other end points of the individual 
arcs considered ar e . fe.eTd' r A gi d ly fixed. "The operations table 
can. "be easily established once the influence coefficients are 
known . 


Care must fee exer'cisted in the matter of sign's. The beam 
convention is not suitable for use in the operations table, 
since according to it fhei si gn s of T and N depend upon 
whether T and N act at the right end or the left end of 
the arc. For this reason, the rigid frame convention shown . 
in figure* 6fe will fee used. "'tin this convent foil ■■clockwi se mo- 
ments are positive at either end of the arc. Similarly, a 
tangential force is positive if its directipn corresponds to 
proceeding in a clockwise Sfe'ns'e along the -arc. The radial 
force is positive when it acts toward the center of curvature. 
Moreover, the operations tafele is -Get up fey copsidering the 
effect of each operation upon the constraint rather than upon 
the teas. In other words.the tafele lists the forces exerted 
P’pon thg imaginary supports (or constraints) when one of the 
points is displaced,. ; * ■ • • •* 


In the particular problem at hand, point s A and D 
will certainly remain in the plane pf symmetry when the dis- 
tortions of the frame are caused fey the symmetric 'loading 
shown In figure 5. It follows further from the symmetry 
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that .jthe.- hori zontal ■tq.ngent s &nd D will-re&aiii - ho? i\- 

zontal. On the other hand, the distance between A and D 
will change under the loading'.- It ‘ Suf f i ce s', ther'ef ore y : f or 
describing all possible relative di epldcement s t.o - consider f : t 
thp .beam as rigidly f ixed *at - D and to allowssebfei’ohr' • A *-a 
single ^degree .of freedom , of ,mot-i on;: pamely, thatiofoa Verti-f" 
cal^di splapefnent . ( See for, -instance reference -&‘t) ■■ 1\\ r\ 

o .f I-t r:r --- . ;*/~.5V> ci v'. - • . r :J lc .1^.1?: 2 ”• 

■ In -thei :f i£ st r ow o f - feh§ open a b 1 on s table p(fifttl'«. , 6)i'. , thev’ 

effect of a radial displacement of point 1 A v ti srl-i shed .■ fr.on 
equations (17^. the radial force acting .upon beam.-t:-AB ■■in sr.dc-^ 
tion A, is 3 .34833 .piounds when section .A da di splaced ■••ttfe. "utfklr 
distance (10“ 4 in.-Orln the positive rediAl dirt&cti cui- - .G-dnse’-; 
..quently,, the reaction-’ of. this force; that i s , i r -3 v-3'4833- .pdiradST 
will be-. exefbed -upo-ri .ithe constraint at ■ A . The.'eff feet .of - thid 
..same T jnat io.n, jupou -trhe tfiimed end of the -beam at- cefc-frion r-B< ts':"' 
given by the. ■"In'fluceuc'f-. coefficients f:or . the .fixed ^erfd ,c£o'Ef— u ■’♦a 
tained in equations (22). These influence coefficients :mu»t ■■ 
again be multiplied by, -1 in order to obtain the effect of 
.-.the displacement uporiTthe constraint. The final values are; 

*; .*■'> >1^50 0. Vrbft-H Wf'tftI . J. 

jffjgf s?.j8 .9 22^8’ di. ff <Th^p cured, s-’ jHp « -2.69614 p ounde T^ = 9 67-7 1 pounds 


With these force 
for- the “duff !u e nc 
used in the oper 
sign, as may be 
It may be noted 
at point A cau 
constralht at. A 
force is felt if 
beam shown in fi 


s the change from the bVh-m- VonVedt ibftr f teVed ! ' '* 
e d'bh^f-fY'ofi ent s to the rigid frame cdhveht'i on 
ations table did not 
seen from a cbmparleon 
that a posifclVe"' t'bdial' defl^b’td'bn' erf" Tshe beam 
sed a negative' tkdi'aV forbe' to' act' upVn ■ the 

’ . T- ^ l . L 


entail any changes in 

m b% f i-gb r*eV _ r 6‘a ahfe.. Y 8V'. 

n ■ j i .c* l t. rji. 1. ‘u'L oA ^ , 


. This' co-uld be 
the movable end 
gure 6a is pushed 1 ’ down', ; ** ' 

o si-ii.-'c': i.f*v ‘ ~ 


o'Jcpebt ed ’, * b i h c e an ‘ upward' 
of the elastic cantilever 


The 1 second row of the operations thbl’d’ shows thb’effect 
of a rotation of the section at B. If end section B of beam 
AB ■ is rotated through, a' positive angle) of- 10"" 4 radian, the 
forces and the moment caused at A are given by $he influence 
coefficients for the fixed end ' ( eqdat i o'n’I £ 22 $-)%’■ Oxily^the val- 
ue, of K_a_ is, listed in the operations table , l ' iihde^i^ i s not 
necessary to balance and T^. Because of the symmetry of 

structure and loading, the point symmetr i cally situated. to B 
would always 1 be ; di splaced i'n the same manner as 'B ■ i’s, .and 
its di splacement 'tfculd cause a tensile force and a moment 
which would balance and T^. (See also reference 9.) 

The effect of the rotation of section B of beam BO upon the 
fixed end at G is represented by the influence coefficients 
of equations (23). Of course, all these coefficients are 
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tmM'tJlW 3>y -i-' “before'’* 'Wht e r i n‘g' t'h'em In: t hie toperrai i on ; fflta , b&le‘. 

\ (■*, i\ , - 4 * * » 1 *3 % '' ,**. •* ,% ~ a t ~ . *. .x . ■ r » i ’ j ' * , , f ■ c 

- -Ifl? ha- ‘Effect- of- &■ 'r'ot'sfcl o& of. re'ction -B - Up.ori ..the f ore.esi r_h 
and '-the* mbffi'ent 1 3 - ' "given by 1 comb l'udt.l . o ire: ctf the in'fjkw.-. 
enoe ccef fi'elerft g-‘ f or the movable -.end .«£ ‘.'beam. a' AsB * and" 
It'4nuet not be •'•£ o&giot t'e'&v*' howevfer , 'that in-the calculation... . 
the rigid frame'" 1 cfohvVnt i.on> muat ■ be" iise,d in de.t e raining the 
correct sign of the influence coefficients. According to fig- 
urea- 8a. And' 6b at- the right rend of-' the: -beam' AB- j the'. two' sign 
conventions are' identical; ,c but at the left : ehd‘ of: b.ea-m. Bp: 
only- the signs 'of' the radial’ force are ’the- eaime 'thjoae; ojt; -th,e 
thhgenti hi- 'force hnd the moment are opposed'. «.. (.Of!-., -fi-g’, >7 * : 

Ddfsp'l&c emen't-s and rotations are considered 't.o*!ha“Ve the: rsa-nre- } :: 
proeitiya’ ee'n’s.e' as, the corresponding fore es;- and m.oraent , , .re^Bpftq- 
t-ively. the multiplication by ~1* is ale o-r- carried'. eputt/ ijnr /*. 

order' to obtain- the forces’ and the, moment act in.g: lupxm" the ■■.c.em.- 
stra-int':, the entire... calculation can be given ini.tha. ,f dllqwAnft?; 
li'hesi I- v: .r e -• ’■ .l&.iloi 


l : \ *; ;• . * * - *j r J «■ t *. * ‘ • » p . 

1 . Influence coefficient for 

beam AU •;,! '• o .v-i 7= ~ 

2. Influence: co.e-ff-lci-eifct! for- 

beam ■ BO. - , *- ^ j- . * -i - - 1 e *.* ' r 
-> • • * : f (> % 1 r f 

3. Row (2) multiplied -by -.1, Jb-o ob— . 

tain-; co-e. ff iciest _e c : orr.ejsp ondi-ng 
t-o t p . 0 s.ift iv-e- r ( pt at 4 on, according-. 
t;o- ?'i gldi f ra^n.c? -p op-yjejit Ion , , .j 

rf j V n 1 j. * r 1 . ’ ' * r f S fj * t" ■? ; ; 


i. - ■■■ ”, *■> £ : •: : \ •'! c. i- 1 7 

5. Sum ; of ., r owp-v ( l)--an'd ; - (-d^- 

■ ' ■ : ’ v . : - r - * 


** . , f T ft 1 .* V 

n ■ . • ■ : ' * f »* ■ . X * ;■* *:♦ ? 

j, r ^incp the ^^ilcula 

t a b 1 e.^i , r r i p (i ^ opt . j a 

details a'^.,»ct e giygn,. 


;; r 0 x/ 

' - !■ 


•: f 

- n*. < ft rt 1 


s ■; 

. w ; <•: c- J n ?. 

,;. VJ i n <; {B 

190.634 . 

-pl-9...e73#iS-.. 13.1014 

.• •( f-?.i **■ .♦» »’ ; J 

■ . -f •' 

r , 13,7 ..23 2. 

.. -8 .1026-7 0.,.;, 

. f ^ r* ■ . . * 

. •*. .-I'rC- v -• r . . j- •-! J/ 


' 1 - A v."i". ■•.-«- , " r • 

' •' ^ j 

. ' X CX : 

1 . r -• . • - 1 1 « 

r r. s :- -■* j: 

,.'-137.232 

. . 8 ,10 2.6.?: a •J.-' Ac. 7 

c ’ - • • ' 

*• J l' '* •* ^ . . p* ■ f 

.. i.* r •• 

t ,•> 

l — 

, 13-7.232 - 

.. 8 -.10.26 7, 0 

i ' ; 

• ' r ?'t - ha 

. - 3 27.. 866 •» 

-11,46.97:. u ,13.1014 

« ■ .r 

■ • , i 

* V • f : *; " O 7 

-3 27., 86 6 • 

- 11.4697 -13 .1014 

,1 • 

0 .* ‘ *f/»n v - ; ■ 

:: , ? P 

* ;‘;i ,t V • 

other items 

^iiLtiji? -operations 

the same -prifficdp^e,?) . further 


1 ’ / - sT : - - > •. • iw 


• • r ^ .1 » - * 

: - ‘i I 'L 
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' Solution, of the Problem by . Relaxat i on 

The purpose, of the .relaxations is t ?. £ i nd , such multipli- 
ers for operations 1- to 7. as result in 500 pounds 

and in vanishing values"for all the forces and moments ft's 
to Tq. The value of 500 pounds is stated. instead of 1000 


pounds, since because of the symmetry only one-half of the 
structure;;! s considered. There is no need for balancing the 
forces &hd. moment t s..;- ftj^. , _ T^ . ft. B , and , Tp, since they are au- 
tomatically balanced when the left-hand half of the ring" under- 
goes distortions symmetric to those undertaken with-. -the,, .fight- 
hend half. Moreover, because of the conditions of equilibrium 
of the external forced, ' a"’-r e eidual.. force R B - 600 must be 


obtained automatically at point D when all the other resid- 
ual s’’are' r ;ta ought down to. negligibly small quant i.t ji.es . This 
force' is then balanced by" tone- ha If of the ICOO-p’ou-nd external 
force at D. ' J .• > 1 

' • ; ■ •: f* > ; •; . • .1 ‘ * ’ ' ‘ 

In working out the present problem it was thought advan- 
tageous to begin with a rough guess as to the final shape of 
distortions of the- •frames '. -I.t- wJa-S. c;<^veni.ent ;t,o. consider sec- 
tion D as rigidly .fixed., apd to assume reasonable amounts of 
radial and tangential di splac.ement s f dr ' e~acfr points -A, .y-B-,.;. 
and 0. The amounts chosen' are jgiven below together, with'_the /' 
forces and moments caused by their simultaneous occurrence. 

A 8 may be seen,. ‘the radial and tangential .displacements listed ■ 
correspond to displacements of points B .. "and C upward- and- 
to the left, and :to an upward di splacement !o£ point ' A. The 
rotation of — B was -aseuraed to be. cl'oc.kwi,se.,‘;3.k&t^pf 0 
counterclockwise. The forces and-- moments .w'ere obtained. by 
adding np..corr esponding -values of - the oper&ti ohs . table after 
the- rows.' were mult ipl led" by the • ‘fae-Vors chosen .- .The diet or- ' . 

■tions arjes ‘ ' * • “ - • 

- ■ \ 

v^ as -10*0 units' w B = 3 units -v-g = 50 units ■ u-g- = -T!T0~unl t s 


w 0 =• -0 ."5 units Vq = 50 units Uq'= t- 31 units 


The for.ces and moments resulting -from; the- di st or 1 1 ons are: 


'S« 


N 


B. , 


R ; B . 1 


T- 


B 


ft. 




rnET- -• 


■ 0 - 


98.839 ,.—145-. 5.7. 26"V132 -134. '9' -368 76 ~1-6 T 77— -11 . 57 
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The set of distortions listed can "be considered as a 
"group displacement." It will tie denoted as. operation [ 8] 

In t-he cour.se of the relaxations it is advantageous to make- 
use of ot h'eSr combinations of the elementary operations,. Two 
such combinations were used in the r elaxat i oh’ table (table 7). 
They are listed below. 


[9] ^ ^ ^ . „ o _ ■ ^ 

~100(u B +u c ) -396.771 1310114 -343.52 47$ . 28 ' 524 . 667 88.929 87.75 

Operation [10] ' 

[9] - 120 v A 5.03 239.5 -19. 98 -0.85 524.667 88.929 87.75 

_ /; ; * : 

In op.erat i on [.9] the shape and* the length of the straight 
beam r( beam* 3C ) do not change; -Operation s of this kind are . .. 
denoted as "block displacements." In a block displacement 
part of the structure undergoes a rigid body displacement. 


Table 7 . i-:- S\1 r 8 1 Relaxation Table 


Operation 

• S A : 

^b"‘ 

*b : * 

' - »n i v l V . 

. B - •• 

N ' 

- C •« r • 

RESH 

■SB 

0.3 [9] 

98,839 
'.-,119.0.3 . 

-145.67,, 
. 393.94 

■ 26 : . 132 : 
■ r 103. 06 

-134*9 h 
14:2* 5jBc 

-•368.6 ' r . 
'■ . 1'57 . 4 

*-16.77” 

. 2G.'68-’- 

MW.!*k 

-25 [3] 

■f -2o.ao' 

67.40 

'247.37-' 
-286.74 * 

-76.93 
"1 00,25 

> • r 


-211 . 2 
202,6 .. 

■ 9 • 91 n , 
-1G*. 64 

37.90 ? 
■ ’0 . 

-2 [4] 

47.21 , 
. ■ -7.94 

1 


: .>78.20 

-8.6 • ■ 
0 " 

* 


■ 

mi 

'• 39-. 27 " ' 
0.50 

-13.17 " 
23.95 

_ r. — 

" '16.87 
-2.00 

-16.29 

-0.09 

-8.6 

62.5 

-6.63. 

8.89 

-14. 51 
8.72 

0.15 [5] 


. 10.78; . 
-9.19 

■ 1 


■' 43,9"’ 
-43.3 

mm 

-5.79 

-0.79 

5 [1] 

39.77 

-16.74 

1.5? ft 
44.6 

18.85 

-13.48 

-16.38' 

19.84 

D.6" 

j^g| 

IggH 

y55555l 

23.03 

i’;25 

, 46.2 c. 

' -45.9 

... -0, 25-. '• 
1.61 " 

.'3/43:- 

-1.83 

j 

0.6- 

-8.57 

■ 1.82 
1.13 

-6.68 

0 


: 24.28 

0.3 

1.71 

1.86 

0.23 


l'" ''''''ll 

■ml 

-6.58 

0.15 


■ 

■ ■ 



' 

mmmm 

3.03 

0.08 

-6.43 

4.06 


24.28 

2.0 

mm 

-2.3 

0.6 

mm 

-2.37 
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, Instead of balancing the largest residual force at each 

stage, it was found more convenient to adopt a particular pro- 
cedure. After . operation [8] was comple.ted, an attempt was' 

» - ■ made to reduc.e all the . residuals hut E-A to negligibly small 

quantities, while keeping E^ as large as possible. . The pur- 
pose was to establish a distorted shape corresponding to any 
finite radial., force at A and negligibly small residuals at 
•points B arid. 0. Division of 500 pounds by the radial force 
■ at A so obtained yields the factor by which the complete set 
of operations undertaken must, be multiplied in order to obtain 
the final displacements of the frame under the load specified. 


Table 8.- Check Table 


Operation 

R. 

A 

- 

r b 

T 

B 

R c ’ 

R c 

T c 

-107 v A 

358.271 

-954.671 

288.487 

-424.545 




2. Id w_ 

O 

19.093 

-701.633 

24. 545 ■ 

-28.037 

-131.058 . 

17.340 

0 

25 v 

B 

-67.404 

286.742 

-100.248 

85.880 

-202.567 

16.540 

0 

-72 u B 

-285.675 

943.301 

r 247. 334 

. 2228,. 875 

0 

o • ■' 

-1886.818 

-0.378 Wq 


23.149 

o • 063 

0 ' 

109.003 

1.117 

1.983 

50 ^0.. 

... 

405.134 

33.079,’ 

O' 

-147.811 

-95.103, 

. -44.465 

-7i;i5 Ug 


0 

• 0 . 

-iS6h.543 

37-3.301' 

63.273 

1926.977 

Sum 

24.286 

2.022 

1.592"- 

-2,369- 

0.868 

.3.167 

-2.322, 


In' table 7 “the objective was attained to an accuracy,, of 
13 percent in 'Eine .operations following operation [83. Before 
the relaxation was continued, a check table (table 8) was then 
established. Some slight differences were found in the re"-' 
suits , and the 'residuals of the check table were further re- 
laxed in 11 more operations Cin table" 9). 


■ - . : f. / 
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Table . 9-.- Second. Relaxation^ Table 



After these the greatest residual force at points B and 0 
was 2.5 percent of the radial force obtained at point A. A 
second check, details of which are not given in this report, 
resulted in the residuals listed in the last row of table 9 . 
Since the radial force at point A had the value 
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- . f -v =, ,1 f -4 7 1 p ojin d s , 


*j r> < 


all tiie 


le .operations Jia^l ■ 1 t.o" > b_e*' multiplied" ^by the "fact or 

, ! ' *“’1 *“ *’ A ii / *’ S’ -■ ' \ * » V -.-v 

' •' •' ! - ni.. \ r i .!* *j- . . 

U.* -,6A0-/a 9.471,=' 25.‘67§V; ; 


,3 <C .T* i. r . 


in order to obtain the final deflected shape of the frame. 

The f in.al . def le ct.i on. s are given below: 

I J :! S -J C : I ’ ■ . . 

... ' ' <«. ' -\- 

VA = -0 .'27 ! 6'4 ’HSSli »B * 0 .-000774 radian. v-B *<0.06420 inch 

u® »,:-0.^|49S iScii' *o * -0. 001066'? ‘radian v c = 0.1.3353 inch)(25) 

■ ‘ * •:' r c i ~s n - ~c- 

UQ =..70 .* jaj3 0 |„ inch 

&ay" be ' seen,-. that the.-.first guess for" the' 'deflected 
shape, as given in operation [8], differed considerably from 
the:, final- shape'.'" ;.Ney er^he^Le ss , only a comparatively .small num- 
ber of- relaxations". was:.' required for the "solution ,o-f the prob- 
lem. The- deflected^. -shape is shown in figure 8. Next, the 
'bending moment distribution is calculated. For 'this purpose, 
the tangential fpr.ce. apd the bending moment at point A are 
needed. these are caused by the displacements (and rotations) 
of points A and.. 3, and can be calculated with the aid of 
the influence coefficients given' in equations (17) and (22). 

The calculation follows: 


Caused by: va=10 4 va»*-0.2704 wb=10 4 wb*0. 005774. v;g =10 




N, 


-19.5724 52924 

-3.4352 9288.78 


16.7165 965 

21 .7972 1258 .-57 


-8.92216 
-3 .96771 


Caused by: v^* 0,06420 ub* 10“ 4 up=-0. 18495 Sum 

*A -5728 '21.7972 -40314 7847 inch-pounds 

T A. .,-2547.27 4.31595 -7982.35 17.73 pounds 

Since influence coefficients were used in the calculation 
of the moment and the tensile force at' a, the .values ob- 
tained represent the action upon the curved beam; and the 
signs are in accordance with the beam convention. It must be 
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noted, however, that the signs of the *di splacement 8 used in 
the ^receding calculations also had to he taken according to 
the beam convention. In. the present ca$e, this requirement 
did not cause any change in the signs of the displacements- 
obtained from the relaxations, since the displacements and 
the rotation at the right end of the curved beam, and the ra- 
dial displacement at the left end have the same sign for the 
two conventions. 


With the values of and T^ so obtained, and with 

Rjl - '■500 pounds, the bending m'oment.s were calculated,, for sev- 
eral points along the median line of the ring. The moment 
'distribution is shown in figure 9. . The- handing, snom-euti-, the_ 
radial force, and the tangential force at point D were also 
calculated from the influence coefficients in .order ,fc o. obtain 
a check of the accuracy of the calculations. Care must' be ' ' J/ 
taken,, t-p. ,.t,r±an..sf,Qrjn ..the signs jsf t,he displacements at point C 

"fcO 1 C OH f O f * r ^ ^ ^ +■ Vi a Vi! m A «« Tr A’w f m* ^ ITCVh a /I A ii-4* 4 /ivi a i 1 1 r\ u • 


0 au s e d* 


inf orm, w.it^ beam convention*. '; The ‘compaltationB follow: 

''.•'.c a r? v i * 

fd* by':' w 0 - vvltcil* .s'w e * d..001 , Q,6^7 '! - 1 V 0 *'0.l'&353 

— i-, .... . .. - 1 C . .. — 


— r- — mrt — - — — — — ~" 

US7.649 . f.;.c .'40,3 ^ •• c -9818 

:~\f. -Jri o,'| ,i •• r,c: rh.» 

h - *,£ ?}3 524' •••'.' 1% ! ' f i -,1.11941? .[ 1495,' : 

; r -' ' ■■■': *•.. ir r -i 

Tp' ' § r?88jro i - ■■■ * ft{104v39 ;-l/63?65. r .'-13.85:44 ** 

r * ■ 'c / " - 'i ■ ■ .. - 

Caused by: uq = 10*'^''' up -^O.I83G9 v v Sum .s - 


V 

V 


f : . 


9 . 7 8 6'1 v 'i '17917 

-1.03755''' -19 00 

0.695839 1374.01 


v. 8-501 inch-pounds 
■ ^483 p p.unds . 
-7.04 pounds 


These values show that the accuracy is sufficient for en- 
gineering purposes. The radial force obtained deviates about’ 
3.4 percent from the actual -500 pounds. The bending moment is 
shown as 9088 inch-pounds in figure 9. The deviation of the 
average of the two values from any one of the values given is 
3.3 percent. ■ The; forces Tp aftcl.-.,, ,Ta should add up t o zero. 

In reality, their sum i e 1 Q.-SS 'p ou.nrd s,- which is about 3.1 per- 
cent of 500' pounds ■”■>£. 'that is - , .of one-half the applied load'. 
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The Growing Unit Procedure 

In the growing unit procedure the unite of the ring are 
combined gradually into units of increasing size until finally 
the entire ring becomes a single unit. In each step of the 
procedure two or three simultaneous linear equations have to 
be solved. The growing unit procedure has been developed, 
since in many problems the relaxation procedure is too slowly 
convergent. This is true particularly when the shape or the 
loading of the frame does not allow a good guess regarding 
the shape of distortions. 

When this procedure is applied to the egg-shape frame 
loaded by the equal and opposite forces at A and D in fig- 
ure 5, it can again be assumed that section D is fixed in the 
plane, and section A is only free to move in a vertical direc- 
tion. The first step in the procedure is then to combine 
beams AB and BO into one unit AG. This can be accom- 
plished by determining ub , vb , and wb so as to cause N B > 
Rg, and T B to vanish when in turn a unit radial displace- 
ment, a unit tangential displacement, and a unit rotation are 
undertaken individually at points A and C. 

In the present problem section A moves only radially. 

The forces and the moment acting at B when v A = 10 -4 inch 
are given in the second, third, and fourth columns of the 
first row of the operations table (table 6). The forces and 
the moment caused by the unknown displacements and rotation 
at B can be obtained from the seeond, third, and fourth rows 
of the operations table. The requirements of the vanishing of 
the forces and the moment at B can then be written in the 
form: 


-327.866 W B + 11.4697 v B - 13.1014 u B = -8.92216 
11.4697 w B - 4.00991 v B + 3.4362 u B = 2.69614 
-13.1014 W B + 3.4352 v B - 30.9566 u B = -3.96771 


Use) 

- 


Solution of the equations gives: 


w-g = 0.003449289 x 10 4 radian 
v B = -0.612145 X 10“ 4 inch 
u B * 0.05878271 X 10 -4 inch 
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The displacements foundP ''gi valise 1 '-t o fox ces '..and: a moment at 
0, which can he calculated from the values in the second, 


third.;’ 1 : 
at i on s 

and ,: fQurt'h 
table: ■' ,1 

rows 

arid the* last three ‘columns of the 

• ic - ■ • * - f 

oper- 



- 

'' " * - ^ 

4 . 74-^7 6 7 n-ch-p oun d s 


f 



*c 3 

-0.377034 3 pound 

> - 

( 27) ‘ 


r >o 


T c ss 1 . 5404'48'^p 6und s • 


In order, to he able to balance out the forces and the 
moment at 0 - as given in equat i ons , ( 2? ) , section C must be 
di splaced -apd- rotated . . Duping these motions; however, no'uri- 
balance s . mu et -appear in section B. In the next Bt ep' - of , the • 
calculations, therefore, section C will be"in turn rotat ed, •- 
di splacedc radially, s and displaced tangent ially, and' in each 
cas? Value s of t - w B , U B will he' calculated 'so as to 

cause- Hg, Rp.f.and T B to vanish. When section, C is ro- ”* 2 

— 4 * " i " * • 

tated through ^.' positive angle of 10 radian, the operations.,., 
table gives: 


r ■? 37-.- 8 66 . w B + 1.1 '1,4697" v^' 1 ' 13 .1014 u B «* 
.11.4697 w B - 4 ,ji0991 v B + 3.4353 u B » 

- • J / t „ r , 

.-,'13,10.14 w B + 3 . 4$ 5 3 v B - 3 0/9566 u B a 

• ‘ ' * r * * * •' . l . 


61.242 



8.10267 h (38) 



w J .* 

.1 a 


It may be seen * that • thb left-hand side terms of equa.Vi ops,, ,( 28^ i .' ; " 
are identical with those of equations (26). Solution of' eqyip-‘ ,, 
tions (26;) yields ’ 


W B = 

.. ;. e 


• -^4 * , ; JL • s "\ 

-0.28372447 x 10 radian 


v B = -3.0160675 x 10 4 inch 


■' o 


} 


(29) 


u B = -0.2146168 x 10 4 inch 


j 


These displacements cause forces and a moment to appear in 
section 0. Their values may.be calculated with the aid of the 
operations table as ‘hfe-fore: 
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■ t -.p-V-' ■ • • 

ufjj-- _ 41.814055 i nch-o ound s 


•ii 


E.q J = -4.2943 0 pounds 
Tg = -5.624205 pounds 


. > 


(30) 


In addit i on *t6r'ce s and a moment are caused at C "by the unit 
rotation of 0 itself (while sections B and D are held fixed). 
These quantities may;be. taken directly from the operations 
table : . .. , . _ 

-• v . f X A ~ 


N = -288.3 6? inch-pounds ,. 
R 0 = -2.95622 pounds 
T^- = • j- 5 . 2466? pounds 


vo 


Consequently , a unit rotation at C, together with the mo- 
tions at B required for Ng , Rg. f , Tg t.o vanish, gives rise 
to the following- force'®, and- moment 'at ‘ C: 


o - 


t.) 


R c = ’-246.55294 inch-p.Qup.ds ;q 


; r • 


R~ = -7.25052 pounds > (31) 

; ’ : = ’ -1(5 8 f 08? 5 p ounds 


j .a -.-ion : ” . 

The effect of a unitradial and a tangential displacement 
of C, . respectively, upon t.he'i'fpr cV& r and the moment at C 
can be' calculated', in an analogous manner . Again, the left- 
hand side members- of the equabioft.® a&^‘” identical with those 
of equations (26) land (28) 
the values 


)iand (28). The ri ght-hand- side members hg.ve.__ 

a ' r i ! * -i * ‘ **••**" * •* 'I “ •* ’ 


-8.10267, -0 v 6 6158,.,a-ndr y :.0 -,- 

no •:>’ ■ .. _■ . , 


when 


when C und 
of the equat 


und.ergoe.sf ai uh.'4t : --=rad'ial displacemq-nt and ' ;r 

- ' : " ! _ ’ . . . _ ^ _ 

0, 0, and -26.2058 

. . r . .. : vo' O'’' 

q~ --- v . — 

ergoes a unit tangential displacement . Solution 
ions gives o: ; ;r - - ^ 


ii. ox: 
i T.U 
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w B = 0 03369.378 x 10 4 ; ,radian 


u B = 0.2753198 X 10 . 4 inch 


(32) 


when 


v B = 0.01629283 X 10 inch 
undergoes a unit radial displacement, and 


w B = -0;01'.639'8'4 ! ff x 10 4 radian 

u B » 0.7726881 X 10 -4 inch 

" 7 - ri , • i ■■ : 


(33) 




0.9366781 x 10*" 4 inch 


when C undergoes a unit tangential displacement. 

The forces and moments caused at C by these displace- 


ments are' 


(V ! - f Si 


N c a- -4.29430 inch-pounds. 


T c — 0.426967 pound 
R c = 0.4551558 pound 


0 


(34) 


when 0 undergoes a unit; padip.1 displacement , and 
:• <• ,,^.Q = -.5.6 24305 inch-pounds 

■ ' • * * * •* ' l * * • -I ■ * t 

"V' : -S(j‘ * '0\4"#&i67 pound 

i • <■ • . <:/■■■: ;i . . - > -ir.r: ^ ■ * r „ ; * 

. . * '24,. r 546'4- pound's ’* ; 

■ r s « ; I ■ f : .. i /• r * **• • • \ ; ?■ 


'(36) 


when 0 undergoes a unit tangential displacement. 

■ ■ , ••• .0 . 

The effect of the motion of 0 upon the forces and the 
moment at 0 (considering ~ JB. and -33 ;■ &*. r-igid-ly. -fixed )o i «■ . ,■ 
again taken directly from the operations table; 


*0 = "v?. 9,56,22 inch-pounds 

- ’ ’ ? . . ' * ■* • ' rr f ;j 

Eq = -1.90206 pounds r ■ ' ? 
Tq * -0.88929 pound 
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when C undergoes a unit rddial df S£I&c§fteftit' , : And. ' 


t .Hg. ; = ±S:.%4&67 inchvpounds 
Rn. = -Q. 88929 pound 

*5.'. V '-l *■ " .. . ^ 

Tq a -27.0833 pounds 

* ‘ .1 . 

when C undergoes a unit tangential displacement. 


Altogether, the. forces 4 n< ^ moments at 0 . become, wtyen . 
0 is displaced radially and 'tangentially , respectively , 
through a positive unit distance: . ~ “ 


Nq = -7.25052 inch-pounds 

Ra ■= -3! .44S8942 pounds 
• ■ ' .I - 

Tq s -0.462323 pound 

• ‘ !- J 

« • . . . t t : 

in the case .of a radial di spl secernent K .and 


(36) 



. -i t tr r- JT q = --1 0 . 8 7'08 7!5 lunch-pounds 

\ .<?'• 5 cnr * ’ 'll '■■■?■ i ! X s s :/*r * ■ <.' 

:mi -•-'■ST'c a. hQ 1 . 4^23:2^3'' ^.dlindh c . 4 »» A 


, ^ 

(37) 


Tq = -2.5369 pounds 

b l i Of' • - * ?■ s. > x c! \ -A .1 o 1 -1 v i* . 


in the' case^.of a , ' t ji p g g n.t i ad d i j p Ija perne ht . r ', H 

Equations (31), (36), and (37) give- ‘t'fe.e^ fb'rfee-te and .the 
moment caused in section 0 by the three possible motions of 
section O '-*, t ’■ T ire 'sb- ^ a lb. eV cor ire spoftQ/ 'i V* bu ch" &i Vb orYfc ions of 
the bin g-- -k's d'p* n'bV' 'ent ail uhbh'labc'e-&'-'Tb'r c-e%- Vt^lab'toent s any- 


where along the ring except in the end sections A and D. 

The values of • thq r mptiohs.’- g.tx Gi ? jn&ss.t. -hows; b_ev det ermined in 
such manner as to balance out the forces and the moment caused 
at 0. by the.'unit radiaf’di splacgmept " of- A,% These unbal- 
anced forces are listed in equations (27). "The requirements 
that Kq, Rq i and- n Tg • taust ,vq.ni sh' can be stated in the form 

of the following three simultaneous equations; 


•i n 


. " , T 246 . 5p94:vf^7” f§0§£**J c | Q, 1 70 j * S* *4 , 7.4 8 7 6 ? 

-7.25052 :> ^ 0 '^ 1 ;^46894§ T: t c n I'''' 6'. 162323 V uq" = ‘ +0*39*6343^(38: 
-10.870875 w 0 - 0.462323 v 0 rr 2.5369 u c = -1.540448 
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The solution of these equations ie 


a '■ :■ - i c » ■ :: .• t ' ■' - •• ' 

Wq = 0.003849891 X 10“ 4 radian 


v c = -0.4976019 x 10' 

a 0: 


>-4 


inch 


(39) 


u G = 0.6814022.x 10" 4 in.ch 


J 


The se* &-r e - the actual motions -6-f. section C when 'A is - -" 
displaced radially through a positive quit distance, and— when 
the motions- of 1 ’ BJ.iara ' sti-eh as’-t'o "r'esult in vani shihg* f or ce s 

and , momeht •: at • . B . . ■ The total' motions of B can he - obtained 

by adding up the values corresponding to the following- fbur 
items: 

1. Motions of B when A ie displaced radially; given 

-. by the solution^. of equations ( 26 ) 

2. Motions of B when ,.C....is rotated; these are the 

values listed ‘in equation's (29) miltiplied by the 
factor 0.003849891 taken from equations (39) 

J ; * - > 

* t ' J 

3. Motions of B when C is displaced radially; these 

are the. values- Dusted, in- eqqat 1 on s (32) multiplied 
by the factor -U.4976'0r9 taken from equations (39); 
and . . 


4. Motions of B when 0. . i-s displaced tangentially; 

these are the value s‘ I'i'st ed in equations (33) mul- 
tiplied by the factor 0.6814022 taken from equa- 
tions (39) .• •- s ■> • <*. . o ! 


ment s 


■■ " i\ 


Add it i'-o'n' ’of the individual ? rotationp^ rh^a r l‘ di’splace- 
, 4h'd' tangential displacements, respectively, "yields,: 


r. f. . J» £ C * ■ .1 —.A ■-* J O ^ ji. 

~0^02149^4 x .%0 r^dianjj 
-y B ». -0.4g3^?45 •xa.O.tj^h.. a ,A._ 


o r. v 

> (40) 


i r. . 

n 


> s: 


. -ti i':'.' ;• fail..--. ■ -■-•■ 

• 1 tLj = 0 . 68 8^055. ,vX IQ ingh -« r ; 

• ■ ‘ = J - - h 


It is nosgiblq -ft-pw- to det ermine ' tha^radial force At’ A c 
corresponding to the distorted shape of the ring just calcu- 
lated , r^-Th.i.s ;; f vd'SDp'ftnd-s u-pon- t.'he Motion a ’''of point# '"'•'A 




'»!*. r» jf 


- <s V 


i - 


,c- 


i r* 3 r. - 


Or 


riV 
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and B. According to the first entry in the .ope'rat i ons- table 
(table - ' 6) the radial force at A caused by a unit radial . , 
displacement of A 1b, provided B is held fixed, 


R, . = -3.34833 
AA 


pounds 


If now A is fixed in its displaced position and B is 
moved into the position indicated' by equations (40) , the ra- 
dial force at A caused by these motions can be calculated 
from the entries in the first column of- the operations table 
( table 6 ) : 


R aB = -0 .0214944 X 8.92216 + 0.234245 X 2.69614 

+ 0,6881035 X 3.96771 - 3 . 169978 'pounds 

The total radial force' at A - is the sum of the pr seeding, .t wo ^ 
forces; 

S-A .= * AA .+. H AB =- -0.178352 pound (41) 

In order to balance out the 500-pound external force carried 
by one-half of the ring, all the distortions of the ring and 
force R a in equati on.. £41) must be multiplied by the factor 


f '=>■ 1B00/0, 178352 = -2803,4448 . ..... (42). 

Multiplication of the values given in equations (40) and 
(39) yields the final ,di splasement s and rotations; 


II 

> 

0.006025836 radian 

V B = 

U B = 

-0.19290602 inch . 

W C =* 

v 0 * 

0.1395000 inch 

U C = 


0,065669293 : inch 
-0 .JOOl'079 29'6’ -radian 
-0,. 191027.345 inch 


\ (43) 


At the same time the radial di splacejnent- o f point A is 
naturally 


. v A =.-0,28 03444 8 inch.- _• 

With. these final values 'of the displacements and rbta- 
tions the f orce s .and . moment s at point s 'A > -and-' : I> can be ; ? 
calculated from the influence coefficients. 
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These calculations were, • carp i ed out in the same manner 
as shown ip .conne'ction' wi’th the- evolution by relaxations. The 
re feu-1 t : s, 'are ’ ■ . ' * ■ - - - 

.. b. ■■ •' 

* H. = 7970.2 inch-pounds 

A • - 

T^ = 13.24 pounds 

. •„ •;} *. • ’ ' S ■ 8868- inch-pounds 

i > ' " 'S_ a -49 9 :8 "y ounds 

^ ! 1 » • 

T^ a -12.52 pounds 

The accuracy of these re suit s, i s . bet t er than that ob- 
tained by Relaxations. If Np is calculated as the sum of 

Nj ' and 70 T A , the result: namely, 8898 inch-pounds, dif- 
fers, from. , N-q. = 8868 inch-pounds by about 0.3 percent. The 
value of Rjj is practically equal to one-half the external 
load. The sum of .. ¥ A ; -add ' ' TjJ deviates from zero by Less 
than O.g percent of one-half the applied load. , 

r . . %s ^ ' . • * ’ , . ... .. . ? , * i : 4 ‘ • ' ■ . .. c . ,, t <: * ’ " 

■ ■ . . * Other Solutions of', the ffobiem 

As was mentioned in connecti-oh ’.With' -the ^pr obi em solved 
in the first section af'-'this paper, namely, the pjrob.lem of 
the Arc consisting of two quarter, pjrcles, ’ the "operations 
table of any, problem, can •. be i nt erpr et ed as a system of simul- 
taneous • linekr'" equati ons . The first two equations obtained 
from the operations table of the present problem- (.t-able 6) 
are listed below: 

-3.34833 v A + 8.92216 w B - 2.69614 v B + 3.96771 u B = 500] 
8.92216 v A - 327 .,86 6 WB + 11.4697 v B - 13.1014 u B 

. - 61 , 242 w' 0 + 8.10267 v c 

In the matrix cal cuius, .matho'd-s are developed which are help- 
ful in solving large numbers of simultaneous linear equations. 
In the appendix such methods are explained The. eqttat i ops . .. , 
repre senti.ng , the . eg-g-sh^pe’. ring 1 oaded" Viljih -,t wo equal 


>(44) 


= 0 
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and opposite forces were" also solved by the matrix calculus. 

The Jesuits obtained, are. a.s, f.pllpws.: 

-r. t r : I r •. -r. u ■« •. _ _ ■ „ _ _ _ 


; . .V, 


r 1 ■ - r .• ,i - 1 •’ •• • . 

= -0.280,349 inch 


v B = 0 . 0656-7005 inch , 

«j • «• . i . f 

v 0 te. 0.139 5P3-- Anph*. . 

= T 0. 192910' inch' 


!*•> i 


-0". 006 025d07^raJ ian 

■ V , ?!"• r ; . • !- • V’ *i“ 

.. -w«- =•■ “0; 00-107944 • radi ah ■{ ■ . 

. 0 ;:J • k- fUs) 

•* .» Uri s * r *-© 1 9 1 6Z1 inch ‘ r ’ 


I. 


• U B 
■ l-’ 


.■ if, .j 
•i . . 


f-r * . 

* 1 


... J;;; 


-rr •* ^ !"• 


It may be seen that these result s ar.e’’prp.ctj.caiiy the 
same as those given in' equations (4&). It follows from this 
agreement of t he . di splacement s that , t,he .bending moments; the 
axial forced, h£d shear forces must also agree. The three 
- ost int er e st in JV^hlue s are given below: 


m 


N 


A ;.-='.7970.2 inch-pounds 


"A 


N 


D 


-■12 .-57 : . pound's ' 

• - 

- , 8S44' Inch-pounds 

i -■ V ' : S C ’ • 


> 


(46) 


■ . j • * f f. • ‘ r-A *.“•>} 

The solutions by relaxation, the procedure of the grow- 
ing unit, and the matrix calculus were all based on the oper- 
ations table. An independent. che.p,k„ wa.s. p^t^inped by v , solving 
the problem' by the semi-graphip . p?^o-'cedur'e_ de.scr ibed- in. r,e f 
erence 7. ' . ’ . . : 7. " ~ ?"•*. i 


-r V 


In this procedure four . imaginary pins. were ln^erJ;e^L‘ ipi.-- .- 
the ring, and the bending moments at the four pins were con- 
sidered as the statically indeterminate quantities, The. 
unknowns were combined 1 int o suitable symmetric ' dhfi' ant i sym- 
metric groups, and their magnitudes were calculated with the 
aid of the minimum strain energy principled details of Che 
calculations are not given here, only the valijps. obtained for 
•^A* A » and are Quoted: 1 ' * ‘ • “ r - 

N. = 7840 inch-pounds T, = ' 8 . 86 -pound s ‘-N^ =k i *84 6(3- inSK— pound s (47 
A A -,*:r » v «v .* * i t ■ * •• .pc 


The deviation of the value 
in equations (46) is about 
tion is about 4.3 percent. 


of N a in equations (47) from that 
1^6' p er c'erif; ~wtth N^ the devia- 


■ -* f K f 

r - • A ». - - 
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The agreement’ between the results of the graphic and 
the matrix analysis is astonishingly good taking into consid- 
eration that in the graphic analysis extensions and shearing 
deformations were neglected, while the matrix calculation was 
based on the operations table which was established with due 
regard, to these types of distortions. For the purpose of 
gaining further insight into the effects of, extensional and 
shearing deformations, the influence coefficients were recal- 
culated neglecting the effect of these deformations. This 
amounts to setting l/Y equal to zero in ail the etjuations 
of reference 5. The operations table calculated with these 
values of the influence coefficients is table- 10. 


In ,thi.s operations' table section 0 is displaced tangen- 
tially whenever section B is displaced tangent ially , and. the 
two di splacement 6 ar e always equal in magnitude and sense. 
This is a consequence of the assumption that the ring, in 
this cane the straight nortion BO of the ring, does not 
change its length during the deformations. Because of- thiB 
rigidity of the straight bar in its o-v/n direction, it cannot 
be determined whether a tangential reaction would be taken 
up by the constraint at B, or by that at 0. Consequently, 
the tangential force Tjjc is listed in the operations table, 

where the double subscript BO indicates that the reaction 
may he a t B or 0 . 

* ’ ' * ' ' '.*• it " • r 

It may be £een -that operations table 10 differs greatly 
from operations table 6.-- Nevertheless, the matrix solution 
of the system of linear ' equations comprising the opexrat i o-n s . 
table was found to be very similar to that given earlier in 
equations (45). The present solution is as follows: 


v A = -0.25186 inch' “ 
V B 0.068543 inch - 

Wq « 0 . 0009 1890 ; radian 





'0,0068334 radian 
u 0 » *-0. 17491 inch)(48) 


Vq a o; 13340 incb„. 



The moments and forces were again calculated as before. 
Some Indicative values are as follows: >. 


N i 


T 


A 


N, 


b 7898.9 
» 1‘3 .134 
* 8319.3 


inch-pound s 
pounds 
inch-pound s 



(49) 
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The .difference ,] 3 etweeii,,Qs^re 8 p.ondiiig values 'in equations (49) 
ana ’ ^4 a v ^bout 0.9 percent; ih.ithe case of’ N A ,' and about 

6.3 petfcemts in the case 1 'of - iN-tf. fr,l jC comparison o’f all the five 

*r .*,*.** •» .. ? * O * ‘t S- • ’ 2 .. ... 

solutions of the p,r oblem _i s givep below: *. 

-r v . » •* j • *1 - * ** ’ ,r r ‘ ’ ' * ‘ ’. - _• 

. ; ■' .* j’ 1 • ‘ '*■* - r 9 1 '** " . . 


Relaxation Growing unit- Matrix Graphic Second matrix 

: j(ihexten signal deformations) 


N A (in.-lb) 

7SU-7 

7970.2. ■ /• * ; 

7970.2: 

■ 7840 

■7.89s. 9" 

T A (lb) 

17-73 

13.24^' f 



•13.134 

Mp(in.-lb) 

8501 

8L868:'. 

1 

' 8844 ' " ' 

4 r . 

s46o : , ; : 

,8819.3 

f ■ 1' 

%)(lb) 

-483 

-499. s 

-499. 8"' ’ 

-500 

”■ ''-560.0 

•J . . 

T D (lb) 

-7.04 

-12.52: -• 

-12.51 

-8.86 

-13.152 

v» ■. 

. . . 

1 



1 - * 1 ■■ I • *■ J. 

' M U<> ‘.G*cr 


The foment di stribution. *i-s. shown in figure 9 for the' so- 
lutions" 'by relaxation, ma^trices, (including the effect of 
shbaVi-rfg- and eitensional deformations) , and graphic analysis. 
The '-chives are very taujch' alike. - • ' 

•:1 ■ : * 

■ » ««*> . 

,i ILLU STRATI VS EXAMPLES 

. •; ■ -• ■ , t ; • ... .......... . 

r • The f Rr obi em' of . Li e t r ibut ed Loads 

r v» . ^ ’ • • • i l '. . - t . 

The first example t aken,,up & in ; ,£his chapter dealB with 
the calculation of the bending '.moment distribution in a cir- 
cular ring frame of a fuselage. -the .loading of which consists 
of two equal and opposite f or cea - wh.i oh ja-dd’ up tb”a moment in ~ ' 
the plane' of the ring, and of, ..the balancing ■‘shear 'flow trans- 
mitt ed- t-b^ 'the ring from the s^e.et .covering, of ■ -the' -fuselage . " ~ ' 
(See yfl'g.v '10 . ) It is customary to assume that the"shear flo.w" 
can be_ .csa'l culat ed from the formula q =* T/2A , ’ - wherfe q v ‘‘is’ ' ' 

th ! e--she&r'-' flow in pounds per inch, T the torque "in inch- 5 " 
pounds'.', e.n-d A the area of the circle in square inches.* " ■ 
Recent investigations Dy Vfignot, Combs, and Ensrud (reference 
10), and by Hoff (reference ll) showed that in sections of 
the fuselage adjacent to rings loaded with concentrated for- 
ces the shear flow differs considerably from that predicted 
by classical the.ory, ^unlqss the ring has an unusually large 
bending rigidit-y. Because of its basic assumptions the clas- 
sical theory holds orhly at a distance from the load applica- 
tion which .is of the order of magnitude of the diameter of 
the fuselage. 
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■ ’-* 5 ** "T '** t 

It is not the purpose of the ~ pr ppent . r .eport to investi- 
gate the ‘ shear flow in the thin wa.li_ of a monocoque fuselage. 
The report deals only with' the c&icul : ht!ion~ of • the' bending 
moments in rings the loading of ■■ whi-.ch i s.-.appoi fi ed . It will 
be assumed, therefore, that the shear flow in the present ex- 
ample is constant and has the value predicted by classical 
theory. 

The procedures developed in the preceding chapters can 
be easily extended to include the effects of distributed 
loads. The basic idea in these procedures was to break up 
the ring into a number of units, and to clamp rigidly the end 
points of the units. The clamps were then removed partially., 
at the different points in turn permitting a change in the 
value of "one displacement coordinate at a time. -When. one or 
more of the units are. loaded with distributed loads, as in 
the present case with the constant shear flow, the units can 
be considered as curved beams the ends of which -are rigidly 
clamped to the constraints. The statically indeterminate re- 
actions ( i ncludingreact i on moment^) at the constraints can 
be calculated with -tK'e aid ■ of the classical beam theory. 

Forces and moment q "equal in magnitude and opposed in sense to 
these reactions ar'b the loads transmitted to. the constraints 
because of the distributed loads of the curved beams. These 
forces at the constraints do not differ in any respect from 
the concentrated external loads of the rings discussed in the 
preceding chapters. They can be taken into account, therefore, 
in the manner shown before by relaxation, the growing unit 
procedure, or t>y solution by matrices. The graphic method of 
reference. 7. takes care of distributed loads directly. 

The prp^'em at 5 hand- ‘reduce s , therefore, to. the calcula- 
tion of the ■^'taticall-y’^hdeterminate reactions „of curved 
beams. I n\ reference 5' "'the' calculat i ons were carried out for 
beams having arcs of a cifroelte- for their median line the load- 
ing of which consisted of a? suni f or.m shear flow. It is be- 
lieved that any ring can be f xyarpr es gated by arcs of a circle 
with an accuracy sufficient f or ’.engineer i ng applications. 
Similarly, any variable shear 'flow over an arc may be replaced 
by a number of constant shear .flows •, acting upon portions of 
the arc. 

t I . esiis* 

\. •» > *2 

Torsion of.' a Cirpplar Bing 

r : ■; ■■ - 

Basic data .- Figure 10' Shows the circular ring of a fuse- 
lage apd its loading. The latter consists of two 100-pound 
vertical loads of opposite sense, and the balancing shear flow 
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ransmltted to the ring, from the-., sheet co'veping. The interc- 
ity ofthe shear flow' 'is. - m 


t 

eity 


q = T/2A = (100 x. .20 x ,£)/( 2u.x 20 s ),-.•* 1 .1253954 1 pdunds ' per inch 


. t 


* ' (50) 


The cross-sectional area and.. th? moment of inertia' of the 
section shosifn An figure 10 are' ’ - ~ — — r— 


A s 0.250 square inch 


I = 0.1666 inch' 


The latter is tia-ken about the horizontal centroidal axis. 

The ratio f of the effective shear area to the -total cross- 
sectional area of the section %as calculated from equation (n) 
of reference* %:* * * ~ - -. j : *- .* • r 

'■••• • £ « 0.160 


Because of the antisymmetry only- one-half of the ring 
need he conqider^d-,.3 * The half'tihg is" Subdivided in two units, 
one extending from point A on the axis of symmetry to the 
external load (point B), the other ,£rqni i; the load to point C 
on the axis of symmetry of 'the ring. '‘ The angles subtended by 
the arcs are 45° and 135°, re spect^ve,^ ■ .The parameter 

. .. -j *• 4 -- *• . - . . 8 . -t 3 l ; - < 

.J. Atf/i' 1 

has the -iralue s " 370'. 258 and 333^.32, respectively, for arcs . -v. 
AB and. BO . ;< With - . * • * • - — . *.- 

^ _ . t A . ■ * » - . 7 i* 

45° = 0.7 8539 8175 Jr adiah and .135°. 2.3561945 radians- • 


and K ..are calculated’ as '’jLn the preceding chapter: 




K = 1.0120785 


, a • -i v 


V = ‘0. 9912535 

* -• 


These values are valid for both arcs, 


t i j ? - 


■ s 


I nf lu erice d o e f f i cl eat !> ■ The Influence coefficient s 
were calculate^... with the fe.id-’'of ’ ref er ence . 5. ifc- the'"- same man- 
ner as was done' for the. egg- shape ring of .the ; preceding 
ehapt er. 'Since the . cr oss section " of thq ^^iag". is constant , 

SI is a constant fact or -in every influence ; coe ffi'CfSttt’. It 
is convenient to assume its value as 10 s pound-inches 2 
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rather than to use the actual 
is taken as 1CT 3 inch and the 
the influence coefficients of 


value. If tji.e unit displacement 
unit rotation as 10~ 3 radian, 
arc AB become: 




-nh M = 281.950 

^ r M = 

-42.385 


49.079 

1 


rr M « 11.178 

any s» —29 .96§ 

nr ? = 

-20.559 

-4.733 

nty * 

' 52t896 
64.675 

>(51) 


O f* l 

• /•*..* tr 4 1 




'rbm * 6.633 

• 'it t. /I b l _ r 


-22.441 

■» i 

t tj =* 

51.516 

j 

. 


'■i'> t .*.■ *•.».. , » 



* ' " w * 




• The, i nf luence c 

oef/icients of arc 

BC are : 



■ -• - 

■ffhjj * 157.899 

nr M = 

-10.942 

£t M =* 

t *«v 

1.563 



rrj,[ =» 1.160 

" ’* «> • / ' * « • » 

n 

X 

r * 3 

-0.419 

tl M =• 

0.382 

^(52) 


: nH|i' 56.5117 

nrp « 

-8,. 8.42 

n"?j> '= 

6 .63 2 

r r - r-» 


•*' 1.116 - 

rt,j, * 

-0.524 


O'. 068 5 

» " 

• .. A ■ 

_ * 


Operations Table 


The operations table can be' established now according to 
the principles discussed earlier. Because of the antisymme- 
try points A and 0 cannot move vert i cally ....( See refer- 
ence 9.) Consequent ly , v^ and vq are not listed in the 

operations table ( tabl.q r ll ) , Moreover, there is no need for 
balancing out and ‘fig in the - f elaxati ons since ' what ever 


values are obtained .for them during the relaxations of_ the 
right-hand side of the ring, valtie's b£ the same magnitt.de and 
opposite sejse would be caused by the antisymmetric distor- 
tions of the left-hand side of the ring. It wa& also found 
convenient to move point C as well as A during the relax- 
ations. This is a deviation from the procedure followed in 
the relaxation of the egg- shape ring. 

* v f ’ ” if f ' ' i •• • fj i . . ■ 

Effect ■ of' the shear, flow. - The next step in the calcula- 
tions is the 1 '.d.e,t ermifnati on of the effect of the distributed 
shear. .The reactions*. of the shear flo.w at .the constraints , 
•act ing- Up'Oh' 'the curved' beams and taken, according to the beam 
convention 'o-f ; .signs shown in figure: 6a, can be calculated 
from equations (33) of reference 5. Tor arc AB the varlues 
of the determinants are: 
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v 1 


A q - 0‘:li'48'63e& i r F 'S ,r 2;&932?2? 

-A- = 0.00058&'7g§l'-'Sr 3 -- = 5.3093649 

q 

A r = 0.00893852? qr 3 r = '4 . 0228506 

r q 

.-A* = 0.05622000 qr 3 = 25.3078929 

“ q. • 






(53) 


Di vi si on ' of the numerator determinants hy, the denominajtoi’ de- 
terminant gives the moment, the radial force, and .the , tangen- 
tial force acting" upon the, "beam at the, support toward whjch' 
the shear flow is directed: ' ' . 


n v = -1.8350724' inch-pounds 


1 


= 1 .3904153 pounds 
t = -8.7471509 pounds 

H. 




-X54) 


The corresponding values for arc BO are: 


■3M 1 ’r.8 


<A q = 1.26215684 r 


26. 2431368 


-An = 0.1489809 qr 3 = 1341.3002 
q ^ ci«* 

A r = 0.7079546, qr 3 «.318.6$fS8 
*1 

-A t = 1.3379626 qr 3 = 602.294^1 

vc 

...» . r r 

r — 

n^ = -53.135241 inch-pounds 
r^ = 12.624880 pounds ^ 

t_ = -23 .859745 pounds 

q * 


'•(55) 




(56) 


.. As mentioned, before , these forces and moments represent 
the reactions upufi'-the' heams-% 'forces and the moment act- 
ing upon the constraint at section A have the values given 
in equations (54) multiplied hy -1 . If the signs; are now 
converted' to aglfree^ ’wltfc.E the rlgi’4 ’ frame convention of figure 
6h, the final values "become: 
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*A - -1.8350734 inch-pounds 

E a = -1.15904153 pounds " > 

T a = -8 .747 15 09 j pounds'* 


(57) 


In section B the reactions upon heam AB can be ob- 
tained by multiplying by -1 the values given by equations 
'(54) since the shear flow is directed away frojn point B ., 

The .signs of th& reactions so ; calqulat ed are i,n. accordance .. 
with the "beam co'rive'nt i on . The different st epq. --needed for ob- 
taining the effect of beams AB and BC upon .,the constraint 
at B are summarized below: 


; f E B 

‘ i ( ia.-lb) . .-;(lb) 

Reactions upon AB, beam. ■ 
convention, shear flow 
away from B (from equa- 
tions (54)) 1.8350724 .-1.3904153 


T 


B 
(lb) 


..8^,7471509 


Action of beam AB upon, 

constraint, beam conven- ~ 

tion -1.8350724 

rx 1 , 

Action of beam AB upon r 

constraint, rigid frame ~ 
convention -1. 8350724 


1.3904153 -8.7471509 

1.3904153 -8.7471509 


Reactions upon BO, beam 
convention, shear flow 
toward B (from e qu a- , 

tions (56)) 7 * -5(5.135241 ‘ l r fe . 624880 -33.859745 


Action of beam BO upon ' 
constraint, beam con-,.,, 
vent i on " ^ 


'5&V13524 1 


r.i " 

' -rI2. 624880 


23 .859746 


•’Action of .beam BO. upoji.-r . . . s 

- 'constraint,, rigid frame. ' 

' •convent! on -53,135241 -1.2.. 624.8 8 0 “ -23'. 8 597450 

” . .• T ■; - . - • ' - ■■ — . ' - ] ‘ ' • ; 

'T'ot.a-'l. Unbalance s .due to 

shear flow -54,970313 -11. 2344647. ? -32 ,'6C[6’8959 
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Ihe 100-pound vertical downward external load at. B 
also contributes to the tangential and radial unbalanced 
forces. If this is resolved in. the radial and tangential • 
directions, the result is 

= 70.7106.78 pounds’ ■' T.^ * 70.710678 pounds 


Hence the total unbalances are* 


JJ B = -54,970313 
R B ** 59.4762133 
T B = 38.1037821 


inch-pounds 
pounds \ 

p ounds 

' J ■' 


(58) 


In section C the reactions' are the values' given by equa- 
tions (56) multiplied by -1 since the shear flow is directed 
away from point 0. The action upon the constraint is ob- 
tained by one more multiplication by —1. The change from the 
beam convention to the rigid frame convention does not entail 
any changes in these values. Hence the final values of the 
unbalances are. .... 


K 0 * -r-63. 13 5241 inch-pounds 

r 

Eg i= 12.524880:'pnundrs 

- ' *. , i 

Iq .= -23.859745' pounds 

» * ! - .r t \*.z ' 

: ! i . r ■"* v ry • • - 

Solution by relaxations. - The i .relaxation was carried out 
in table 12i It may be noticed . that ^it was more advantageous 
to carry out operations on the' ibnge'r , than on the shorter 
arc. Moreover, sequences of the rotations wq, wb» w a were 
found ueefuli „ .... 

’ ' • * ? . , 

Rroia the total motions, of , A end B the radial force - * " 
at A can-,, be calculated. -jPh© motion's "are , as computed from 
table 12! . , . ‘ ’ “ '■ " -* ■ - *• : 

*' r •*> *f * »» -4 ? .. - r p . L 

w A = 1.15 X icf 3 radian 
,u A .. -4..8J x icf 3 inch 

w B = -2.96 x 1CT 3 radian 

V B 3 U B 3 0 
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• ~ ~ '-m i ■ ■ 


These notions md^t' be mult ipl i e d . by .the respective 'influence 
coefficients as lV given is equati ons . ( 51 )■ . Bhould , remem- 
bered, however, that in the calculation of the influence eo- 1 
efficieints the beam convention was used, bo that the forego- 
ing mot-i one mustcje. transformed to conform with the beam 
convention. Reference to figure 6 discloses that the sign-s 
of w A and u A must be changed, while the sign of w-g 

remains, unchanged. Consequently, 


w A = -1.15 x ; 10 
W-g » -2.96 X ' 1 0 


,-ar- 


-3 


radian 


radian 


+4.87 X 10' 


“3 


v B =-u B * 0 


inch 


With these values the radial force becomes 

R a = li'IS .'X 43.38,5 ;- -4 .87 X 20,559 + 2.96 x 4.733 » -37.3699 

\ 3 l _ po uagei; 

This l 
radi a 
( 54 ) ? 

R a = -37.3699 + 1.3904153 » -35.979 pound 

; j a : 

.1 

• SAnce the tangential force and the moment at A have 
the value zero, the bending moment distribution along the 
ring can be calculated with the aid of the elementary methods 
of statics. The moments are caused by R A and the shear 

flow along arc AB, and by R A , the external load, and the 

shear , fl o'v ) al o'if^ arc BO; The actual calculation is not 
shown. h,er'e . ' "T.tT --i'ill be taken up in connection with the ana- 
lytic.al s i 9 i ,lutd t bn ‘ of the problem. The bending moment distri- 
bution i B' siidi^n ’in' figure" '12 . ’ 


force, caiugad^by the 


1 f or 
n ■ 


ce cstfi Sed ■ by- ; the 
Ofrdther. l 


ft.lt ogdther, 


distortions, has to be added to tH6" ' 

shear flow and given in equations 

• * s. <f 


Solution by the growing unit procedure .- The problem was 
also, worked ^o\it' by the growing urti't method. ^ Inrthese calcu- 
lations point ' r C was considered as rigidly ..fixed. This fix- 
ation means only that the three possible rigid body displace- 
ments of the ring in its plane were, eliminated. The changes 
in the shape of the ring caused by*:the loading were in no way 
restricted by the assumption. 

A V " / «■ t* , r; 

As a first step, section A was rotated through an angle 
of 10~ 3 radian. The effect of t hi ; s. r-et«.t,}.'.on upon B is the 


v 
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appearance of the unbalances listed in the third to fifth 
columns of the operations table (talkie ll) .. ~ ’ * 

*T ‘ ■* / * _ **■■■- ** r P"?f 

Heit ; section. B was displaced and rotated, while. seo- ? 
tions A and 0 were held fixed, in such a manner as 1 to balance 
out the unbalanced forces and moment at B. -The necessary ' . _ 
mot'i onS-'"of .. B were- calculated from the following three si- 
mult anbo-us linear equations: • o • • 


•s ■ -439.849 w B + 31'. 443 v B - 50. .642 U B 
31,443 W B — 12.338 v B + SO.. 14> -Ugt 

"• . :*t * •* v •» | i 

-50.642 w B + 20.14 v B - 52.618 u B 
The . soluti pn of these .equations is 


= 29 . 966 

t ■ , 

.4 .,733 ’^■[(60) 

= -64. 6^5 



Wg = 16:1113569 2 x 10 3 . radian^ ,. t 
v B *= 4.035013 x lO" 3 inch l 
u B = 2.880754 x 10“ 3 inch 

. J 


(61) 


The effect of these - motions upon - A can be computed from 
the operations table: 


N a = 29.966 X 0.11135692 - 4.733 x 4.035013 + 64.675 X 2.880754 
= 170.552 inch-pounds ' ' 

T a = -64.675 X 0.11135692 - 22.441 X 4.035013 + 51.516 X 2.880754 
= 50.653185 pounds 


The effe-ct of.' thfe -unit-.- r.efeatijon • of A upon A (when B is 
fixed) is listed in the operation's table: '• ' 

^A = -281,95 inch-pounds T^ = -49.079 pounds 


Altogether., : 1 . 5 " ' • " - " • 

- ■■ s ' T r r f C: • 1 • o ■'* 

N a -281.95 + 170.5 
T A = 50.653185 - 49.079 '= 1.574185 pounds^ 


f i 


= -381.95 + 170.552 = -111.398 inch-pounds, 

7 a 'n. - v (62) 
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These a:.*e' the inomeAt ■ and • the tangential- for ce obtained at A 
when A is rotat ed. through a ' unit, angle ; • and B sores in 
such a manner as to cause vanishing unbalances at B. In 
other words, -N A and T A - of .equation- (62) are influence 

coefficients of the expanded unit arc AG. - 

Thei influence coefficients of ai*e‘'AC cor re spending to 
a tangential displacement of A can ;: 'b&' 1 btSkihS'd" ih a simi- 
lar manner. When section A is displaced tangentially through 
a di stance. of 10“.® inch,. the unb^l^nce s^at B ape the items 
listed in the ' thirdV*’ fourth , and fi'fth’ columns of’ the second 
row-of the ^p^rations table.. These unbalances vanish if the 
motibns of 'B'“ satisfy the f oil owi ng‘ equat i on s : 


-439.849 w B + 31.443 v B - 50.642 u B = 
31.443 w B - 12,338 v B + 20.14 u B = 
-50.642 .w B + 20.14 x B — .52,818 u B *= 


-64.675 

22.441 : y~?63)' ; ’ 
-51.516 


The sdlut ion is 


w B a 0.022242916 X 10“ 3 radian 


-. » ■- V B 


-0.53023 016; x- 10 


u B = 0.7546986 x 10 ’ inch 


10“ 3 ;i*nC'h vs “:y 
.-3 . ... ' 


J 


(.64) 

o 


- t n - «? coo -Cl . 

The effect of these motions upon A can be computed from 
the operations table: ?r ;• 


a -j L<i - 6J0-. - f. *' ?■ ^ X i f . .* ■ . 

N a - -29.966 X 0.022242916 + 4.733 X 0.53023016 + 64.675 X 0.7546966 

i Ji!-*. ’i 

= 50-653135 inch-pounds 

T a =064.675 x 0.022242916 +''22/44l"-x 0; 53023016 + 51.516 x O.75469S6 

; f ^ ; ; I ,l > .. .a ■ * f 

= 52.216509 pounds 


The effect of the unit tangential displacement of A upon 
itself is again taken directly from the operations table, 

N a ss -49.079 inch-pounds T. * -52.296 pounds 
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Altogether. 


• I i* - J i £ ^ 


i ■> 




■50 . 65£i8'5 ^9'. 079 ; = 

52.216,5^9 - 52.296 = 


1.574195 Inch-pounds 
-0.079491 pound 


(65) 


V • , . •* • 0 /■ /*'-■■■ . ^ 

These ia-re the influence xoeffici ent s' Vf arc AC sought. 

'• . *r . 

In the next step 1 , the external loads acting upon 3 
are balanced out by moving section B. The necessary motions 
can be calculated from three linear equations the left-hand 
members of whicii are idehti ca‘i“ trith' those of equations (60) 
and (62). The right-hand members, are the unbalances of equa- 
tions ( 68) : 


54*9 7033.3 . 


-59 .47621 
■=38 . 103 782 . 


The solution is 


w B «= 0.2815947 X 10 -3 radian 
v B <= 16.7320519 X 10 _3 .±‘nch 
u B = 6.8574788 X 10~ 3 inch, „ . 


( 66 ) 


The effect of these di splacement s'- upph '■ A ' is: 

= -29.966 x 0.2815947 - 4.733 x 16 .'7320519 + 64.675 x 6.S5747 SS 

= - 335 . 2764 , inch-pounds : ■ 

T a = 64.675 x 0.2815947 - 22.441 x 16.7320519 + 51-516 x 6.25747SS 

= -4.001964 pounds 

. , '■ £ ■ ; 

These are the unbalances caused at A by balancing out the 
unbalances at B, To these must be a&ded^the unbalances 
caused at A by the shear flow-as giveh'by equations (57). 
Altogether, the unbalances at A are: 
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1J A = 355.8764 - 1.8350724 = 354.04133 i nch-p ouri'ds 
T a » -4.001964 a ...7471609 <= -12. ,749 114 pounds 


(67) 


The unbala’rtc'fe s* quot e'd In equations ( v 67; must then be 
balanced out by moving A which is now considered as the end 
point of .the, esp^nd^d unit. AO. The— requirement of. vanishing 
unbalances at A can be written with the aid of the influ- 
ence eoe f f i cl ent s . jof .equations ($- 2 .) and (65) as 


■i ') :> » ii 


-1111398 w A + 1.574185 u A * -3 54 .T>4133 
1.574185 w A - 0.079491 u A = 12.749114 


■A • ■ 


( - 


The solution of these equations is 

w A * 1.2660489 x 10 3 radian 
u A = -135.31193 X 10~ 3 inch 


(69) 


The values contained in equations (69) completely de- 
fine the deflected shape, of the- ring. Equations ( 6 l) and 
(64) contain the data necessary for the calculation of the 
motions of B due tp the displacements of A;'’ 

w-g = -1.2660429 x 0.1113569k - i^. 31193 x 0.022242916 = -3.1507157 
V B = l. 2660429 X 4..Q.35013 + 135^31^93 X, 0.53023016 = 76.2560744 
ug = 1.2660429 x 2.220754 - 135.31193 x O.75469S6 = -92.4724222 

' t * ■ *** 

To these motions those undertaken when B was balanced . must 
be added. The latter are stated in equations ( 66)1 Alto- 
gether, x • r^v . ■" 

W-g = -3.1507157 + 0.2815947 = -2.869121 X 
v-g = 76,. 8560744 + .16 . 7C520519 * ,93.. 5881246 
u B = 7 9.8 ^4724888 '+ 6 .8574788 = ,,-91.615010 


-3 "'i 

10 radian 


x -IQ": 3 -., -inch J> (70) 


x ,.l;0" 3 .;inch 

- j" • : n 'J m - ' 4 


: . ;i i 
• « ; i 1 / ** 

* % V ■ . — 
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Equations (69) and (70) represent the complete specification 
of the distorted shape of the ring since section 0 was con- 
sidered ae rigidly fixed. The value of R& can. now "be cal- 


culated easily with the aid of the influence coefficients 
listed in equations (5l). It must he remembered, however, 
that the influence coefficients are always calculated in ac- 
cordance with the beam convention. In this convention the 
motions of. end section B of beam AB have the same signs as 
in the rigid frame convention. The motions at A,- however, 
must be transformed to reads 


- 3 


W£ = -1.2660489 x 10 radian 


Ui = 135.31193 x 10"" 3 inch 


(71) 


Hence, 


Ha = 42.385 X 1.2660489.- 20.559 X 135.31193 + 4.733 X 2.869121 
+ 6.633 x 93.5881246 + 22.441 X. 91.61501 = -37.93446 pounds 

It must be remembered that an. unbalanced radial force of 
-1.3904153 pound was applied at A because of the distributed 
shear, as stated in equations (54). If this Is added to the 
value calculated above for Ra* the following final shear 
force in section A is obtained: 

Rj^ = -3 6.544045 pounds (72) 

Knowledge of the shear force' in section A permits the calcula- 
tion of the bending moments in the ring. The calculations are 
not shown- here . They are .discussed, in connection with the 
analytic solution of the problem. 

The forces and the moment. i,n section 0 'are now deter- 
mined since they can’ serve as a check of the accuracy of the 
calculations. If the influence coefficients given in equa- 
tions '(52) are used, the motions of section B as stated in 
equations (70) must be transformed to conform with the beam 
convention: ^ 

w B = 2.869121 x 10~ 3 radian 

■ . • - Vg = 93 . 5881246 X' i0~ 3 inch 

u s * 91.615010 x 10" 5 inch 
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Then, 


N c = 66.6117 X 2.869121 - 8.842 X 93; 
a -67.77655 inch-pounds 


58 81246 ‘+ 6.632 X 91.615010 


i i ■. 


R c = -8.842 X 2.869121 + 1.116 X 93 . 5881246 •- '0 . 524 X. • 91 .615010 . 

* • ■ ' ■ : ’ / ■ ! ■ 

= 31.069316 pounds • » ' • a •* ‘ * 

T 0 =,6.632 X 2.869121 - 0.524 X 93. 5881246 + 0.0685 X '91.615010; 

. ■ 

= -23.736486 pounds 


To these values must he added the forces and the moment caused 
at 0 "by the shear flow in bar BO. Equations (56) contain 
these Quantities. The signs, however, must he changed since 
the shear flow is directed away from section 0. Addition gives 


E c 


-57.77655 + 53.’l36241 = 
31,069316 - 12.624880 = 
-23.736486 + 23.859745 


“'i 


-4,641309 i nch-pou-nds 
18.44444 pounds ' j>*'(73) 1 

0.123259 pound 




The terms Nq and Tq r ghpuld he zero because of the. anti-”' 

symmetry. Since the maximum bendlhg mon^e.nt in the .ring is 
-482 inch-pound, the deviation fro'm zero is only about 1 per- 
cent of the maximum. The deviation of Tq from zero is 


about 0.1 percent of the applied load” of ‘ 100 pounds. The 
check of jRq can he made as follows; - Rq is -54.98849 


pound.- This quantity represents the total .shear transmitted ■ 
thr.ou-gh the plane , of symmetry of the ring. ; The shear is act- 
ing upward ..be cau se of .-.the negative sign. . The vertical com- 
ponent of the shear flow act ing .upon the right-han'd side -.of- _ 
the ring is 1.1253954 x 40 = 46.015816 pounds. This is also 
directed' upward.; Altogether, the upward force' upon -.one-half 
the ring is • ... ... :**• - -• 1 r-f 


54.98849 + 45,0158.16 = 100. 004, 306 pounds".' 


r 


r* »«» 

i 7 




It should he equal to t.he 100-p,ound. downward external force 

1 ’ : T - \ . 1 ”, • \ * V . * , V. • 
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Matrix solution The problem was also solved "by the 
method of matrix transformations as shown in the appendix. 
The motions obtained are listed below:. 


w^ = 1.2646 X 10 3 radian i 
u A = -136.3985 x 10“ 3 inch 
Wg = -2.8709 x' 10~ 3 radian 
v B = 93.6280 x lO” 3 inch 
u B = -91-.6842 x 10“ 3 inch 

i 

Section C was again held fixed. 


(74) 


The displacements according bo equations (74) are in good 
agreement with those found by the growing unit procedure and 
listed in equations (69) and (70). 

Analytic solution .- In the analytic solution again the 
lowermost point of the r.ingj i« ; consi-.dered fixed. (See fig. 
11.) In the section at the topmost point the unknown tensile 
force T, shear force T, and moment .N must be applied. 
Because of the anti symmetry ; 1 however , H s T = 0 . Hence the 
bending moment is caubedn only by V, the external concen- ?, 
trated load of 100 pounds, and the shear flow. The contribu- 
tion of the latter can be calculated with the aid of figure 
lib . . . ... . 

j . . r / 

The shear force acting upon the inf inite si'mal element of 
arc subtended by the angle d| is qrd£. The moment dM 
caused by this Bhear f^cjs at. .. i_s " / * - 1 1 r 

dM = aqrd I * r[l - cos ( cp - £ )] qrd £ 


The moment caused at-,r cp by ’’fell the shear 

|=0 to £ = cp is consequently 


<P 

Mqj = qr 2 [l - cos (cp - £) ] d = qr s (c£ - 


The total moment is 


flow from 
sin cp) (75) 



6,4 


NACA T3J No. 998 


... J |M l = -4r a '(0 - 'sin‘^) + Vr • flin cp/ 

. _ •. n t 1 ‘ " 

.... ,■ I I 7 ir ; ‘ 

when 0 < cp < 46° 

. c ; < : . I x - 

M a ss cl r 2 (op sin cp) + 7r sin cp- + 100 r(sin cp - 

"T r,t *- •*» ' 

when • , • • -..^S 0 < cp < 180°' 


... -M 
w : v a * 


0.707) (76) 


.vr 


The strain energy of bending is 
1 


4 6 ° • 


.180 


U = 


2SI . 


M, 3 rdcp + — 

1 253 1 • 


Mg rd cp 


4 5 


According tp Castigliano !-s .principle ; , ' , ' 1 ’ ' ' 

\ 

.... dU / dV * 0 ’ 

; : • : \ ;• icr -• '1 1 • ’’ ‘‘ - 

The differentiation may. be carried Qiu.-.'fcr b’ejforef the integration: 


dU. ' ._r_ 

4V EX 


4 5 ' 


.. , " •Xf<jr ! v . s * ' • ’ 

, '1 8 0 . . . 1 .. , „» >•.. 

-■dMi _ -'t T- af4 


, «* + - san a ? *' 0 


If the operations i n^i c.at pd , e , carried out, the following 
equ&t i en~ I s obtained: ’ ‘ ’ ,r J 


Since 


qr + V + (200/tt)[(3tt/8) - (l/4) - 0.707] * 0 


q = ; (1'00 X 2r x ,0.707)/(2ur 3 )’ ' 


the shear force becomes 


V = -36.5766 


(77) 
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^ ’ : ‘ .?-*? ; i. * ' w 

It may "be seep, ; t)iat this value agre'eB w'el'l wi’tfy ?those .db~ ' ** 
tained for “ i n t.he earlier sbiubi onls . Substitution of' ; 

q-, r* ' and t X . -in- equat i one (76) jLeads*to the-, final ■ expre s- 
siohs f op" t^e 4 - binding moments: •' * “ " 


M, = 450.16 (cp - sin cp) - 731.^53 sin cp*= ’450.16 cp 

: ■ .*n ' " ■ ' 


when 


0 < cp < 45 


■-^l-r81, : 69 ' 4 : in cp 
■ f i ■■■ T ■ 


o ■ 


(78 ) 

: r ‘ 


M 3 = 450.16 (cp - gin cp) - 731.53 sin cp + 2000 (pin cp n *; 0.707) 
= -1414 + -450.16 -cp + 818.31 sin cp 

* i 

when .. ... -v ' ‘ cp < 180° ' 

... , r 5.\“ ■ — 

.e • ■ y . . . 

Th'h^ bending' moment diagram is shown in figure 12. 


( 79 ). ^ 
ztr.* 

S * ■' • V* 


r $ 


' J * ' 
lor ! 


r sietn . of '■t’h.S Braced Circular ' Ring 

- ' ■" . ;-.it . • 

In actual airplane fuselages the rings are often inter 
nally braced when large concentrated loads are introduced 
into them. Such a braced circular ring-"- Is shown in figure 
13. Since the addition of the brace causes six more stati- 
cally Indeterminate quantities to appear in the calculation 
the work involved in the analytic solution of the problem o 
the bending moment distribution becomes very, muqlj more in- 
volved than It was in the case discussed l ; h the preceding 
article'.' Because of the .antisymmetry the actual work in- 
volves only the determination of one unknown quantity in th 
ring proper and three unknowns caused by the braces. Ihe~ 
increase in the.work/is theh'aV follows 


a, 

f 


.. -r- . T 


‘Calculation of 5 moment diagrams instead of 2 


Calculation of 14 definite integral's* instead of"/_2J^ t l 

n • • ~ */.U v. c* v-fp-e 

SojJ.ut.ion of 4 simultaneous equations- for 4 unknpwnjs., . , *■<?.: 
„ „ .instead of getting 1 unknown from 1 equation ° ‘ . tAj'J 


In t.he methods 'proposed in this report, the additional 
work consists only of the calculation of the influence coef- 
ficients of the braces. Because of the effect of the braces 
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the numerical values’ id 'the operations table will.be diffdr-c:' 
ent , but the number of -entries in the operations table .willj >t 
remain the same. The ■ s'olut i ; on of the problem by r elaxat i oq>.9io 1 .-4 
or by the growing unit procedure will, therefore, not in- 
volve any more work than in the case of the ring without the 
brace. ~ t 

The length of the diagonal 1 b L = 36,955182 inches. 

Its cross-sectional area is A » 0.125 square inch, and its 
mail mum moment of inertia I * 0.020833 inch 4 . Since in the. 
calculation of the Influence coefficients of the ring in 
equations (51), BI was assumed to be 10 s pound-inches 3 
while I was 0.1666 inoh 4 , the calculation^ will be con- 
sistent only, ■..If; 1 v'BI ■ of the bra’ee' is assumed 


\EI a (0.020833/0.1666)10 * 125,000 pound-.inehe s 8 

r ** . 

This means, of course, that the deflections to be calculated 
would be the actual ones if the material had a modulus of 
E = 6 X 10 s psi Sfiith aluminum ' alloys, or steel, the values 
must be multiplied by a constant factor. 

The ratio C A* /A ■ i e again 0.16. Consequently, 

** * 8194.1131 and the product *Y | s 1311.0581., The influ- 
ence coefficients fo£ the .'movable end of the straight bar 
were calcu.lat e.d, »f jjoir - je-quAti ons (25) of reference ' 5:. 

..... . • • ' : ; .: • 


,v> 




nn 

rr 


M * T2 : .;43’79 ; *Q.. 5*419.4 r ; - .nt M - =• 'O' > • ; 

M = SOya 9 ! 49’ 

“ l " . r» r 1 V I' & ii * S* J * r r h 

> n * *t 4 r i ;* ■ - •* • • •' ' * ' ' , , , 

... *♦**>*'*’ r ^ ^ ^ 5 *•/ • ^ ■' "■ * 

The inf luehfee'- '&’■?>'« f fl ci^’At'^! f ( 9r.' t&e\f ixejflU'.endr.: wete; determined ' 
with thb ;! add of ' e'^h^i Ons' Xas^r/of.-re.ferenEei'St-i ■;' * ,::ir " • ... ‘■ At 

*-'• 1 j{ '.i i ■; : ■• 

: a\- o ^ 1 • 


nnp = -6 . 67292 
rr™ « -0.029451 6 s 


nr F = 0.544194 


rt 


F 


fj a ■’ 


: . 1 


• i*- 1 *" * ^ ' ' l' ( 81 ) 

0 $ 1 - = 20.294 9L . . 


a r -■ 


Since in the { ,,pp^erp.J>lop.s jt.^bsl^. p«.int&' ' B' ■ k^d C ape;,',” 
moved tange-htialljr radially.: -.t.o.^t.ive fei'k=cie--, &: lt' i s’ neces- 

sary to convert*' the above influence coefficients intq val^s 
that correspond. ,$.o ,thq se-^i splacejfient is . Shs^apdfe', a' 

unit disijlaftPmfe^^^^B* ..tangent ieilly-.torthk cl'fcle 7 dQVnyar| 

a - .r.- o -Jrtt? ' - - • 
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to the. right is equivalent t o- a- tangential' (longitudinal) dis- 
placement of --sin 22.5° for the straight bar combined with 
a radial (perpendicular) displacement of- coi 22.5°' for the 
straight bar, when the beam convention of figure 14b is used. 
With 

: sin 22.6° = 0.38268.- '■ cos 22.5° = 0.92388 

r • - 

the following forces and moment are obtained at B, acting 
upon the straight bar and designated according to the beam 
convent ion: 

N sb = -0.544194 x 0.92388 =’ -0^5027? inch-pound 

R sb = 0.0294516 X 0.92388 = 0.02721 pound 

■r SB = -20-2949 x 0.382&8 = -7.76845 pounds - -- 

where the subscript SB designates that the forces and the._ 
m6ment r are basedupon the syst'em of ' coordinates of the 
straight bar;.. They must be ’convert ed to Conform with the , " , ' L 
system of coordinates based on the circle before they can be 
entered in the operations table. If the subscript C3 re- 
fers to the system o.f< ; ,coor;di-i?at e 8 of the circle, and use is 
made of the diagrams of figures 14c and 14d, the following 
values are obtained: 

■ ««. • • < , * i >’ •" / L -■ ‘ 

Nq B = Ug B = -0.50277 inch-pound . 

E ob — sin 23.5 — Tg B 22.5 . 

= -(0.02721) X 0.38268 + (7.76845) X 0.92388 = 7.1667 pound 

J *^CB * a SB 008 22,5 . ,®,8B ® in . .. .. . . 

S (0.02721) X 0. ; 92388 + ('? . 76845 T f x’ Cf. 3^26 8 * 2.9980 pounds 

» ’ . - ■ : .. ■ - \ 

Th - e forces ah'd r the moment just calculated- represent the action 
upon the ’structure. Multiplication by -1 cphverts them into 
the effect' of the’ structure up tfn' tlier ’const rain’t s , tahe’n in ac- 
cordance w'ith ther rigid frame 'convention based on' the circle. 
These values must* be added to those listed in't'&e operations 
table of the circular frame wl t'h out ■ bracing'' (table 11). The 
sums have been entered in the new operations table (table 13). 
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The •" effect ' of k ■ radial displacement and of a rotation..]*. ^ 
can he 'found in a similar manner. As a furthelr. example for 
the resolutions, the effect upon 0 (the fixed end) of a .*../ 
unit ; r.ptati on of B (the movable end) is nov considered.' 


A unit clockwise rotation of section B of .the. circle is 
the positive unit . rotat i oh of section B of the Straight’ bar 
as shown in figure 14b. The forces and the moment caused at 
C by this rotation are given in equation (8l): . ;• 

* - • 

Ng B = -6.67292 inch-pounds : : 


3*83 .» 0,544194 pound 




These quantities ar^" the forces and moments acting upon the 
straight bar at C and taken according to the beam conven- 
tion shown in figure 15b. A positive moment Ng B ah. C is 
counterclockwise] Hence, the moment . taken in accord- 

ance with the-rlgid frame convention of the circle has the 
opposite signs.;-,,'} _ ‘ 

= -Hg-g = 6 . 6729 2 inCh*-pQund s 


¥ 



The radial force R SB is resolved with the aid of figure 15c: 

RqB = ”®SB 22.5° = -0.544194 X 0.38268 » -0.20825 pound 

* ‘ J *’ ~ r • ’ ‘ - C' * 

T ob a -R sb cos 22.5° = -0.544194 X 0.925818 = -0.50277 pound 

-■: : r .*■ ••• •' ■ • ’• ' * * 

The values of 3 q b , Rq b , •* T q B - must be multiplied by -1 

in order to obtain the effect of the structure upon the con- 
st ralnt s . a^s,. •reqp.ij'ed; ;f or. ..the operations. table. 

The operations table for the braced circular ring is 
pr e sen£ ejl,. a^ . t abl 9 . 13., 1 1 _, V£ p -p bed • relaxations con- - 

tainedj ipl tapie 14 As, : may be seen-,’. the !' rb laxat 1 on b were car- 
ried jjut ].'in fi t he same manner,. as in Jhe.cabe of. the ring- with-' 
out brfici'. jig.,' The absolute magnitude s-. of 1 the; de fleet i ons 'were' 
found ' t o.'be , much . pmailer than befpye . r. ’For ’ the ghedi* ’ force at 
A the v.s.lue 6 .4964 p ound% -rwa s ■••obt aine d •; ~ * 4 : ’ 
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The distortions were determined alpj by the growing unit 
procedure. The individual steps in the calculation are not 
given here since, except,, fpr the numerical values, they were 
identical with those p're'sehted in ’the pr e-cedihgi-'apt i cle . 
Section 0 was again assumed to he fixed. The distortions 
were found to he: v .~ 


w A = OY'5748'68 59 x*. e 10;~ * radian 
U A = 6.499450 X 10~ 3 inch 


: ’-' n J = -0 . 04252815 x _ 10. 3 radj.ap > (82) 


.... = 0.21585307 X 10 -3 inch 

— 3 


. rUj T 

■ ^ fv r* is i l 


7.4627458 X 10 

,4 » 


Inch. 


' ' T ' 



Comparison of the values given in equations > wi th 

those contained in equations (69) and (?0) reveals the great 
stiffening effect ■of'fche .trace*. . 

w 1 ‘ ■ r. . 


The value of ’the shear, -.fop-pie t ^Ln section A was found to 
he -6.460 pounds. This value deviates from ttfatf obtained hy 
relaxations hy 0 . 6 .per cent.,' ' .T-he- "bending ..moment diagram is 
shown in figure 16. It may he seen that ^the D in-troducti on of 
the -brace s , re suit ed in the reduction of the maximum hending 
moment from -482 inch-pounds- , 


nds- : ..feo a . n 58^1. ( 4nch-pounds^ 

, %. , ~ “• *■ i 1 1 r. “ 

. ...•JTOCCC heda ; .( (in 1 a-. * v , r 

The problem was- al s o-. s olv^d by the . meth pd. of ‘ 

. « . » • * - . . : . ~ a Art n _ T 4 ‘ V * ^ K * 3 A 


matrices 


The resulting- distortions-', differed .Icantly- f r om 

those i’i st ed L * in' equ'a-t i ohe (:63 i • c in 1 | 

“ ; . .. ’* * . - Z irJSXlf 3 j_ . . ' 

, c . *» *( 3 f * O «. It ^ . . . , . 

■ M ,.T -r, i.~ i- * 


T or si.o-h: <o:f- the Sggt-SJiap.e Ring 


•t T 



i 

Th 

Her ;ih thi &" report when "bh^, hendfing ,njQ,m.ent.s. " caused "'by "tv/d 
equa^l f and opposite- forces were coated. jAd" . in t.he "case 
of '"the"’ 'twist ed circular rin|,- it f|,a : 

flow transmitted from the sheet, cpv 
the ring is constant and can he calculated* 


assumed' "ths^t "the shear 
esr^n^'rpf” the, "ftiselage to' " 


from the "f ori&ula 


- a •) 


- frjc* 


q. r 


iiS'..:.orr. iv <rj >=-/.* -.J 
T/2A°' : <i/: * . *U STif-r-J"* .. « 

■ t r . ‘.axle*-:*. * ?t;~s.cx 

• fcO'J V Jv r R. ■ • 3 C - : ».s 
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In- the present case the tor^e ; . 

, T = 100. x : -30. «. ; 3.000 inch-poundB : . . 

* - -r . * '' - - ' • r t 

while the area included’ by the ring - * • 

■ .#•;*« ’ 

A= 2757.642 square inches 

Hence , , 

3 000/(2 x 27.6.7 . 642) = f 0,. 54394298 pound per inch 

When the operations table-- was set' up, the shearing and ex- 
tensional strain energies were neglected. This is equivalent 
to setting 1/7 = 0, a'e ! wafi J Explained when the loading con- 
sisted of two equal and opposite forces. The angles subtended 
by the uni.t.s. .considered are: 

’ •* • ** ‘-.it. . , . . 

' 1 l ; » .. . * 

Arc AB 0.848.06211 radian 

... Arc. . ■ 0.52137719, radian^ 

Arc.. . D$ * • ' 1. 7721543 radians- 

• 1 - ■•«■■■»' s 

* • * • • . " * * ‘ •; * f 

The Influence coefficients were calculated 'And the oper-,' 
ationB table was established according l, f*o "th'e -prccedur'e d'i s- 
cussqd-. I n the earlier parts of. this report. The operations 
tabi,e...is reproduced a% table 15 . , served as the basis for 
the calculation of the di st or't'i’o : n|f .by t’lie growing unit method. 
The distributed shear forces were' taken r int o' ^account as 'ex- 
plained in connection with the torsion of the circular ring. 
Point .10 was considered as rigidly fixed. ? Because of the 
ant i symmet ry, sect ionJJi. could not move radidlrXy. 

f u c - J, *\ * * ‘ • •* i . , 

Solution dv - the growing unit method .- The 'application of 
the growing unit,' .procedure to thi6 particular prob J.em ' con- 
sisted of three .qiai'p jpart s . In the first, the expanded units 
were established;. iri*'the second, the unbalances were trans- 
mitted from sect i on ^ ' B , 0 , and D to sections A and'E; and In 
the last part, thq unbalance s at A were eliminated , 

The bending moment diagram was calculated and is pre- 
sented In figure 18. The dotted line represents the bending 
moments obtained in reference 7 .by the graphic method. The 
agreement is reasonably good. 
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Torsion- of the -Braced^ Fgg-Shape Ring 

' ' . f i ■ » ■* . t : ‘ ‘ ' *1 ' T * ■ 

Figure 19_ shows the braced” egg-shap-e' rin^ ; and-, its anti- 
symmetric loading ■ d! he. -cross -sep/t i.on of the brace was taken 
the same as in the - case'-of. the" .braced circular riiig.’ - The 
length of the diagonal is 64.98365,5, inches. - The- influence, 
coefficients were calculated and the -operations "table' was es- 
tablished according to- the principles discussed earlier. 

; ■ * * * . 

The data Pont ained in the Qpe.rabi'oijs : table (table 16)- 
served as -a ’basi s 'f or the calculat i qp of the deflections by 
the grofting-’uait -procedure . The results 


. ■.* ^ 


w A = -0.7703992 x 10" 3 radian . w Q , 0.,.1?370249 x 10 3 radian 

u A = -12.059232 X 10~ 3 inch-'- ' v» ' = '9 .243786 x 4 l<T 3 inch 

: *f ' '■ - r. ■ - - _ - . _ 

w B = -0.00896880 x 10“ 3 . radian u 0Di “ \ -K’3‘0'1079 x ID -3 inch-. 
v B = 1.9469763 x 10" 3 inch w D * -O‘vl'9693 24 X 10“ 3 -, radian 

u B = -12.20075 X 10“ 3 inch v D = 10.421443 X 10 -3 inch. ■ . 

The motions of sections A and B' in- commotion with th,e. 

influence coefficients permit t J h'e carcul'ahion .of the shearing 
force R* in section A: • 

. ’ ••• ■ t . .. 

4 * * »> ’• ’* 

R a = .'6 .3Q15 .p.ounds . 

.t .-*•« hi,.- 

The motions of B and the influence coefficients of. the 
straight -bar -BE suffice fqp ; \the calculation of the forces 
and moments acting on bar jSS' When’ these are known, the 
bending momeat diagram of the^ rljig can be determined. 

«i - ■ - , . • - 

“'The 'bending moment diagram. is p feb eh t'i'd*' in" figure . -.20 . 

It may bV> seen: bhajfc the br acing. was 1 ' Ver^r’ 'bA'h'e'l’iclal- in reduc- 
ing the 'bending jmoment s . •* 


GONGLJJSI OHS 

‘ ‘ - r -t . . 

■ ; •• ... * 

The subject of this red'irt is the numerical calculation, 
of the distortions of and thb’bending moments in rigid fr.a.mes 
(fuselage rings) of arbitrary "shape upon which known epncen-T" 
trated and distributed loads- are acting. The calculations 
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consist of twp major parts: the settings 'up of the equations 

characterizing the problem, and their solution. 

The equations are set up by dividing the frame i p t c . a 
number of .-beams eaqh having a c on st ant 'radius of curvature 
and a constant (often zero)* distributed 'load. The iJif-lue.nee 
c oe.fi'i c i-erntt s of the individual beams ca r n' "be easily comput.ed 
. ,,wit&...the' add of' the formulas, tables, /a"hd. graphs presented in 
r ef eren'c'e- .'5 . TTa'ey' 'can be combined into quantities that are 
entered in the operations tab 1-0': accord in'g to the suggestions 
made in'bhis report. The operations table, together with the 
external -ioa'aB, completely defines the problem, and can al- 
ways be established without difficulty. 

On the other hand, difficulties are enc ount er ed « in the 
s olut.1 on <':df -thWet of linear squat i one ..defined by the opera- 
tions table. The, fundamental reason for the' di ff i cult i e s is 
the fact thh’t in" most cases the bending, moment , the normal 
force, and the shear force in sections of the' frame are ob- 
tained as- small differences of large quantities, 

Three methods of ' solving the linear equations are ;p.r«'- 
sented. The first. s blje r eolut i on by 'Southwell 1 s method of 
systematic relaxations' which is a procedure of step-by-step 
approximations. Four numerical exampl.e.B give an Indication 
how the f ndivid.^ailv. .pt'ap.a 1 in the relaxation procedure may be 
arranged, so as- to u achieve a rapid convergence of the proced- 
ure. Iifmost cases, however, it is not easy to find the 
proper succession of steps. Because of this, the procedure 
of the growing unit was devised. In, this procedure, the in- 
dividual beams are combined into units of increasing length 
until finally the entire frame becomes a single unit. In 
each step of th.e .cal culat i on not more than three simultaneous 
linear equations musi.be solved. Finally, the third method 
is the' direct solution of the set of linear equations by 
means of*the matrix calculus. The latter two procedures are 
straightforward and can theoretically always be carried out. 
The work involved, howsver, ■ increases with the number of sub- 
division's chosen, and because ,of the small differences of 
large. ; nu-mber b occurring in the computations the use of cal- 
culating machines may become indispensable. 

It :L s of interest to compare t'he numerical procedures 
here suggested with known methods of analysis of the bending 
moments :Ln fuselage rings. Fully analytic methods.. have been 
published only for circular and elliptic rings - for re-i '• 
stricted number of conditions of loading. .Rings of ..ar.bit rary 
shape were discussed by, Lundquist and Burke .(reference 8), 
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Hoff (reference 7), and du Plant'ie.r .(reference 8), rings with 
additional internal tracing elements by Lundquist and Burke 
(reference 6), Simple ring problems can be solved by any of 
thege methods 'Vf^th comparative ease, sb that the procedure's • 
sugge st ed ' in thiB report show no particular advantages ’in 
these cases. The numerical procedures of this report become 
advantageous whefi ; 'one or mo^e int erhai’’ bracing -element s are ; ’ 
incorporated info' the 'frame . The additional elements ■'entail 
a slight increase in the w'ork ; of setting up the operations 
table, but the solution of the equations does not necessitate 
any additional ‘Work. If the eau&t f 3ng are solved by reia$a- • 
tifizil X the number of steps necessafy 'may ‘ even become .duller . 
Since’ : iii v th'e ” earlier methods ' t?he\‘work‘ involved . in the' calcu- 
lations increases t^a'pidly when redundsLn't bracing element's are 
added t o . the _fr^me the advantage of using the pr ocedur ^qf ^ 
this re'port ' incr £a£as '‘with each adddd bracing element ‘ 

- - . ' ■ " V - - ■ - >! 

■ .. 

Polytechnic Institute of Brooklyn, 

Proo-klyn, Hew York, February 1945. 
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APPEND. IX • 

DIRECT SOLUTION OF A SYSTEM OF LINEAR EQUATIONS 


As was shown in the body of thi 3 report, the operations 
table represents a set of simultaneous linear ' equations, ' 

The equations can be solved directly by the method 'O'f elim- 
inations. Tne only difficulty lies in the great number of 
simultaneous eauations and the corresponding great number of 
unknowns contained in the operations table. Because of these 
the numerical work'is far too cumbers one *u» : les s the dpdra- ’■* 1 ■ 
tions are carried out in a systematic manner. The systematic 
solution of sets of linear equations is one of the subjects 
of the matrix calculus. The elements of this discipline 
needed here are very simple. An attempt is made, therefore, 
to present them in. an easily understandable manner. 

First a set of linear equations will be solved by the 
methods of high school algebra. Then It will be shown that 
the same calculations can be carried out in a more concise 
manner 'through the usd of the matrix calculus. 

As an example , the following three equ-at i ons containing 
three unknowns .will, be' .considered,* 

• 2 x + y + z <s 8 - • (la) 

J » ‘ - -**’**?' * - " 

V • . ..! X + 2 y +.. X ;=*;9 (lb) 

x +: y .+ S', z . - 7 (lc) 

i 

Divide equation (la) by 2:: 1 

x + 0,5 y + 0.5 Z a 4 (2) 

I 

Subtraot (2) from (lb): 

... . 1. S’ y + 0. 5 z ;* 5: • • (3a) 

Subtract' (2) from (lc): 

, t • 9 • 5 y + 1 • 5 z » 3 ■ ■ ' - T - « : (3b ) 

.: r .... - . - • * ■ ' . r is- * * - .d-.r s* 

I.t may;.(be seen, that equations ( 3 ) d<? -raojj ^Qntain- . X ' 
any more, j Next y ■ can be eliminated fro.a equati on ; (3b ) if 
first equation (3a) is divided by 3 : 
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0*5 y + 0.166 z =* 1,66 (4) 

and subsequently equation (4) is- subtracted from equation 
(3b) J 


1.333 .z .= 1.333 . 


(5) 


Equations ('§)» (3a)', and (la). c t an now ,he ^easily solved for 
the unknowns. Obviously, equation .(5) yields 


z s*. 1 . 


( 6 ) 


Substitution of equation (6) in .equation' (3a) and solution 
for y. gives 

. j' ■. ' ■ : y a 3 ( 7 ) 

Finally equation' (la)' yields after substitution of equations 
(6) and (7) : ... 

■i- ■ ■ J 1 • • • '■ x ’ '• ' . ' 

* 'v=. ‘ ■ x * 2 ' (8) 

w . i ^ ‘ •, ' • * , » • • — . • * " * 

The same pr obleto • V’iil now be solved using the matrix no- 
tation. The first achievement of. the matrix calculus is 
econQmy-of writing;’ As long as the first column 'always con- 
tains the x terms, the second the y terms, and the third the 
z terms, it is necessary to rewrite j**,; y., and z. The si- 
mu 1 tqrqeaus equations are characterized by the set of coeffi- 
cients arranged in an ordered arr^y, ^ Tfau,s the left-hand 
membe^p, of equation (l) can' he represented by the set 


(9) 


It is customary to enclose the set of coefficients with- 
in brackets as shown in eauat-ioji A, ‘ o’r the set of coef- 

ficient's'. ’in brackets, is known as a "matrix." The nine num- 
bers in the brackets are the "elements" of the matrix. 


The second important device iij matrix calculus is the 
operation called "matrix multiplication." This operation is 
in many respects very different front, what is known as multi- 
plication in arithmetic and. algebra. It would be better per- 
haps to coin a new word for this operation in order to avoid 
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misund'er et andin'g. This would necessitate, however, ‘the memo- 
rizing of some artificial expression. As things stand now, 
the operation is known as "mul tipi icati on“ in mathematics, 
and engineers have to conform with the usage. They should 
remember, however, that matrix multiplication means the 
complex operation shown below, and not what someone might 
infer because of the familiar expression "multiplication." 


The product of two matrices is again a matrix, the coef- 
ficients of .which have to be calculated according to. the fol- 
lowing rules . . ...... 


a i 

b. 


where 


8-3 


d x d a 


e. 


e 


a 


f. 


d. 


S i 63 
h, h 


k x 
3 1 


a 


€3 

*3 


(10) 



s 

a i d 1 

+ 

a a ®x 

+ 

a 3 f X 

e 2 

s 

a i d a 

+ 

a s ®a 

+ 

a,f _ 
3 a 


s 

a X d 3 

+ 

a a e 3 

+ 

a 3 f 3 


s 


+ 

^a e x 

+ 

* 3 fx 

h a 

• 

=7 

*x d a 

+ 

^a e a 

+ 


i. 

s 

c x d s 

+ 

1 - 

c e 

3 3 

+ 

°3 f 3 


(11) 




The rule shown in equations (10) and (ll) can easily be 
extended to square matrices of any order. The "order" of a 
square matrix is equal to the number of rows (or columns) in 
the matrix. Two examples of matrix multiplication are given 
below : ‘ 


1 3 

x 

5 V 6 

• ' ■ 

(1X5)+(3X ! 1), (1X6)+(3X^) 

• 

8 12 

2 4 


1 2 


(2X5)+(4X1), (2X6)+(4X2) 

CS 

14 20 
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1 

0 

0 

0 

1 

1 

1 

1 

— 

1 


\ 

1 

1 

— 

1 

1 . 

1 

0 

- 4 

;0' 


•2 ' 

1 

0. 

0 


3 

2 

1 

1 

0 

0 

1 

1 

X 

0 

0 

X 

2 

p 

a 

7 . 

0 

0 

2 

2 

0 

0 

0 

1 . 


0 

0 

0- 

2 


i 

0 

0 

0 

2 . 


It should, be noted that, in general, the "matrix prod- 
uct" changes' if the factors in the product (that is, the two 
•matrices) are interchanged. Tor instance • ' 


r 

1 , 2 


r 

3 0 




X 


s 


1 1 


2 1 


5 1 



M — 



' - 


_ 


r -) 

3 0 


1 2 


3 6 


X 


a 


2 1 
__ _1 


1 \ 


_ 3 5_ 


i 


In equation (l2a) the matrix 

3 0 

2 1 


. (12a) 


(12b) 


is said to he "premultiplied," in equation. (12b) "postmulti- 
plled" by the matrix 

"l 2~ 

[l i_ ' 

With the aid of the .mult iplicati on rule of matrices,; 
equations (1) njdy now b4;W r i'fct8n in the following form: 


O. y, z] x 


2 




[8, 9, 7] 


( 13 ) 
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The “rojjr matr ix^ , [x ,* V;* z^" * i*' identical- vitkv.the square' 
matrix ' ' ' • •* x -i' .:■ .. ■-- a - \ . 




Application of the multiplication rule to the left-hand, mem- 
ber of equation (13) yields 


* *• *1 vJ ■*? 


2 11 


0 0 O x l 2 1 = 


(2x + y + z) (x + 2y + z) (x + y + 2z) 




The right-hand member of equation (14) can be written as a 
row matrix 

[<2x + y + z), (x + 2y + z), (x + y + 2z)] 
so that equation (13) is equivalent to 

[( 2x + y + z), (x + 2y + z) , (x + y + 2z)] = [8, 9, 7] (15) 

liquation (15) expresses the equality of two matrices. Matri- 
ces are considered e quad if, and only if, all their corre- 
sponding elements are equal .. I,f one equates corresponding 
elements of the matrices .on, the" t.wo* sides of the equality 
sign of equation (15), the following algebraic equations are 
obtained: 

■ 1 ~ fx • ' * 

2x + y + z = 8 

x + 2y + z = 9 (l) 

X + y + 2z = 7 


which are indeed identical with equations (l). 

It should be noted that when the matrix equation (13) is 
transformed into three simultaneous algebraic equations, the 
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matrix is read "downward," not from the left to the right, 
in agreement with the rules of matrix multiplication. For 
Instance, the matrix equation 


[x, y, z] x 


12 3 
5 4 2 
2 3 2 


Cl, 0, 3] 


is equiral ent % o u ! t^e" ’thVe e algebraic equations 


„x , + 5 y 


2x + 4y 

l 

3x + 2y 

and not to 

• * nat> x + 2y 

5x + 4y 

* X v + 

. 1 v 

2x + 3y 

Thi 3 . di.st in’pt.i on. cannot be noticed in the case of equa- 
tion (13), since there the square matrix is "symmetric," ' ■ 1 

which means that elements in it are equal if they are eym-.. \ 
metrically, situated, to the principal diagonal. The "priaei-'i^j 
pal diagonal" of a square matrix "is the diagonal passing ... fr.om 
top left to bottom right. 

The ^following results of matrix mul t i nl i ca t i on are worth 
noting! 


+ 2 z * 1 
+ 3 z = 0 
+ 2z - 3 


+ 3 z = 1 
+ 2 z = 0 
+ 2 z - 3 


+ V ■>' !• '• I ,» 


■ ' 3 C V 


I 


3 


10 0 
0 10 
0 0 1 


, -r .i ' ■' 


(l) The matrix 



1TACA TIT No. 998 


81 


is called the "identity matrix" of the third order. The 
matrix 


.1 


2 


1 0 

0 1 

_ 


is the identity matrix of the second order. An identity ma- 
trix of any other order is defined in a similar manner: t 

is a square matrix having 1 for its principal diagonal, ele- 
ments, and 0 for all the other elements. Now, application ot 
the rule of matrix multiplication shows that 


7 4 1 


1 

0 

0 

■ rH 


7 4 1 

5 5 7 
12 3 

10 0 

X : 

‘ 

0 10 
0 0 1 

”333 

£S 

5 5 7 

12 3 
L . -1 

3 3 3 

0 10 

X 

2 0 1 

* 

2 0 1 

0 0 1 


12 3 


1 

H 

W 

W 

1 


It can he seen that any 'square matrix remains unchanged 
upon premultiplication or postmultiplication hy ■the identity 
matrix of the same order. , . ...... 

(2) Multiplication yields 


12 3 

1 

■110 


13 3 

0 2 0 

X 

0 10 

= 

0 2 0 

3 3 1 

— ' T- 


0 0 1 


3 6 1 


In other words, g o st mult ipli cat i on hy the matrix 
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. . "l 1 0 

0 10 
; 0 0 1 

left the first and third columns of the first matrix un- 
changed, and transformed the second column into the sum of - 
the original first and second columns. In a similar way , ■ it 
cap, he seen that 


2 0 0 


1 0.5 0 


2 10 

8 11 

X 

0 10 

= 

8 5 1 

4 3 2 


0 0 1 


4 5 2 


or in other words, po st mul tipli cat i on by the matrix- 


1 0.5 0 
0 10 
0 0 1 


left the first and third columns of the first matrix un- 
changed, and transformed the second column into the sum of 
the original second column and one-half of the original 
first column. Also, l - • 


2 0 1 
0 3 2 

X 

10 0 
0 13 

=3 

2 0 1 
0 3 11 

i 

to 

o 

i — ! 

i 


0 0 1 


1 '0 3 


Hence, po st mult ipli cat i on by the matrix 

10 0 
0 13 
0 0 1 
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left the first and second 'columns of the first matrix urf- . ■ 

changed, and transformed the third column into the,, sum of the 
original third, and three times the original second columns. 
It may he stated, therefore: ; . 

■■ v i • i - ... . , t r : ■ 

Postmultiplication hy a matrix in vfhich all principal 
diagonal elements are unity, the element in-'the . -pth row and 
qth ■ column is ,,k, and all the other elements are zero, is 
equivalent to an addition to the’ iith column of -the -original 
matrix of k times the pth column of the original matrix. 

The effect of several nqnvanishing elements off the 
principal • diagonal can . he understood from the example; 


» 

H 
- 03 
H 

I 


110 


13 5 

0 2 1 

X 

0 12 

e 

0 2 5 

2 0 1 


0 0 1 

i 

2 2 1 


Here, p o st mult ipli'eat i on hy the -matrix . 

•: l l. 6 • 

0 12 

ooi 


left the first column unchanged, - added tb ! .the sechnd column 
the first column, and added to the third c'oluton .the do,u.h.l.e r of- . 
the second column. The knowledge gained from these examples 
ma ^' utilized in the transformation of a set of linear 
equations , ' * £ r$ * ‘i * * ■ [■? 


In order to solve the matr isf -equati on (13) for the un- 
knowns x, y, and z, hoth sides of the equation are first 
postmultiplied hy the matrix - :..-. 


V $ r 



;C 16 ) 
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PostmM.lt ipli cation means' that., , M x is written as., ,a, factpp 
behind * the member ft' of the equation: ■ > ■ 


2 11 


*-*l 

1 -0.5 -0.5 

=' [8, 9, 7] X 

1 -0.5' -0.5 

.1 2 1 


0 . 10 


0 1 0 

1- 1: 2 

. 

0. 0. 1 


0 0 1 
— w. 


( 17 ) 


It was v ,uho.wji i,n equations (.12) that the order of- factors in a 
mat ri x py oid.uqt i s of great importance. The result of a post- 
multiplication differs, in general, from tthose of a ’’premulti- 
plication" in which latter the same factor M x would he 

written he for e the members of equation (1.3):. The rule for the 
construction of M x ie as follows: 

(l) Fill in the principal diagonal; that ia, the diagonal 
of th-e matrix running from top left to bottom 
right, with the digit 1; 


(2) Fill in the remaining places -In* tire’ fi'r'st' row with 
the corresponding values of the original matrix 
divided hy the first element in the first row of 
the original matrix (in the present example 2) 
and hy -1; 


(3) All the remaining eleraputs of M x are to he filled 
in with 0. 


■ I.f the matrix mult i*plicat ( i , qpi I e -qa-rried' out-, 
i'ng result is obtained ‘ * - / 


[*• y, z] X 




5, 3] 




■ 1 1 , * 5 ; ,-0 «,6 

• : ; .(.'. I '!:' 

1 0.5 1 .& ■: 


:> 4 ■ v J 


Oil 


th-e .follow-- 




This equation can be written in an equivalent row matrix form, 
l!f- j's'o desired! * l Of- M 


' X ' ■ ■ 

C(3x + y + z), (l.ly + 0.5 b), (0 ;5y + l,5z)] = [8, 6, 3] ( 19 ) 
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The matrix equat i on . i s . again equivalent to three algebraic 
equati one 

' 2x +' y' +' z = 8 

» ■ 1 . 5y + 0.5z =5 (SO) 

0 . 5y + 1.5z = 3 


It will he recognized that the first of equations (20) is the 
same as the first of equations (l), while the second and the 
third of eouations (20) are identical with equations (3a) and 
(3h), respectively. In other words, postmultiplication hy 

achieved exactly the same transformation of the system of 

equations as did the algebraic operation that originally re- 
sulted in equations (3). 


The matrix -operationB can now be continued by .the appli- 
cation of one more p o st mult ipli e r M s defined as 



( 21 ) 


Po st mult ipli cat i on- of equation (18) is indicated as follows: 


Cat, y. z] X 

2.0 0 

X 

10 0 

= [8, 5, 3] x 

1 0 0 

* 

l’ 1 .5 0.5 


0 1 -0.33 


o- 1 -0.33 

. :. « • 

1.0.5 1.5 


. 0 0 „1 

..'V* 


.60 1 

t 

— - _ 


-• ‘ ” ‘ - 

.. 

r> * - 


When it is carried out, the result is: 


[ y, z] 


x 


2-0 0 
1 1.5 0 
1 0.5 1.333 


= [8, 5, 1.333], 


(33) 
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Equation (23) may "be written in the row-matrix form: 

[(2x + y + z), (l.5y + 0.5z), 333 z)] = [8, 5, 1.33 3] (24) 


This matrix equation is equivalent to the three simultaneous 
algebraic equations • 


2x + y + z = 8 

1 . 5y + 0.5z = 5 ’ (P 5) 

: „ . 1 .333z = 1 .333 

Comparison reveals that the first two of equations (25) are 
identical with the first two of equations (20), Consequently, 
the first of equations (25) is the same. as. the fiyst of equa- 
tion (l), and the second of equations’ (25) is the same as 
equation (3a). The third of equations (25) is identical with 
equat i on ( 5 ) . 

■f ■! t 

Thus, it has been shown that the two datrix multiplica- 
tions Indicated in equations- (17) and (22), respectively, re- 
sult in the same set of linear equations as those obtained by 
the use of methods of high ^school algebra. The characteris- 
tic property of the set of equations (25) is that the last 
equation contains one, the middle one two, and the first one 
three unknowns. Because of this property the equations can 
be easily solved one by one, as was done in equations (6) to 
(8). The square . mat ri x in equation (23) 1 b just a convenient 
representation of the left-hand member of equations (25). It 
falls Into the category known as "triangular (l) 11 matrices be- 
cause all itB elements are zero on one side of the principal 
diagonal . 

It is, therefore, the purpose of the p 0 st mult ipli cat i on s 
to transform the original matrix into a triangular matrix, 
which can be Bolved then i>y elementary a,lgebra without- any 
difficulty. 

The procedure as developed here may seem to ba lengthy. 

A considerable part of the material presented, however, was 
purely explanatory and can be omitted from the actual calcu- 
tions. The essential features were: 

(l) Representation of the equations by the matrix of 

the coefficients of the unknowns 



NACA Ttf V o. 998 


87 


...2X1- • 

• • 121 

. - 112 .. 

and by the row matrix of the right-hand side 

[8 9 7] 



(b) 


(2) Construction of the p o s t mult ipli er 

fl -0.5 -0.5 I 


M 


1 


0 1 0 


( c) 


L° 0 1 i 

according to the rules . „ * : 

I 1 , r f i. - 

(a) Elements of the principal diagonal ; have 1 the ' value . 3T;~'' 

(b) Elements of the first row off the principal,. diagonal 

are equal to the corresponding tel'ek'e'nt s 'o'f the 
original matrix divided by the principal diagonal 
element of the first row and by -1; . ' 


(c) All other 61dments are 0. ■' 

(3) Fostmultlpli.cati on by M L of* both the square matrix 
and the row matrix: ' •••' 


— 

2 11. 


■1 -0.5 -0.5 

= 

• 3 O’’ ■ 0 

12 1 

X 

0 10 

• 1 ' 

1 1.5 0.5 

112 


0 0 1 


1. >0.5 1.5 

— — 


: . * • _ 


_ _ 


[8, 9, 7] x 


1 -0.5 -0.5. 
0 ’ 1 0 


>£8* 5 ,'-3-3 


(e) 


0 


0 


1 
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(4) In all the later operations undertaken, the first 
column and the first row of the right-hand member of equation 
(d) remain unchanged. It is, therefore, permissible to omit 
them from the further considerations and to center the atten- 
tion upon the two-by-two square 11 eubmatri 1 " obtained by 
striking out the first column and the first rows 


1.5 0.5 

O'.' 5 1.5 


(f) 


This submatrix has to be postmultiplied by the matrix M* a 

constructed exactly according to the prescription contained 
i n ( S ) s „ 


M* 

a 


1 -0.33 

0 1 


<«> 


It may be noted that M* s can be obtained from M 3 (see 

equation (2l)) by striking out the first column and the first 
row . 


Postmultiplicat ion results in 


|1 


i 0 . 5 • 


0.5 

1.5 


1 -0.33 


1.5 

0 

0 1 

X 

0.5 

1.333 

_ 

. 

_ 




(h) 


(5) In a similar manner the first element can be omitted 
from the row matrix of the right-hand member of equation (e). 
The submatrix so obtained is 


[5, 3] 

It must be postmultiplied by 


[5, 3] x 


,1 *0.33 


: : [5, 1.333] 


<J> 


(k) 


(6) The complete transformed square matrix consists of 
the first row and the first column of the right-hand side 


r!H 
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member of equation, (d), and of the right-hand side member of 
eauatior (h): • . " 


s ■ ‘0 0 

1 1.5 0 

1 0.5 1.333 



(7) Tire complet e. transformed row matrix consists,' of the 
first element of the: right-hand side member of equation (e), 
and of the’ right-hand side member of equation- (h): 


[ 8/ 5, 1.333] 


Cm) 


(8) The complete transformed matrix equation can be 
written as a combination of (l) and (m) as follows: 


[ x, y , z] x 


2 


0 


0 


11.5 0 


[8, 5, 1.333] (n) 


• .jjL 0..5 1.333J 

This equation is identical with equation (23) and can be 
solved as explained below equation (23). 

The operations outlined in the -preceding eight points 
can be arranged concisely according. to the following scheme: 


8 9.7 

... 2.1 1 

12 1- 
112 

-0.5 ' -6.5 

t. . - - 

1 o" 

0 1 

5 ' • ■■ 3 


1.5 0.5 

. -0.33 

- : 0.. 5 ■ 1.5, 

. 1 

to 

„$0 

• 

rH 

* 

- ' 

1.331 
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The first line in 'the scheme represents 'the original 
row matrix (b). Below it is arranged the original square 
matrix (a). To the right is the p o st mult ipli er (c) except 
that its first column is omitted. This first column is not 
needed since it would enter only in the calculation of the 
first column of the matrix product’ which need not be rewrit* 
ten since it remains unchanged as explained under (4). 

The values under the first full horizontal line are ob- 
tained by p o stmult.iply i ng the row matrix [8, 9, 7] by the in- 
complete-, p o st mult ipli er . This po st mult ipli cat i on is carried 
out exactly according to the rules of matrix multiplication 
given earlier. Because of the incompleteness of the post- 
multiplier only two elements are obtained in the product. 

The first is * 

[8 x (-0.5)] + (9 x 1) + (7 x 0) = 5 

tho second » 

[8 x (-0.5)3 + (9 X 0) + (7 x 1) = .3 

Since the omitted first column of the postmultiplier contained 
the elements 1, 0, 0, multiplication by it would have given 8, 
the unchanged value of the first element of the row matrix. 
This need not be rewritten. 

The' square matrix under the row 5, 3 is computed in a 
similar manner by postmultiplying the original square matrix 
by the incomplete p’o st mult ipl i er . Multiplication of the row 
2, 1 , 1 would give: . 

[S X (-0.5)] + (1 X 1) + (1 X 0) = 0 

[2 X (-0.5)] + (1 X 0) + (1 X l) « 0 

These two zeros need not be listed. Multiplication of the 
second row yields 

[1 x (-0.5)] + (2 X l) + (1 x 0) « 1.5 

[1 x (-0.5)] + (2 x 0 ) + (1 x r) * 0.5 

These two values are list-ed in the first row below the dotted 
line. The elements in the second row are obtained in a simi- 
lar manner : 
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[1 x (-0.5)3 + (1 X i) + (2 X 0) 

[1 x (-0.5)] + (1 X 0) + (3 X l) 

r - 

The column, to the right is the second column of the poet 
multiplier (g). The first "Column of this po st mult ipli er is 
again omitted as unessential for the computations. 

The last two numbers in the scheme are computed by multi 
plication by the incomplete second po stmult ipli er in the same 
manner as described in connection with the preceding opera- 
tions. 

It must be emphasized that the scheme given does not con 
tain a complete presentation of the matrix multiplications. 

It is Just a convenient short-hand reproduction of the compu- 
tations given In more detail under (l) to (8), Familiarity 
with this scheme of operations makes it possible to reduce 
quickly, with the least amount of time wasted In writing, a 
given set of linear equations to tfh*e diagonal form. 

An example Is now giygn., wlthqut explanations, for the 
reduction of a set of f our' Pquafc i on s . 

r 

The equations arei 


4x + 

2y 

+ 

z 

zs 

n 

3x + 

3y 


+ 

2w .= 

16 

X 


+ 

2z + 

w = 

n 


2y 

+ 

ft 

Z + 

i 

y - 

li 


The computations can be arranged as -follows-: 


*•! IS 

0.5 
= 1.5 
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This matrix equation is equivalent to the four algebraic.’" 
equations: ■. ' 7 -” > r? 


4x + 2y + 


*'11 


<Sy- *-• O.Sz ■ \4 3v 


10.5 , 


1-.625 - +'T.5w ■ 10.875 *. 

•. * .1 ' ! r 0 V ' 5 

- 2.385V. = -9.55 ’ 


The last one of the equat i one , give s immediately 

, t r ' * 1 1 ‘ 

, t * 

. ' W = ‘4 

Substitution in the' preceding, equation yields 


Hence ; 4 


1.625b.* 10.875, - 6-=- 4.875 


z = 3. 


Substitution of v and z in the, second equat io-n results in 


Consequently, 


2y =' 10.5' 1.5 = 4 


• y 3 '■ •. ^ 


• ■ r>,‘ . 


ffinally, there is obtained,. 


and thus : >* 


4x a .11. - 3 - .4 = 4 1 


X = 1 •- -■ * '* ' 

: : . * ' . v , , - r * s > 


' 1 J( fi “V *• 


■. : ' v* - ' . . , - v , i i,* ; •■j.r.r./f - : 

' When the --number • of simultaneous linear equat i one t i s . ,- f ■* 

large, the ,cpmputati onal work is greatly .reduced ; by th-e Use 

of a calculating.. machine . A calculating-machine may become 

necessary also for reasons of accuracy, when the computations 

. Z * ‘ - ... ... 
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happen t o . involve, small differences of large number s . An ex- 
ample is. now given for the solution of seven equations' with 
seven unknowns. The equations represent the problem of the 
egg-shape ring with symmetric loads. 

In this example, as in all the matrix computations of 
this report, the. work of writing was further reduced by the 
adoption of the following scheme. Instead of writing out the 
entire matrix multiplier, 'only the first row is put down, and 
even the first element of the first row is omitted. It may 
be remembered that, all the ' element 6 omitted are either equal 
to unity, or to zero. 

The first row >in table 17 contains (-1) time’s the exter- 
nal loads of the egg-shape ring with symmetric loads. The 
following seven rows, under the dotted line, are identical 
with the> rows of the operations table (table 6). The next^ 
row, under the solid.' line’ and included in par entheseV, ‘ pre- ' 
sents the multipliers. They are (-1) times the ratios ob- 
tained by dividing the elements of. the first row under the 
dotted line, except the first element of this row, by the 
first element of the first row. The first ratio, 

-(8.92216)/ (-3.34833) = 2.66466, is listed in the second 
column (in the column of 8.92216). The other ratios follow 
in the third, fourth, and so forth, columns. 

- , * * i ■ 

The row .und.q^ the- ratios "contains the transformed exter- 
nal loads. The value in the ith column of this row 1 b ob- 
tained by multiplying the , first .element, in the first row of 
the table by the multiplier in the ith column and adding to 
the product algebraically the value of the element in the. 
first row and the ith column. Thus, for instance, in the 
second column the value 133,233 was calculated as follows: 

2.66466 X 50 + 0 = 133/293.' . • •' 

It may be seen that in this manner /only,, six transformed ex- 
ternal load elements are "obtained from the original seven 
external load elements. 

Following the same rule, ,th.e ..original seven-by- seven 
matrix is transformed into the six-by-six matrix contained 
between the second dotted line and the second full ope.s., For 
instance., the element i nf it he 'fourth' '"col umn k'nld t !h e,' .thirteenth 
row of. the entire table was". "d'omput'ecT'as-' ‘foil owe': . 

■ 3/96771 X T;i8'4'5B : 30. 956%' = -86.3549 “ 
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This scheme is followed until the external load elements 
are reduced to a single element, and the matrix is also re- 
duced to a single element. The leading columns represent then 
the equations of the triangular system'. For instance, the 
leading column of the third step may he written as 

-1.77854 v B + 0.204688 u B - 8.96572 w c + 0.77576? v 0 = -38.3834 
The leading column of the sixth step ist 

-1.034845 v 0 - 0.654970 u 0 = -*19,2443 

The full set of these equations can he solved as shown 
earlier . 

As a second example of the procedure actually used in 
the calculations, table 18 is presented. It contains the 
transformation of the matrix of the circular ring with anti- 
symmetric loads. 

More information on matrices, written from the standpoint 
of the person Interested in applications of the matrix calcu- 
lus to problems in physics and engineering, may he found in 
"Elementary Matrices and Some Applications to Dynamics and 
Differential Equations 11 "by &. A. Frazer, VT. J. Duncan, and 
A. R. Collar (University Press, 1938). The procedure used in 
the present report can he shortened further if advantage is 
taken of the fact that all the matrices discussed in this re- 
port can he set up in a symmetric form. Such a simplified 
procedure was suggested by M. H, Doolittle in "Method Employed 
in the Solution of Normal Equations and the Adjustment of a 
Triangulat 1 on , " U. S. Coast and Geodetic Survey Heport, 1878, 
pp . 115-120. A discussion of the Doolittle technique may he 
found in an article by Paul S. Dwyer in the December 1941 
Issue of the Annals of Mathematical Statistics, Vol. XII, No. 
4, pp . 449-458. A more modern approach to .the solution of 
simultaneous equations arising in engineering is given in the 
paper "A Short Method for Evaluating Determinants and Solving 
Systems of Linear Equations with Real or Complex Coefficients" 
by Prescott D. Grout Transactions A.I.E.E,, vol. 60, 1941 

P. 1235. ... _ 



to 

CD 


fable 6. Operations gable 


Operation 

e a 

■h . 

. \ ' 

V 

% 

% 

J ~ 

i- *■ 

T c • • ' 

EjJy a _ 10-4 in. 

-3.34833 

8.92216 

-2.69614 

3.96771 



! 

.. ■ . u 

[2]wg=10“ 4 rad 

8.92216 

-387.866 

11.4697 

-13.1014 

-51.242 

8.10267 

0. 

[3]v b « 10“ 4 in. 

-2.69614 

11.4697 -4.00991 

3.4352 

-8.1026 1 ! 

0.66158 

- 

Wug^lO" 4 in. 

3.96771 

-13.1014 3.4352 

-30.9566 

0 . 

0 • 

26.2668 

» r- ■ 

[sjw^wlO" 4 rad 


-61.242 

'-8 .1026' 

r o 

-288.367 

-2.95522 

* 

-5 ^24667; 

< i 

[6]v c »10“ 4 in. 


8.10261 

r 0.66158 

0 

-2.95622 

-1.90205~ 

-0.88929 

. n 

[?] o c =10“ 4 in. 


0 

0 

26.2058 

-5.24667 

- .88929, 

-27.0833 


• 1 ’ ' V 


The notation adopted for the displacements is as follows: u 

stands for a tangential, v for a radial displacement, w for a rotation, fhe 
subscript refers to .the point that moves? 
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Table 10. Operations Table 


Operation 

e a 

e b 


t bc 


R C 

v^=10“ 4 in. 

-8.3726365 

31.372640 

-7.6032219 

10.097316 



i» B »10~ 4 rad. 

31.372640 

-479.93253 

33.93750 

-40.987910 

-81.64964 

10 

Vg-10~ 4 in. . 

-7.6032219 

33.93750 

-9.1891329 

9.4690625 

-10 

0.8164964 

n B =n^=10" 4: in. 

10.097316 - 

-40.987910 

9.4690625 

-13.515682 

-8.855104 

-1.3664455 

\Up»10~ 4 rad. 

- 

•81.64964 

-10 

-8.855104 

-354.21850 

-5.791741 

v c =10~ 4 in. 


10 

0.8164964 

-1.3664455 

-5.791741 

-fc. 6516009 
. 



co 

-j 
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Table 11 » Operations table for circular ring 


Operation 


|[ll w A “10“ 3 rad 
[2} u A “10" S i:a. 

[3] v< B “10‘ 3 rad 

[ 4 ] ▼ ' -10-3in, 

B 

fs] u B ,10“ 3 in. 

[5] w c =“10’ s rad 
[7 3 u c *10" z in. 



N, 

in , 4 lb . 

lb. 

N 

B 

in. -lb. 

E 

B 

lb. 

t b 

lb. 

N 

C 

in. -lb. 

T 

C 

lb. 

• 

-281.95 

-49.079 

-29.966 

-4.733 

64.675 




-49.079 

-52.296 

64.676 

-22.441 

51.516 



« 

-29.966 

64.675 

-439.849 

31.443 

-50.642 

56.5117 

6.632 


-4.733 

-22.441 

31.443 

-12.338 

20,14 

"8.84a 

0.624 


64.675 

51.516 

-60.642: 

20.14 

-52.618 

6.632 

0.0685 

• 



56.5117 

8.842 

6.632 

-157.899 

-1.563 




6.632 

0.524 



0.0685 

-1.563 

-0.32a 
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Table -12. Relaxation Table for Circular Ring 


Operation 

N 

A 

in. lb. 

T 

A 

lb. 

V 

in. -lb. 

lb. 

t b 

lb. 

N c 

in. -lb. 

lb. 

External 

-1.835 

-8.747 

-64.970 

59.476 

38.104 

-53.135 

-23.860 

loads 








0.42^63 



23.76 

3.71 

2.78 

-66.4 

-0.665 


-1.84 

-8.75 

-31.22 

63.19 

40.88 

-119.5 

-24.5 

-76.2[7} 



-505 

-39.9 

-5.21 

119.2 

24.6 


-1.84 

-8.75 

-536.2 

23.3 

35.67 


0.1 

-2,o[s3 

60.0 

-129.4 

880 

-62.8 

101.2: 


-13.26 


58.2 

-138.16 

343.8 


136.87 

-113.3 

-13.2 

l[ll 

-282 

-49.1 

-30.0 

Mam 

64.7 




-223.8 

-187.3 

313.8 

-44.2 


-113.3 

-13.2 

-6003 



-398.0 

-31.4 

WBSM 

93.9 

19.3 


-223.8 

■HUM 

-84.2 

BS 

mmm 

-19.4 

6.1 

-0.35[6] 


■m 

-19.8 

X£l 

KSifl 

55.4 

0.55 


-223.8 

-187.3 

-104.0 

-78.7 

195.2 

36.0 

6.65 

-0.90[33 

27.0 

-58.3 

396 

-28.3 

45.5 

-50.9 

-5.96 


-196.8 

mmm 

292 

-107.0 

240.7 

-14.9 

0.69 

O.lCQ 

-28.2 

■m 

-3.0 

-0.5 

6.5 




-225.0 

-250.5 

289 

-107.5 

247.2 

-14.9 

0.69 

-o.iCel 



-5.6 

-0.9 

-0.7 

15.8 

0.16 


-225.0 

-250.5 

283 

* 1 

246.5 

0.9 

0.87 

-0.06[33 

1.8 

-3.9 

28.4 

■Efifl 

3.04 

-3.24 

-0.04 


-223.2 


BW 

mmm 

249.54 

-2.34 

0.83 

0.05 Cl] 

-14.1 

mmm 

loom 

■EEa 

3.23 




-237.3 

-256.9 

309.9 

-110.5 

252.7 ? 

-2.34 

0.83 

-4.87[ 2^ 

239 

255 

-316 

109.3 

-251 




1.6 

-1.9 

-6.1 

-1.2 

1.77 

-2.34 

0.83 

Repults of 








a check 








table 

1.637 

-1.944 

-7.246 

-1.37 

1.969 

-2.792 

0.412 
































































Table 13. Operations Table 


Operation 

h a 

In. -lb. 

T i 

lb. 

in, -lb. 

[ lj w A "lO**®rad. 

“281.95 

-49.079 

-29.966 

[2] v -10"3in. 

A 

-49.079 

-52.296 

64,676 

[ 3 ] w *10" 3 rad. 
B 

-29.966 

64.675 

-453.287 

[ 4 ^ Tg B 10" a in. 

“4.7334 

-22.441 

31.2348 

[5j u^lO^in. 

64.675 

51.516 

-50,1390 

[6] v» c “10' 3 rad. 



49.8588 

[7] ylO'3 in. 



7.134970 


Braced Ciroular Rin| 


B' 

lb. 

-4.7354 

-22.441 

31.2348 

-29,6649 

12.9752 

8.63406 

7.709565 


lb. i 


64.675 

5 : 
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Table 14. Relaxation Table for Braoed Circular Ring 


Operation 

n a 

in. -lb. 

t a 

lb. 

n b 

in . -lb • 

R 

B 

lb. 

lb. 

K c 

in. -lb. 

T 

C 

lb. 

External 

loads 

-0.4t63 

-1.835 

-8.747 

-54.970 

-19.92 

69.476 

-3.46 

38.104 4 
-2.85 

-53.135 

68.6 

-23.860 

0.424 

' 

~8.3t.73 

-1.835 

-8.75 

-74.89 

-59.15 

56.026 

-64.05 

35.254 

-25.1 

15.46 

8.80 

-23.436 

27 . 6 

~0.375t3l 

-1.835 

11,25 

-8.75 

-24.25 

-134.04 

170.0 

-8.02 

-11.70 

10.1 

18.8 

24.26 

-18.70 

4.2 

-2.67 

0.l[l3 

9.41 

-28.2 

1 

36.0 

-3.0 

-19.72 

0.47 

Si 

5.56 

n 

-0.24t4l 

-18.79 

1.14 

mm 

KEfli 

IBM 

35.4 

-3.1 


1.5 

-1.85 

-0.05U3 

-17,65 

1.50 

-32.5 

-3.24 

25.5 

22.6 

-12.12 

-1.56 

32.3 

2.5 

3.49 

-2.49 

-0.35 

-0.356 

0.064L13 

-16.15 

-18.05 

-35.74 

-3.14 

48.1 

-1.92 

-13.68 

-0.3 

34.8 

4.15 

1.00 

-0.706 

0.1L43 

-34.2 

-0.47 

-38.88 

-2.24 

46.2 

3.1 

mm 

mnmm 

38.95 

1.30 

1.00 

0.86 

-0. 706 
0.771 

0.0ltl3 

-34.67 

-2.82 

-41.12 

-0.49 

49.3 

-0.3 

-16.95 

-0.05 

40.25 

0.65 

1.86 

0.065 

-0.795^23 

-37.49 

39,1 

-41.61 

41.6 

49.0 

-51.5 

-17.00 

17.8 

40.90 

-41.0 

1.86 

0.065 


mi 

0 

wmt 

0.8 

-0.1 

1*86 

0.065 

Results of 
a oheok 
table 

1.5217 

-.05656 


-.072270 

0.01066 

1,0078 

0.00159 



















































Operation 


w^-lO^rad. 

v = 10" 3 in. 
A 

Wg=10~^rad. 

Vg^lO-^in. 

u *=10“ 3 in. 

B 

w *10 -3 rad. 
C 

v *10” 3 in. 

C 

u ■lo'^in. 
CD 

w =*10 -S rad. 
D 

v^lO^in. 


Table 16. Operations Table for Toraion of Bgg-Shaped Ring 


h 

in. -lb. 

t a 

lb. 

"b 

in. lb. 

*5 

lb. 

t b 

lb. 

*C 

in. -lb. 

R 

C 

lb. 

T 

. CD 

lb. 

»0 

in. -lb. 

lb. 

-523.185 

-217.240 

178.054 

-82.8465 

234.496 






-217.240 

-175.295 

234.496 

-79.8920 

174.431 






178.054 

234.496 

-1606.64 

-291.974 

-1496.08 

367.170 

-278.790 

1314.66 



-82.8465 

-79.8920 

-291.974 

-248.18 

-660.021 

278.790 

-185.110 

744.755 



234.49- 

174.431 

-1496.08 

-660.021 

-2955.38 

1314.66 

-744.755 

2778.29 




• 

367.17 

278.790 

1314.66 

-1246.76 

403.048 

-1278.84 

-81.3496 

10 



-278.790 

-186.110 

-744.755 

403.048 

-211.246 

738.001 

-10 

0.816490 



1314.66 

744.755 

2778.29 

-1278.84 

738.001 

-2781.37 

-8.85510 

-1.36345 






-81.6496 

10 

-10 

0.816496 

-8.85510 

-1.36645 

-354.218 

-5.79174 

-6.79174 

-2.6516 


10 


0.816496 


-1.36645 


•5.79174 -2.6516 






Table 16. Operations Table for 


Operation 

k a 


h b 

"b 

T 

B 


in. “lb. 

lb. 

in. -lb. 

lb. 

lb. 

w =10" 3 rad. 
A 

-523.185 

*217.240 

178.054 

-82.8465 

234.496 

u “lO^in. 
A 

-217.240 

-175.295 

234.496 

-79.8920 

174.431 

■w B -» 10 “ 3 rad. 

178.054 

234.496 

-1614.32 

*292.08 

-1495.94 


-82.8405 

*79.8920 

-292.08 

-255.880 

-665.458 

Ug'lO^in. 

234.496 

174.431 

-1495.94 

-665.458 

-2959.23 

w c “10“ 3 rad* 



367-17 

278.790 

1314.66 

v -10" 3 in. 
C 



-278.790 

-185.110 

-744.765 

« c =10" 3 in. 

w “lO'^adj 
D 

v ='10"3 in* 
D 



1314.66 

744.765 

2778.29 


I 


1 of Braced Egg-Shaped Ring 


“c 

ill #,-lb • 

\ 

lb. 

T 

CD 

lb. 

n d 

in. -lb. 

R 

D 

lb. 

367.17 

-278.79 

1314.66 



278.790 

-186.110 

744.756 



1314.66 

-744.765 

2778.29 



-1246.76 

403.048 

-1278.84 

-81.6496 

10 

403.048 

-211.246 

738.001 

-10 

0.816496 

II 

-1278.84 

738.001 

-2781.37 

-8.86510 

—1.36645 

-81.6496 

-10 

-8.85510 

-364.218 

-5.79174 

10 

1 

0.816496 

-1.36645 

-5.79174 

-2.6516 
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Table 17. Matrix Solution of Egg-Shaped Ring With Symmetric Loads 





























rHI T 


NACA TN No. 998 


105* 


Table 18. Matrix Solution for Circular Ring With Antisymmetric Loads 


\ 

.835_d7_ 

281.950" 

-49.079 

-29.966 

-4.7334 

64.675 


8. 74716 
-49.079' 
-52.296 

64.675 
-22.441 
51.516 


% 

54.9703 _ 
”-59.966 

64.675 
-453.287 
31.2348 
-50.1390 




- 59.4262 

”-4.7334 

-22.441 

31.2348 

-29.6649 

12.9762 


-38^1038 
" 64.6^5 
51.516 
-50.1390 
12.9752 
-65.6152 


(-.17407 
8.4277 
-&77W2W - 
69.89118 
-21.61706 
40.25804 


-.10628 

64.77527 

6§~.iT9TlT " 

-450.1022 

31.7379 

-57.0127 


-.016788 

-59_.5070_ 

~-Sl7 61706 

31.7379 

-29.5854 

11.88943 


,2293846) 

-37.6829 

"4S.-55Fo~ 

-57.01275 

11.88943 

-40.7797 


(1.69744 

68.23776 

-3W8,4W7 

-2.7934 

7.29590 


-.694073 
-63.6709 
= 2»7§3C ” 
-18.9050 
-8.00098 


0.920124) 

-29.92837 

~7.2'S5S0 

-8.00095 

-3.73731 


(-.00826334 
-64.23409_ 
-18". 88X95 
-8.06119 


.02155606) 

-28_.4574 

-e.osirs 

-3.58004 


(-0.4269241) 
-1^03431 _ 
-.13852 
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Figs. 1,2*3 



I00 # 


FIG. I. CURVED BEAM. 



FIG. 2. ARC AB. 


FIG. 3. ARC AC. 
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Figs. 6,7 



C FORCES ACT UPON BEAM) (FORCES' ACT UPON CONSTRAINTS) 

FIG. 6 SIGN CONVENTIONS FOR FORCES AND MOMENTS. 
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Fig. 8 



FIG. 8 DEFLECTED SHAPE OF EGG-SHAPED RING 

WITH SYMMETRIC LOADS. 



MOMENT X I0" 3 (IN. LB.) 


9 



FIG. 9. BENDING MOMENT DISTRIBUTION OF EGG-SHAPED RING WITH 

SYMMETRIC LOADS. 


31 

cT 
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Figs. 10, II 





NOTATION FOR ANALYTIC SOLUTION OF 

RING PROBLEM. 
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Fig. 12 



FIG. 12. 

BENDING MOMENT DISTRIBUTION IN CIRCULAR RING. 


Figs. 13,14 
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BRACED CIRCULAR RING WITH ANTISYMMETRIC LOADS. 



(a) RESOLUTION OF TANGENTIAL DISPLACEMENT. 

(b) BEAM CONVENTION FOR STRAIGHT BAR. 

(C) RESOLUTION OF RADIAL FORCE. 

Cd) RESOLUTION OF TANGENTIAL FORCE. 

FIG. 14. 

RESOLUTION OF DISPLACEMENTS AND FORCES. 





FIG. 15. RESOLUTION OF FORCES AT FIXED END. 
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Fig. 16 



FIG. 16. 

BENDING MOMENT DIAGRAM OF BRACED CIRCULAR RING. 
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